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PREFACE 


This  report  was  prepared  by  the  Douglas  Aircraft  Company,  Long  Beach, 
California,  under  Contract  DNA  001 -75-C-0216  and  documents  the  overall 
program  descriptions  and  method  of  analysis,  the  input  and  output  data 
descriptions,  the  program  operation  and  a  sample  problem.  This  work  was 
performed  under  Program  Element  NWE  D  62704H,  Project  N99QAXA,  Task  Area 
E500,  Work  Unit  04  and  was  funded  by  the  Defense  Nuclear  Agency  under: 

RDT  &  E  RMSS  Code  B342075464N99QAXAE50004H2590D.  Funding  of  this  effort 
was  also  supported  by  the  Air  Force  Weapons  Laboratory  under:  Program 
Element  62601F,  Project  8809,  Task  03,  Work  Unit  40.  Inclusive  dates  of 
research  and  development  were  May  1975  through  June  1976. 

This  work  was  also  performed  under  Program  Element  NWE  D  62704H, 

Project  N99QAXA,  Task  Area  E502,  Work  Unit  01  and  was  funded  by  the 
Defense  Nuclear  Agency  under:  RDT&E  RMSS  Code  B34207T462N99QAXAE50201H2590D. 
Inclusive  dates  of  research  and  development  were  August  1976  through 
August  1977. 

Volume  II  of  this  report  details  the  unsteady  aerodynamic  procedure 
and  Volume  III  contains  the  Fortran  listing  of  the  program. 

J.  A.  McGrew  was  the  program  technical  director  for  this  task. 

The  technical  development  was  performed  by  J.  P.  Giesing  and  T.P.  Kalman 
with  the  assistance  of  Dr.  W.  P.  Rodden.  The  programming  effort  was 
carried  out  by  T.  P.  Kalman  and  H.  H.  Croxen. 
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SECTION  I 


INTRODUCTION 

The  VIBRA-6  computer  program  is  similar  to  earlier  versions  of 
VIBRA  (Vehicle  Ine^stic  Bending  Response  Analysis)  which  were  programs 
designed  to  calculate  the  structural  response  of  aircraft  exposed  to  a 
nuclear  blast.  Specifically,  the  present  program  is  an  extension  of 
VIBRA-4  (Ref.  1)  for  response  analysis  of  arbitrary  wing-body  configura¬ 
tions  at  subsonic  speeds.  Aerodynamic  interference  is  accounted  for 
among  nonplanar  lifting  surfaces;  e.g.,  wing,  stabilizer,  canard  and  fin, 
and  among  slender  lifting  bodies;  e.g.,  fuselage,  nacelle,  and  external 
stores.  At  present,  VIBRA-6  has  no  capability  for  response  analysis  at 
supersonic  speeds,  and  VIBRA-4  continues  to  provide  for  that  requirement. 

The  primary  extension  is  in  the  area  of  subsonic  aerodynamic  loads 
with  the  moving  blast  wave  considered  as  a  travelling  gust.  The  aerody¬ 
namic  loads  for  wing-body  interference  are  based  on  the  extensions  of 
the  Doublet-Lattice  Method  for  nonplanar  lifting  surfaces  to  include 
slender  body  theory  and  the  Method  of  Images  to  account  for  interference 
(Refs .  2  and  3) . 


1.  Hobbs, N. P.  ,Zartarian,G. ,  and  Walsh,  J.P.,A  Digital  Compute r  Program 
foi  Calculating  the  Blast  Response.  of  Aircraft  to  Nuclear  Expiation* ,  Mr 
Force  Weapons  Laboratory,  Report  No.  AFWL-TR-70-140,  Vol.I,  April  1971. 

2.  Giesing,  J . P. , Kalman ,T. P. ,  and  Rodden,W.P.,  Subionic  Unsteady  A erody- 
namics  for  General  Cun figure tions ;  Part  II -Application  oft  the  doublet- 
lattice  Method  and  the :  Method  of,  Images  to  Li  f  ting -Surface/Body  Inter¬ 
ference,  Air  Force  Flight  Dynamics  Laboratory,  Report  No.  AFFDL-TR-71 -5, 
Part  II,  April  1972. 

3.  Giesing,  J.P., Kalman,  T.P.,  and  Rodden,  W.P.,  "Subsonic  Steady  and 
Oscillatory  Aerodynamics  for  Multiple.  Interfering  Wings  and  Bodies," 

J.  Aircraft,  Vol .  9,  No.  10,  pp.  693-702,  1972. 


Volume  I  of  this  report  contains  the  overall  technical  and  program 
description.  Volume  II  of  this  report  documents  the  modifications  to 
Reference  2  to  add  the  travelling  gust  field,  improvements  in  the  aerody 
namic  influence  coefficients,  and  changes  in  the  aerodynamic  load  output 
Volume  II  also  contains  discussions  of  the  aerodynamic  modelling  of  an 
entire  aircraft  with  a  simple  example  configuration  as  well  as  the  de¬ 
tails  of  the  aerodynamic  program  subroutines. 

The  solution  method  for  the  transient  response  has  been  changed 
from  the  method  in  VIBRA-4  of  numerically  integrating  the  equations  to 
an  inverse  Fourier  transform  method.  The  early  versions  of  VIBRA  con¬ 
sidered  structural  inelasticity  and  large  disturbances,  but  VIBRA-4 
restricted  the  structure  to  elastic  deformations  while  still  considering 
large  disturbances  in  the  rigid  body  response  motion.  Experience  with 
VIBRA-4  has  indicated  that  particularly  for  large  bomber  or  tanker 
aircraft,  if  primary  structural  failure  occurs,  it  generally  occurs 
before  substantial  changes  in  the  Eulerian  angles.  The  further  restric¬ 
tion  to  small  disturbances  permits  VI8RA-6  to  utilize  linear  response 
analysis  techniques  that  are  also  consistent  with  the  assumptions  of 
the  linearized  aerodynamic  analysis.  The  theoretical  aerodynamic  loads 
are  not  known  for  arbitrary  transient  motions  but  are  only  known  for 
harmonic  motions;  i.e.,  in  the  frequency  domain.  The  appropriate  linear 
response  analysis  method  is  therefore  a  Fourier  transform  method,  and 
Zartarian  (Ref.  4)  has  demonstrated  the  feasibility  of  calculating 
transient  blast  loads  by  the  Doublet-Lattice  Method  and  Fourier  trans- 

4.  Zartarian,  G.,  Application  ofa  the  Vou.bte.t- Lattice  MeXkod  ^on i  Detci- 
mination  o ft  Blast  Load- s  on  Lifting  Suifiaccs  at  Subsonic  Speeds , Air  Force 
Weapons  Laboratory,  Report  No.  AFWL-TR-72-207 ,  January  1973. 


forms.  The  frequency  response  of  the  vehicle  in  a  harmonic  travelling 
gust  field  becomes  the  fundamental  part  of  the  solution.  From  the 
time  history  and  orientation  of  the  blast  overpressure  and  its  follow¬ 
ing  travelling  gust,  the  Fourier  transform  of  the  gust  field  may  be 
calculated.  The  product  of  the  vehicle  frequency  response  and  the 
transform  of  the  gust  field  is  the  Fourier  transform  of  the  transient 
response  of  the  vehicle  to  the  blast.  Its  inverse  transform  is  the 
desired  transient  response. 

The  assumed  linearity  of  the  system  permits  superposition  of  the 
blast  response  loads  with  the  loads  in  trimmed  flight.  The  trimmed 
maneuvering  conditions  considered  are  level  flight  or  a  symmetrical 
pull-up  or  pushover  with  the  velocity  vector  horizontal  at  the  time  of 
blast  wave  intercept  (as  in  VIBRA-4),  and  a  level,  climbing,  or  de¬ 
scending  steady  turn  (VIBRA-4  considered  the  steady  turn  at  constant 
altitude).  A  new  trim  solution  is  provided  that  is  based  on  the  aerody¬ 
namic  influence  coefficients  from  the  wing-body  theory.  Angle  of  at¬ 
tack  and  elevator  setting  are  determined  for  a  given  speed,  load  factor, 
elevator,  rudder,  and  aileron  settings  are  determined  for  a  given  speed, 
load  factor,  rate  of  climb,  and  altitude  in  the  steady  turning  meneuvers. 
The  final  equilibrium  loads  on  the  aircraft  are  then  the  sum  of  the 
elastic  aircraft  trimmed  flight  loads  and  the  incremental  loads  due  tu 
the  dynamic  perturbation  response. 

The  interaction  between  the  response  caused  by  the  blast  loading 
and  the  autopilot  can  be  a  critical  problem  and  has  not  been  considered 
in  previous  versions  of  VIBRA.  This  interaction  can  cause  significant 
stresses  as  the  autopilot  attempts  to  restore  the  disturbed  vehicle  to 
its  initial  trimmed  flight  condition.  The  autopilot  is  included  in  the 
system  by  adding  the  transfer  functions  of  the  autopilot  components. 
These  equations  relate  the  signals  sensed  from  the  vehicle  motion;  e.g., 
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by  rate  gyros  or  accelerometers ,  to  the  motions  of  control  surfaces  com¬ 
manded  by  the  active  control  system.  The  additional  equations  are  writ¬ 
ten  as  transfer  functions  which  are  expressed  as  ratios  of  polynomials 
in  the  Laplace  transform  variable  and  become  ratios  of  polynomials  in 
the  frequency  domain  which  then  may  be  used  to  augment  the  aeroelastic 
equations  of  motion. 

The  internal  structural  loads  (called  integrated  loads)  are  found 
by  integrating  the  resultants  between  the  applied  aerodynamic  loads 
(from  the  gust  and  the  induced  motion)  and  the  inertial  reactions  (cal¬ 
led  equilibrium  loads).  The  structural  loads  are  shears  and  moments 
(bending  moments  and  torques)  and  are  obtained  by  numerical  integration 
of  the  pressure  loads  at  the  aerodynamic  panel  and  body  points  and  the 
inertial  loads  at  the  structural  mass  points.  Internal  loads  can  be 
calculated  for  any  or  all  of  the  components  of  the  vehicle;  i.e.,  the 
wing,  stabilizer,  fin,  fuselage,  nacelles,  pylons,  etc.,  and  at  as 
many  stations  as  desired.  The  internal  loads  may  also  be  converted  into 
stresses  by  providing  a  stress  transformation  matrix  containing  the 
necessary  section  property  data  at  each  load  station. 

An  option  has  also  been  included  for  experimental  correlation 
studies.  The  measured  blast  characteristics  can  be  input  and  the  time 
histories  of  pressures  at  the  transducer  locations  can  be  calculated 
and  compared  with  the  measurements.  The  assumption  is  made  that  the 
model  is  rigid. 

The  concept  of  dynamic  core  has  been  utilized  throughout,  thus 
restricting  the  limitations  of  problem  size  only  to  that  of  computer 
capacity. 

Volume  III  of  this  report  contains  the  program  code  listings. 
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SECTION  II 


GENERAL  THEORY 

The  VIBRA-6  computer  program  is  a  modularized  extension  of  the  VIBRA- 
4  program.  It  consists  of  twelve  principal  modules.  Output  data  from 
four  major  modules,  as  explained  below,  may  be  saved  on  tape,  at  the 
option  of  the  user,  so  that  subsequent  runs  will  not  require  recalculation 
of  the  same  data.  The  twelve  module  interactions  and  the  saved  data  are 
shown  in  Figure  1.  Additional  modules  for  further  extensions  of  capability 
may  be  inserted  with  relative  ease  in  computer  programs  organized  in  this 
fashion,  since  each  module  is  a  stand-alone  subprogram  with  the  inter¬ 
facing  provided  by  the  control  module. 

The  VIBRA-6  program  has  been  coded  using  the  principles  for  a  dynamic 
core  allocation.  The  program  has  been  coded  so  that  the  computer  core 
required  for  execution  of  any  specific  solution  consists  only  of  that 
necessary  for  the  program  (using  an  efficient  overlay)  and  the  core 
required  for  the  data  associated  with  the  largest  module  to  be  executed 
in  the  case.  This  results  in  a  cost  savings  since  the  user  is  charged 
only  for  the  core  actually  used. 

Although  Figure  1  serves  to  illustrate  the  efficiency  of  the  modular 
program  under  the  direction  of  the  control  module,  a  flow  chart  is  neces¬ 
sary  to  show  the  logical  flow  of  the  computational  sequence.  The  flow 
chart  is  shown  in  Figure  2.  The  flow  chart  shows  the  necessary  steps  In 
the  transient  response  analysis  by  modal  and  Fourier  transform  methods. 

The  steps  in  order  Include  calculation  of: 
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Figure  2.  Program  Analysis  Flow  Chart 
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1.  The  generalized  aerodynamic  forces  for  motion  and  the  harmonic 
travelling  gust  field  and  the  physical  aerodynamic  forces  for 
load  calculations. 

2.  Generalized  mass,  stiffness  and  damping  forces  for  the  structure. 

3.  The  Integrated  loads  (i.e.,  the  internal  shears,  bending  moments 
and  torques)  caused  by  unit  aerodynamic  and  Inertial  control 
point  loads. 

4.  The  transfer  functions  and  modal  data  necessary  to  obtain  addi¬ 
tional  terms  in  the  modal  equations  of  motion  arising  from  the 
active  control  systems. 

5.  The  frequency  response  of  the  vehicle  modes  to  the  unit  harmonic 
travelling  gusts  of  varying  orientations. 

6.  The  symmetrical  and  anti  symmetrical  integrated  loads  for  the  unit 
gust. 

7.  The  trim  load  distribution  for  level  flight,  a  symmetrical  pull- 
up,  or  a  climbing  (or  descending  or  level)  steady  turn. 

8.  The  time  history  of  the  blast  induced  gust  field  velocities. 

9.  The  Fourier  transform  of  the  blast  time  history,  the  Inverse 
Fourier  transform  of  the  frequency  response  of  the  vehicle  to 
the  blast  to  obtain  the  transient  time  history  of  the  response, 
and,  finally,  the  Iteration  for  the  critical  range. 

The  general  theory  for  each  of  the  computational  steps  is  outlined 
below.  The  detailed  equations  are  given  later  with  the  descriptions  of 
the  various  modules. 
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The  subsonic  aerodynamic  forces  are  calculated  by  the  Doublet-Lattice 
Method  and  the  Method  of  Images  (Refs.  2  and  3).  The  aerodynamic  influence 
coefficients  are  determined  for  harmonic  motion,  and  the  load  distribution 
is  found  for  the  harmonic  travelling  gust  field.  Then  using  the  aerodyna¬ 
mic  mode  shapes  which  are  read  in  for  each  rigid  body  and  vibration  mode, 
or  determined  from  inertial  modes,  internally,  if  the  h-a  modal  input  is 
used,  and  specified  at  the  aerodynamic  control  points,  the  generalized 
aerodynamic  forces  are  obtained  for  motion  and  for  the  gust  field.  The 
generalized  forces  correspond  to  unit  generalized  coordinates  and  a  unit 
gust  amplitude.  These  aerodynamic  matrices  are  formed  for  a  limited  number 
of  reduced  frequencies  and  are  dependent  on  the  Mach  number  but  are  inde¬ 
pendent  of  the  flight  velocity.  These  data  are  saved  on  tape. 

The  generalized  mass  and  stiffness  matrices  are  calculated  from  the 
mass,  frequency  and  inertial  mode  shape  data.  The  mass  matrix  and  the 
vibration  frequencies  are  read  in  along  with  the  rigid  body  and  vibration 
mode  shapes  as  specified  at  the  mass  points.  The  generalized  mass  matrix 
is  found  from  the  mass  matrix  and  the  modes;  and  the  generalized  stiffness 
is  determined  by  the  generalized  mass  and  the  vibration  frequencies. 

The  Integrated  loads  consist  of  the  internal  shears,  bending  moments, 
and  torques  at  stations  and  in  directions  specified  by  the  user.  From  the 
geometry  of  the  structural  model  and  the  location  of  the  mass  points  and 
aerodynamic  panel  and  body  points,  the  structural  reactions  are  determined 
for  unit  values  of  the  Inertial  and  thrust  forces  and  the  aerodynamic 
forces,  and  unit  generalized  response.  These  Integrated  inertial  loads 
and  the  Integrated  aerodynamic  loads  are  saved  separately  on  tape  until 
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the  Inertial  and  aerodynamic  load  distributions  are  found  In  the  response 
analysis  from  which  the  combined  Integrated  Inertial  and  aerodynamic  loads 
are  calculated  for  use  In  the  stress  analysis.  Loads  due  to  engine  thrust 
are  Included  also. 

The  active  control  system  transfer  functions  are  found  and  the  kine¬ 
matic  relationships  between  sensed  and  commanded  motion  established  with 
Input  data  defining  the  sensed  and  driven  degrees  of  freedom.  These  data 
are  saved  in  core  for  subsequent  use  in  the  frequency  response  and  unit 
gust  load  analyses. 

At  this  point,  all  data  have  been  Independent  of  flight  condition, 
with  the  exception  of  Mach  Number. 

The  modal  frequency  response  of  the  aircraft  Is  found  from  the 
generalized  mass,  aerodynamic,  and  stiffness  matrices  and  the  active 
control  system  for  unit  gusts  of  all  available  orientations  at  a  speci¬ 
fied  flight  altitude  and  velocity.  Symmetric  and  antisymmetric  solutions 
are  formed.  These  solutions  are  obtained  at  arbitrary  user  specified 
frequencies,  with  spline  interpolation  employed  to  obtain  the  aerodynamic 
forces  at  frequencies  other  than  those  used  for  the  basic  aerodynamic 
matrix  calculations.  These  data  are  saved  on  tape. 

The  symmetric  and  antisymmetric  Integrated  loads  due  the  unit  gust 
are  found  from  the  modal  frequency  response  and  the  unit  load  solutions  for 
all  orientations  and  saved  on  tape. 

The  availability  of  aerodynamic  Influence  coefficients  permits  the 
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estimation  of  all  the  stability  derivatives.  Including  static  aeroelastlc 

behavior,  that  are  necessary  to  solve  for  the  trim  condition.  For  the  j 

symmetrical  case,  the  angle  of  attack  and  a  symmetric  trim  mode  permit 

calculation  of  the  longitudinal  stability  derivatives  and,  with  the  velocity, 

ambient  density,  and  load  factor  specified,  the  trim  angle  of  attack  and 

horizontal  stabilizer  (or  elevator)  deflection  are  calculated.  The  load 

distribution  on  the  vehicle  then  follows.  For  the  case  of  a  steady  turn, 

additional  lateral -directional  derivatives  are  required.  These  require 

the  aileron,  rudder  (or  its  equivalent),  and  yaw  rate  modes.  With  the 

i 

l 

addition  of  these  lateral-directional  derivatives  to  the  longitudinal  set, 
the  angle  of  attack,  stabilizer,  aileron  and  rudder  deflections  are  de¬ 
termined  for  a  given  speed,  ambient  density,  load  factor  and  rate  of  cl  in*. 

The  load  distribution  on  the  aircraft  again  follows. 

The  blast  Induced  gust  time  history  is  obtained  for  specified  range 
and  yield  in  the  same  fashion  as  in  the  VIBRA-4  program. 

The  Fourier  transform  of  the  blast  time  history  is  found  numerically 
from  the  Integral  definition  of  the  transform  and  the  assumed  or  specified 
profile  of  the  transverse  blast  velocity.  The  numerical  evaluation  of  the 
Integral  Is  straightforward. 

The  product  of  the  Fourier  transform  of  the  blast  profile  and  the 
frequency  response  of  the  vehicle  modes  to  the  unit  harmonic  travelling 
gust  Is  the  Fourier  transform  of  the  transient  response  to  the  specified 
blast  wave.  The  Inverse  transform  of  the  product  Is  therefore  the  desired 
time  history  of  the  modal  responses.  The  numerical  evaluation  of  the 
Inverse  transform  Integral  Is  also  straightforward  and  Is  similar  to  that 
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for  the  forward  transform  integral.  Superposition  of  the  loads  due  to  maneu-  j  \ 

vering  give  the  final  integrated  load  time  histories  and  stress  time 
histories. 

Numerical  checks  on  the  time  histories  of  the  integrated  loads  are 
made  against  input  allowable  loads  or  stresses  and  an  estimate  is  made  of 
the  allowable  peak  gust  velocity  and  overpressure.  Iteration  then  is 
carried  out  (if  desired)  to  establish  the  critical  range  for  specified 
orientations  and  flight  conditions. 

Two  coordinate  systems  are  used  in  VI BRA-6.  They  are  the  Earth  Fixed 
Axis  System  (EFAS)  and  the  Aircraft  Axis  System  (AAS),  The  AAS  is  used 
for  definition  of  the  aircraft  geometry  and  determination  of  its  modes, 
mass,  and  aerodynamics.  Figure  3  illustrates  the  AAS  system.  The  EFAS  is 
used  to  position  the  AAS  in  space  and  establish  the  burst  to  aircraft 
relations  as  a  function  of  time.  Figure  4  shows  the  EFAS  system.  The 
aircraft  is  initially  positioned  at  the  input  altitude  at  time  zero  (shock 
intercept)  and  out  the  YEpAS  axis  a  distance  Ry  where  Ry  is  the  turn 
radius  (if  the  aircraft  is  In  a  turn).  Time  at  time  of  burst  is  negative. 

Unlike  VIBRA-4,  the  EFAS  is  fixed  at  sea  level  and  the  burst  is  located 
in  space  from  the  orientation  direction  cosines  in  the  AAS  system  and  the 
initial  slant  range  at  time  of  intercept.  The  assumption  has  been  made 
that  at  distances  from  the  burst,  corresponding  to  damage  thresholds  for 
large  subsonic  aircraft,  the  shock  front  appears  planar  to  the  aircraft 
and  the  shock  front  is  moving  at  sonic  velocity. 
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SECTION  III 


MODULE  DESCRIPTIONS 


This  section  describes  the  operation  of  the  various  modules  presently 
in  VIBRA-6.  Contained  within  each  module  description  are  the  equations  for 
the  solutions.  Nomenclature  is  defined  as  used. 

In  some  instances  matrix  operations  are  shown  in  the  figurative  sense 
whereas  the  actual  coding  for  the  most  part  carries  out  such  operations  in 
a  more  efficient  fashion. 
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1.  CONTROL  MODULE 


This  module  serves  to  control  program  data  flow  and  operation.  All 
data  necessary  to  the  modular  execution  and  dimensioning  of  the  modules  are 
read  into  this  module.  Data  to  any  modulehave  three  sources  at  a  given  point 
of  execution:  the  input  data  file  (card  or  disk),  core  stored  data,  or 
generated  data  files  (tape  or  disk).  The  modularization  of  the  program 
allows  the  generation  of  a  data  base  of  files  which  may  be  reused  for  a 
specific  aircraft  configuration  thus  avoiding  the  necessity  of  recalculation 
of  the  majority  of  basic  data  necessary  for  final  solution. 

The  input  data  file  is  divided  into  two  classes:  fixed  data  and  run 
data.  Fixed  data  are  data  which  describe  the  basic  configuration  under  analysis 
and  are  expected  to  vary  the  least,  while  run  data  varies  from  case  to  case. 

The  collection  of  these  data  is  called  the  Fixed  Data  Deck  and  the  Run  Data 
Deck.  While  the  Fixed  Data  Deck  must  be  available  to  the  program  during 
execution,  only  the  data  needed  for  execution  during  any  single  phase  is 
used.  Header  cards  are  used  within  the  Fixed  Data  Deck  to  define  particular 
groupings  of  data. 

The  following  input  data  file  data  comprise  the  Fixed  Data  Deck: 

•  Data  necessary  for  execution  of  the  aerodynamic  module. 

•  Data  necessary  for  calculation  of  the  generalized  mass,  stiffness 
and  damping  in  the  inertial  module. 

•  Data  necessary  for  calculating  unit  loads  matrices  in  the  unit 
loads  module. 

The  following  input  data  file  data  comprise  the  Run  Data  Deck: 
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•  Flight  condition  data- 

•  Data  necessary  to  define  the  active  control  system. 

•  Frequency  data  for  the  frequency  response  module. 

•  Data  for  defining  the  maneuver  condition. 

•  Data  for  definition  of  the  blast  condition. 

Generated  data  consist  of  stored  arrays  which  result  from  execution  of 
four  of  the  modules.  If  these  data  are  available  to  the  program  and  so  noted 
by  the  user,  the  program  draws  upon  the  files  as  necessary.  Core  data  consist 
of  data  which  are  minimal  in  cost  to  generate  and  required  in  numerous  parts 
of  the  overall  program  and  thus  are  not  saved  as  generated  file  data. 

Figure  5  details  the  major  routines  called  by  the  control  module, 
each  of  which  is  the  controlling  routine  of  a  module.  Tables  1  and  2 
illustrate  the  data  requirements  for  each  module  and  the  data  source.  A 
more  complete  description  of  the  data  input,  module  by  module,  is  contained 
in  the  module  descriptions  following.  In  these  descriptions,  input  data  file 
is  referred  to  as  'card'  input  and  generated  data  file  is  referred  to  as 
'tape'  data. 

The  analysis  of  an  aircraft  can  be  carried  out  with  a  single  pass 
through  the  program.  The  suggested  procedure  for  operation  is  to  develop 
aerodynamics  for  a  few  Mach  numbers,  unit  loads  once,  and  cycle  through  the 
active  system,  frequency  response,  and  unit  gust  load  modules  for  as  many 
airspeed/altitude  conditions  as  required.  Then,  cycle  through  the  trim, 
blast  and  time  history  modules  only  for  maneuver  and  blast  conditions  as 
required. 
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Dimensions  core  fixed  data  and  calls 
the  appropriate  modules 

Calculates  the  aerodynamics 
(outputs  aero  file) 


Calculates  the  generalized  mass, 
and  stiffness 


Calculates  integrated  loads  due  to 
unit  modal  amplitudes  and  unit  gust 
forces  (outputs  unit  load  file) 


Calculates  transfer  function  polynomials 
and  sets  up  data  for  active  control 
analyses 


Calculates  generalized  frequency 
response  due  to  unit  gusts 
(outputs  frequency  response  file) 


Calculates  integrated  loads  due  to 
unit  gusts 

(outputs  unit  gust  loads  file) 


Load  and  response  routine 


CTRIM 
TRIM  MODULE 


GUSTDR_ 

BLAST  AND 
TIME 
MODULES 


\  Calculates  steady  loads 
'for  trimmed  flight 


Calculates  blast  time  history 
and  integrated  loads  time 
(histories 


Figure  5.  Control  Module  Routines 
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Calculates  rigid  loads 


Merges  aerodynamic  files 


Figure  5  (cont'd).  Control  Module  Routines 
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TABLE  1 


MODULE  EXECUTION  REQUIREMENTS 
FOR  NEW  DATA  INPUT 


MODULE  R  =  RERUN  REQUIRED 
AFFECTED  FDD  =  FIXED  DATA  DECK  DATA  REQUIRED 
RDD  =  RUN  DATA  DECK  DATA  REQUIRED 


MODULE 

EXECUTED 

FOR  NEW  DATA 

AERO  ;  INERTIAL 

- - 

ACTIVE 

CONTROL 

UNIT 

LOAD 

FREQUENCY 

RESPONSE 

UNIT 

GUST 

LOAD 

BLAST 

RESPONSE 

AERODYNAMIC 

FDD 

R 

R 

R 

R 

R 

R 

INERTIAL 

FDD 

R 

R 

R 

R 

R 

ACTIVE 

CONTROL 

RDD 

R 

R 

R 

UNIT 

LOAD 

FDD 

R 

R 

FREQUENCY 

RESPONSE 

RDD 

R 

R 

UNIT  GUST 

LOAD 

RDD 

R 

BLAST 

RESPONSE 

_ 

RDD 

Note:  Inertial  module  is  always  exercised  with  input  of  run  data. 
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TABLE  2 

BASIC  DATA  REQUIREMENT  AND  TRANSFER 


fed—  mnr 

INPUT  RUU  =  RUN  W 

Fmm - 

TA  DECK 

OUTPUT 

MODULE 

CARD  DATA 

CORE  DATA 

GENERATED 
DATA  FILE 

CORE  DATA 

GENERATED 
DATA  FILE 

AERODYNAMIC 

(FDD) 

AERC  GEOM. 
AERO  MODES 
ORIENTATIONS 
REDUCED  FREQ 
MACH  NO. 

AERO  GEOM. 
D<t>c  D*fl 

F  ’F  H 

V  A,s 

HS  ,Sp£hA 

INERTIAL 

(FDD)- 
MASS  GEOM. 
MASS 

FREQUENCIES 
INERTIA  MODES 
MODE  DEF. 

m 

K 

g 

M<j>, 

MODE  DEF. 

UNIT  LOAD 

(FDD)'~ 

BEAM  GEOM. 
INT'GD  LOAD 
GEOM. 

LOAD  ALLOW. 
STRESS  DEF. 
MASS  SPARSING 
MASS  ORIENT. 
BOX  SPARSING 
BODY  SPARSING 
ENG.  GEOM. 

MASS  GEOM. 

*1 

M*j 

AERO  GEOM. 

w$ » 

PIQ 

PAQS 

PAQA 

THRGNF 

THRL0D 

PINT 

STRESS 

ACTIVE 

SYSTEM 

f  (RDD) 

ACS  KINEMATIC 

BLOCK  DATA 

;  ♦! 

hh 

JM 

FREQUENCY 

RESPONSE 

(RDD) 

FREQUENCIES 

ALTITUDE 

AIRSPEED 

m 

K 

g 

V  ta 

ACS 

KINEMATICS 
MODE  DEF. 

Zk’  ^A 

Fs*  fa 

SPLHS ,SPLHA 

q 

UNIT  GUST 
LOADS 

i _ 

Vta 

ACS 

KINEMATICS 

q 

PIQ 

PAQS 

PAQA 

Fs*  fa 

PINT 

STRESS 

p  p 

V  KA 

STRESS 
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TABLE  2  (cont) 


BASIC  DATA  REQUIREMENT  AND  TRANSFER 


MODULE 

CARD  DATA 

CORE  DATA 

GENERATED 
DATA  FILE 

CORE  DATA 

TRIM 

(RDD) 

MANEUVER  DAT/ 

THRUST 

MODE  DEF. 

m 

a 

'7j  X> 

'S’ 

PIQ 

PAQS 

PAQA 

THRGNF 

THRL0D 

PSTRIM 

PATRIM 

BLAST  AND 
TIME 

RESPONSE 

- TO - 

ORIENTATIONS 
VIELD  &  RANGE 
MAX  TIME 
ALLOWABLES 

PSTRIM 

PATRIM 

PS*  PA 

STRESS 

P 

a 

2.  AERODYNAMIC  MODULE 

The  aerodynamic  module  SDLM  has  been  modified  from  the  computer  program 
of  Reference  2.  The  modifications  to  add  a  travelling  gust  field,  improve 
the  aerodynamic  influence  coefficients,  and  change  the  load  output,  are  dis¬ 
cussed  in  detail  in  Volume  II  of  this  report.  The  basic  data  requirements 
and  output  are  summarized  here  for  convenience.  This  module  serves  to  pro¬ 
vide  the  aerodynamic  forces  due  to  motion  and  the  travelling  gust  and  the 
generalized  aerodynamic  equations  for  solution  of  the  blast  response  problem 

The  geometry  required  is  that  necessary  to  describe  the  lifting  surfaces, 
the  slender  bodies,  and  the  interference  surfaces  between  the  lifting  surfaces 
and  bodies.  The  lifting  surfaces  are  idealized  as  plane  panels  with  no  thick¬ 
ness,  camber,  or  twist,  but  with  dihedral.  Each  lifting  surface  is  subdivided 
into  smaller  trapezoidal  lifting  elements  (boxes)  arranged  in  strips  parallel 
to  the  freestream  such  that  box  boundaries  lie  on  surface  edges,  fold  lines, 
and  control  surface  edqes.  For  coplanar  or  near-coplanar  surfaces,  e.  g., 
a  wing  and  tail,  the  boxes  must  be  aligned  in  the  streamwise  direction.  For 
non-coplanar  surfaces  alignment  is  not  required  if  the  perpendicular  separation 
is  more  than  one  strip  width.  For  intersecting  surfaces  the  box  edges  must 
also  be  located  at  intersections,  as  in  the  case  of  a  wing  pylon. 

The  idealization  of  fuselages,  nacelles,  or  external  stores  as  slender 
bodies  does  not  require  geometric  similarity  to  the  actual  body,  but  it  is 
recommended  that  the  best  elliptical  representation  of  the  body  cross-sectional 
area  be  used.  The  displacements  of  the  body  during  motion  are  accounted  for 

in  the  upwash  and  sidewash  boundary  conditions.  All  bodies,  including  jet 
engines,  must  be  idealized  as  having  pointed  noses;  however,  a  body  need  not 
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be  closed  at  its  downstream  end  and  a  suitable  base  area  may  be  selected  to 
approximate  flow  separation  effects.  In  addition  to  the  idealization  of  the 
body  cross  sections,  the  body  length  is  divided  into  segments.  Each  segment  is 
described  by  its  length  and  width  at  each  end.  Smaller  lengths  are  chosen 
where  the  body  cross  section  is  changing  rapidly  and  in  regions  of  maximum 
interference  with  lifting  surfaces  such  as  occur  near  the  wing  root  or  at 
a  nacelle-pylon  intersection. 

The  interference  surface  is  a  cylindrical  tube  with  an  elliptical  cross 
section.  All  attached  lifting  surfaces  are  connected  to  this  tube.  Within 
this  tube  is  placed  a  system  of  images  of  the  external  lifting  elements,  and 
the  image  system  approximately  negates  the  flow  field  at  the  body  surface 
induced  by  lifting  surfaces.  The  interference  tube  is  segmented  into  elements 
by  means  of  angular  divisions  around  its  circumference  and  streamwise 
divisions  along  its  length.  Since  the  interference  elements  increase  the 
computational  cost,  their  number  should  be  limited,  but  they  need  only  be  used 
near  lifting  surface/body  intersections  and  only  upstream  and  downstream  with¬ 
in  a  chord-length  of  the  lifting  surface.  In  situations  in  which  high  accuracy 
is  not  required,  e.g.,  for  small  external  stores,  the  interference  elements 
might  be  eliminated  altogether. 

In  addition  to  the  geometry  of  the  idealized  aerodynamic  configuration, 
the  aerodynamic  load  calculation  requires  the  normalwash  distribution  at  the 
aerodynamic  elements.  The  aerodynamic  control  points  are  located  differently 
on  the  lifting  surfaces  and  on  the  slender  bodies.  On  the  lifting  surfaces, 
the  normalwash  control  point  is  chosen  at  the  three-quarter  chord  po^nt  of 
the  centerline  of  each  box  and  the  lift  acts  at  the  one-quarter  chord  point. 

On  the  slender  body  elements,  both  the  normalwash  and  lift  (or  side  force) 
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are  determined  at  the  midpoint  of  the  segment  length.  The  deflections  and 
slopes  at  these  points  are  determined  by  a  surface  spline  interpolation  among 
a  set  of  deflections  of  aerodynamic  modes.  The  aerodymamic  modes  differ  from 
the  structural  modes  to  the  extent  that  they  cover  the  lifting  surfaces  and 
slender  bodies  more  completely,  e.g.,  including  the  leading  and  trailing  edges, 
whereas  the  structural  modes  give  the  deflections  only  at  the  mass  points.  The 
increased  area  covered  by  the  aerodynamic  modes  improves  the  accuracy  of  inter¬ 
polation  for  the  normalwashes  and  deflections  and,  hence,  the  accuracy  of  the 
generalized  aerodynamic  forces.  The  aerodynamic  modes  correspond  to  the 
vehicle  rigid  body  and  elastic  modes,  the  control  surface  and  trim  modes, 
and  jig  modes.  The  jig  modes  are  necessary  to  establish  the  basic  pitching 
moment  and  basic  lift  distributions  over  the  vehicle,  and  they  describe  the 
positions  from  the  reference  aerodynamic  planes  of  the  actual  chord  lines  of 
surfaces,  (camber  and  twist),  and  local  body  pitch  and  surface  pitch. 

The  aerodynamic  influence  coefficients  (AIC's)  of  Reference  2  are  useful 
either  in  a  direct  influence  coefficient  solution,  which  avoids  modal  convergence 
problems,  or  in  the  modal  solution  which  consider  many  varying  modes,  e.g.,  modes 
that  change  because  of  design  changes  in  stiffness  or  weight  distribution. 

Since  the  AIC's  depend  on  the  planform  (and  Mach  number  and  reduced  frequency) 
and  not  on  the  modes,  the  generalized  aerodynamic  forces  are  easily  obtained 
from  the  AIC's, 

[y-]  =  oa]t  [aic]  ua] 

where  [<f>A]  is  the  matrix  of  modal  deflections  at  the  aerodynamic  nodal  points. 

The  additional  calculations  required  to  obtain  the  AIC's  are  not  justified 
for  the  VIBRA-6  analysis  because  only  a  few  sets  of  modes  are  considered  in  a 
response  analysis.  The  modifications  to  Reference  2  described  in  Volume  II 
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of  this  report  are  summarized  below  along  with  a  few  basic  relationships 
necessary  for  use  of  the  finite  element  aerodynamics. 

Let  [SPL ]  be  the  spline  interpolation  matrix  that  yields  the  complex 
normalwash  (w)  from  the  deflections  {h^l  at  the  aerodynamic  nodal  points 
(w)  =  [SPL]  {hft} 

The  nodal  point  deflections  (h^}  may  be  regarded  as  either  translations 
or  rotations.  The  modal  solution  expresses  the  deflections  as  a  series 
of  vibration  modes  with  generalized  coordinates  (ql: 

{hA]  =  lV  fq’ 

and  the  modal  approximation  to  the  normalwash  is 
{wq}  =  [SPL]  [(f>A]  {q} 

The  Doublet-Lattice  Method  initially  relates  the  downwash  to  the  local 
pressure  coefficients  by  the  relationship  (Ref.  5  and  6)  through: 

!V  *  |:w  !4V 

where  [Ddlm]  is  a  large  matrix.  The  downwashes,  however,  are  given  in 
terms  of  a  reduced  set  of  coordinates  q  as  shown  above.  The  forces  on 
the  boxes  and  slender  bodies  are  given  in  terms  of  the  pressure  coef¬ 
ficients  by: 

-  ^PTP3 

where  [FDTp]  is  an  integration  matrix  relating  local  force  to  local 
pressure.  These  forces  then  are  obtained  by  the  general  relationship: 

^B'  =  ^DTP3  3°DLM3  fSpl.]  [,^]  {q} 

5.  Alba  no,  E.,  Rodden,  W.P.,  A  Doublet  Lattice.  Method  fan  Calculating 

U*t  Vis  tn.ibuti.oni  on  Oscillating  Sun.  faces  in  Subsonic  Flows,  AIAA  Journal 
Vol.  7,  No.  2,  February  1969. 

6.  Kalman,  T.P.,  Rodden,  W.P.,  Giesing,  J.P.,  Aenodijnamic  Influence  Coeff¬ 
icients  by  the  doublet  Lattice  Method  fan  In-ten/fening  Nonplanan  Lifting 
Sun  faces  Oscillating  in  a  Subsonic  Flow,  IRAD  Final  Report,  DAC-67977. 
November  1969. 


where  the  inverse  is  not  actually  carried  out,  but  the  forces  due  to  unit 
q  obtained  by  solving  for  the  number  of  right-hand  sides  corresponding  to 
the  number  of  modes  required.  The  generalized  forces  are  given  in  terms 
of  the  above  forces  by  the  relationship: 

=  [$ait  [SPLH]  fFB) 

where  [SPLH]  is  an  interpolation  matrix  relating  local  force  to  aerody¬ 
namic  nodal  point  forces.  The  stored  generalized  aerodynamic  forces  due 
to  motion  then  are  given  symbolically  by 

lZ)\  =  (*a]T  [SPLH]  [FDTp]  [D^f1  [SPL]  [4>a] 
and  the  stored  local  aerodynamic  forces  due  to  motion  (necessary  for  in¬ 
tegrated  loads)  by: 

[M>]  =  [fdtp1  tw"1  [SPL]  [*A] 

The  gust  loads  are  calculated  for  each  orientation  of  the  aircraft 
and  the  blast.  The  unit  normalwash  wg  induced  by  the  blast  at  a  point 
on  the  aircraft  is  given  (in  Vol.II)  by 

wg  (x,  y,  z)  =  GQh  exp  (  -  i  2  R/c) 

where 

Co  =  the  amplitude  of  the  harmonic  blast  material  velocity 

->  .  -> 

o  =  n  vg 

n  =  the  unit  normal  vector  to  the  lifting  surface  or  slender  body 
vg  =  i  cosa  +  j  cosB  +  k  cosy 
a,3,y,  =  the  blast  orientation  angles 

kr  =  0JC/2U 

ai  =  circular  frequency 
c  =  reference  chord 
U  =  aircraft  velocity 
9.  =  x  cosa  +  y  cosR  +  z  cosy 
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R  =  U/(V  +  U  cosa) 

V  =  velocity  of  blast  pressure  wave 

y 

In  general,  an  arbitrary  burst  location  results  in  an  asymmetric  loading 
of  the  flight  vehicle.  Since  any  asymmetric  condition  can  be  found  as  a 
linear  combination  of  a  symmetric  and  an  antisymmetric  condition  (assuming 
linear  systems),  both  symmetric  and  antisymmetric  aerodynamic  gust  analyses 
must  be  performed.  If  the  gust  normalwash  on  the  right  side  is  wg(x>  y»  z»  y) 
and  on  the  left  side  is  wg(x,  -y,  z,  -7),  then  the  symmetrical  normalwash  is 

wgs  =  (1/2)  [w  (x,  y,  z,  y)  +  wg(x,  -y,  z,  -y)] 

and  the  anti  symmetrical  normalwash  is 

w_,  =  (1/2)  [w  (x,  y,  z,  y)  -  w  (x,  -y,  Z,  -y)] 

9  y 

The  symmetric  and  antisymmetric  normalwashes  can  be  used  as  outlined 
above  to  obtain  sets  of  symmetrical  and  anti  symmetrical  local  forces  {Fg}  and 
generalized  forces  {r9g}. 

{Fg}  =  ^FDTP-*  ^DLM^  {V 

{^g>  =  C^a3T  CSPLH3  {Fgl 

The  aerodynamic  data  required  for  the  blast  response  analysis  consist 
therefore  of  generalized  aerodynamic  forces  for  symmetric  and  antisymmetric 
modes  of  vibration,  and  for  symmetric  and  antisymmetric  harmonic  gust  fields 
having  a  number  of  orientations.  The  orientations  depend  on  the  burst  locations. 
A  standard  set  of  thirteen  burst  locations  is  shown  in  Figure  6.  Other 
orientations  may  be  chosen  but  must  be  specified  separately  by  the  user. 

The  generalized  forces  are  required  at  a  sufficient  number  of  reduced  frequen¬ 
cies  that  interpolation  will  lead  to  accurate  values  at  intermediate  frequencies. 
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flgure  6.  Blast  Orientations  for  Aero  Module 


The  spacing  of  the  reduced  frequencies  necessary  for  an  accurate  interpolation 
depends  on  the  blast  orientation;  an  estimate  is  given  as: 

£kr  =  irp/ u i  '  P  =  (1  +  PL  cosoL/PLi/c 

where  l  is  the  slant  distance  from  a  point  on  the  aircraft  to  the  blast  plane, 
as  the  blast  plane  passes  through  the  aircraft  origin.  If  the  maximum 
distance  l  is  used  the  estimate  for  Akr  will  be  conservative  for  all  points 
on  the  aircraft. 


The  maximum  reduced  frequency,  k  ,  that  is  necessary  to  consider  is 


max 


estimated  as: 


'max 

where  l2  is  approximately  2.5.  In  many  cases  the  k  value  is  above  the 

rmax 

practical  limit  for  the  Doublet  Lattice  Method  (which  is  2  or  3).  For  such 

cases  piston  theory  is  used  to  fill  the  gap  between  the  DLM  krmax  ar|d  the  one 

required  by  the  above  formula.  The  application  of  piston  theory  is  discussed 

in  Vol .  II,  Section  1 1 -4 ,  and  is  based  on  the  following  relation: 

■Xp  =  (4/M,)  W/V 

where  W/V  is  the  dimensionless  downwash. 
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3.  INERTIAL  MODULE 


The  inertial  module  calculates  the  generalized  mass  and  stiffness  from 
the  specified  inertial  mode  shapes  and  point  mass  data.  Input  data  which  are 
used  for  core  allocation  consist  of  the  number  of  masses  for  the  half  aircraft 
analysis  and  the  number  of  symmetric  and  antisymmetric  modes  used  (which  must 
agree  with  the  aerodynamic  module  input.)  Array  data  input  consists  of: 

•  The  point  mass  data  for  the  right  side  of  the  aircraft  and  the  point 
coordinates  in  the  AAS  system.  Though  the  analyses  are  carried  out 
for  only  half  an  aircraft,  full  centerline  mass  point  data  are  loaded. 

*  Inertial  mode  shapes  and  associated  frequencies  defining  the  motion 
in  three  orthogonal  axes  at  each  point  (identified  as  PHIX,  PHIY, 
and  PHIZ.) 

*  An  array  called  AM0DN0  which  defines  the  type  of  modes  by  location 
in  the  modal  array.  The  input  requirement  is  set  that  all  symmetric 
modes  precede  all  antisymmetric  modes  and  the  order  of  inertial  mode 
input  within  these  general  groupings  must  agree  with  the  aerodynamic 
mode  ordering  sequence.  The  rigid  body  modes  must  have  zero  frequency 
input  and  input  modes  defining  any  jiq  modes  must  be  present,  but 

may  be  of  zero  deflection. 

•  The  structural  damping  of  each  mode,  which  should  be  zero  for  all  but 
the  elastic  modes. 

The  inertial  mode  shapes  (as  distinguished  from  the  aerodynamic  mode 
shapes)  define  the  mass  point  motion  in  each  mode.  The  motion  of  each  point 
must  consist  of  three  orthogonal  modal  deflections  but  they  need  not  necessarily 
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be  deflections  parallel  to  the  AAS  system.  (Input  data  In  the  unit  loads 
module  are  provided  to  establish  the  alignment  of  these  deflections  in  terms 
of  AAS  system  motion.)  The  modes  need  not  be  aircraft  free-free  modes,  as 
the  analysis  procedure  will,  in  effect,  orthogonalize  the  modes  in  the  event 
that  constrained  modes  are  used.  A  distinction  is  made,  however,  between 
modes  used  to  describe  the  symmetric  and  antisymmetric  motion  of  the  aircraft. 

There  are  two  formats  for  the  mass  data.  The  first  simply  consists  of 
the  concentrated  masses  at  each  mass  point.  The  second  consists  of  the 
sectional  mass  data  for  each  section  when  the  modes  are  specified  by  the 
user  in  terms  of  sectional  translations  and  slopes  (called  the  h  and  a 
input).  The  first  format  is  a  subset  of  the  second  although  the  two  have 
been  programmed  separately.  The  sectional  mass  matrix  Is  shown  in  Figure  7 
for  the  maximum  number  of  eight  degrees  of  freedom.  Figure  8  shows  the 
local  coordinate  systems  of  the  mass  data  and  modeshapes  for  typical 
components.  The  mass  matrix  itself  is  not  input  to  the  program.  The 
mass  properties  are  input  and  the  program  forms  the  mass  matrix  intern¬ 
ally  from  the  mass  properties. 

The  generalized  mass  matrix  is  given  by 

Cm]  «  [n]T  [M]  [n] 

where  Ui]  is  the  array  of  all  inertial  shapes.  The  off-diagonal  terms  of 
the  generalized  mass  matrix  coupling  symmetric  and  antisymmetric  modes  are 
zeroed  since  the  assumption  of  a  vertical  plane  of  symmetry  in  mass,  stiff¬ 
ness  and  geometry  is  made,  and  thus  this  coupling  is  identically  zero  for 
free-free  modes  but  not  arithmetically  for  the  one  half  aircraft  analyses. 

The  generalized  stiffness  is  given  by 
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local  coordinate  system 
right-hand  side  of  aircraft 


Figure  7.  Sectional  Mass  Matrix 
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Figure  8.  Component  Local  Coordinate  Systems 


ki  *  r-11  wi2J 

The  generalized  mass  and  stiffness  matrices  are  saved  in  core  along  with  the 
input  structural  damping  and  mode  definition  array  for  subsequent  use.  The 
array  EMPHI  which  gives  the  Inertial  forces  at  the  mass  points  for  unit 

generalized  response  is  calculated  by 
[EMPHI]  *  [M]  [♦!] 

for  use  in  calculation  of  the  Integrated  loads  due  to  mass  point  motion. 
Figure  9  shows  the  principal  routines  of  this  module. 


Controls  module  operation 

Reads  input  data  for  module 

Calculates  generalized  mass,  stiffness 
and  zeroes  appropriate  arrays 


Figure  9.  Inertial  Module  Routines 
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4.  UNIT  LOADS  MODULE 


The  unit  loads  module  serves  to  generate  the  matrices  required  to  obtain 
integrated  loads  at  all  required  locations  on  the  vehicle.  Its  primary 
functions  are  to: 

*  Find  the  spatial  location  and  orientation  of  each  of  the  local  beams, 

the  network  of  which  defines  the  load  integrations  for  integrated 
loads  and  the  spatial  orientations  of  the  integrated  loads. 

*  Calculate  the  integrated  loads  due  to  unit  inertial  loads  at  the  mass 

points  and  due  to  unit  aerodynamic  loads  (motion  dependent  and  qust) 
at  the  aerodynamic  box  load  points  and  slender  body  load  points. 

‘  Calculate  the  integrated  loads  due  to  unit  generalized  responses  in 

the  modes  for  loads  from  mass  point  motion  and  aerodynamic  box  and 

slender  body  motion  and  thrust. 

Figure  10  details  the  major  routines  of  this  module  and  a  brief  description 
of  their  functions.  Input  data  required  for  this  module  consist  of: 

’  An  aero  file. 

'  The  number  of  local  beams  describing  the  possible  load  paths  for 
internal  loads,  the  AAS  coordinates  of  the  beam  end  points  and 
their  component  definition  numbers. 

'  The  number  of  integrated  loads  desired,  the  local  beam  number  each 
load  is  associated  with,  a  load  and  component  code  for  type  and 
orientation,  the  AAS  coordinate  of  each  integrated  load  and  the 
maximum  allowable  positive  and  negative  load  for  each. 
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LOAD 


Dimensions  arrays  and  controls  module 
operation 


INPUTL 


Reads  in  required  data  for  unit  load 
matrices  calculation 


Defines  load  paths  and  spatial  orienta¬ 
tions  of  integrated  loads 


f  GE0MMB 

V  MMBMLT  , 

Z'"'" 

1  THRUST  ! 


Forms  integrated  loads  due  to  unit 
inertia  forces 


Forms  integrated  loads  due  to  unit  thrust 
and  generalized  and  symmetric  forces 
due  to  unit  thrust 


Checks  for  centerline  integrated  loads 
and  forms  symmetric  and  antisymmetric 
components 

Forms  integrated  inertial  loads  due  to 
unit  modal  amplitudes 


/^~GE0MAB 
V  ABAMLT  J 


Forms  integrated  loads  due  to  unit 
aerodynamic  panel  loads 


GE0MSB 


SBAMLT  J 


AR0L0D 

AR0PHA 


Forms  integrated  loads  due  to  unit 
slender  body  aerodynamic  loads 


Forms  integrated  aerodynamic  loads 
due  to  unit  modal  amplitudes 


Figure  10.  Unit  Loads  Module  Routines 
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•  the  number  of  stresses  required  (if  any)  and  the  matrix  relation  between 


stresses  and  integrated  loads. 

•  the  number  of  mass  point  groups,  the  first  and  last  mass  point  number 

and  local  beam  number  the  loads  from  each  qroup  of  mass  points  first 

enter  and  direction  cosines  (TLAMV)  relating  inertia  motion  to  AAS  motion. 

•  the  number  of  aerodynamic  box  groups  and  the  first  and  last  box  numbers 

and  local  beam  number  the  loads  from  each  group  of  aerodynamic  boxes 

first  enter. 

•  the  number  of  slender  body  element  aerodynamic  groups  and  the  first 
and  last  body  element  and  local  beam  number  the  loads  from  each  aroup 
of  body  elements  first  enter. 

•  number  of  engines  and  engine  thrust  location  definition. 

Routine  LOAD  dimensions  and  controls  the  operation  of  this  module. 

INPUTL  reads  the  required  external  input  data  which  consist  of  beam  network 
geometry,  integrated  load  definition  data,  stress  calculation  matrix  (if 
required),  mass  group  definitions,  aero  box  group  definitions,  and  aero  body 
group  definitions.  Group  definitions  of  mass  and  aero  points  are  used  to  form 
sparse  matrices  and  to  define  which  beam  each  group  is  associated  with. 

Routine  BEAM  defines  the  load  paths  by  checking  connections  of  the  beam 
network  and  calculates  the  orientation  in  space  of  the  beams. 

A  local  beam  coordinate  system  is  used  to  define  the  orientation  of  the 
integrated  loads.  The  location  in  space  of  any  beam  is  determined  from  the 
end  points  of  the  beam  defined  as  the  inner  and  outer  ends.  Load  integration 
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always  passes  from  the  outer  ends  to  the  inner  ends  of  the  local  beams.  The 
connection  of  these  beam  ends  then  defines  the  load  path  for  determination 
of  the  integrated  load  which  is  the  summation  of  all  equilibrium  loads  'out¬ 
board'  of  the  integrated  load  station. 

The  spatial  location  and  orientation  of  a  beam  segment  is  determined  as 
follows: 

•  Let  the  beam  segment  be  initially  oriented  parallel  to  the  AAS  system, 
with  the  beam  Y  axis  along  the  AAS  Y  axis.  The  final  orientation  is 
determined  by  rotating  through  dihedral  (positive  about  the  X  axis) 
then  sweep  (negative  about  the  dihedrally  rotated  Z  axis).  The  beam 
segment  is  now  properly  rotated  and  its  end  points  are  now  as  given  by 
its  endpoint  coordinates  in  the  AAS  system. 

•  The  direction  cosine  and  transfer  matrix  relating  the  two  coordinate 

systems  (XR  the  reference  or  AAS  system  and  X£  the  beam  local  system) 

%  % 
are  given  by 


where  XJR  are  the  coordinates  of  the  inner  end  of  the  beam  segment 
given  in  the  AAS  system  and  x$,  us,  Xp,  are  the  cosines  and  sines 
of  the  sweep  and  dihedral  angles  respectively. 

The  AAS  coordinates  of  the  outer  end  of  any  beam  are  given  by 
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which  is  required  in  the  load  calculations  to  insure  that  any  equilibrium 
load  which  is  inboard  of  an  integrated  load  station  is  not  included  in  the 
integration.  The  above  rotation  arrays  are  generated  in  BEAM. 

The  integrated  loads  are  defined  by  specifying  which  local  beam  each 
is  associated  with,  a  load  code  which  identifies  it  as  one  of  the  six  possible 
shears  and  moments,  the  aircraft  component  it  is  associated  with,  its  spatial 
location  in  the  AAS  system  and  the  maximum  positive  and  negative  values.  The 
location  in  space  of  the  integrated  load  need  not  be  on  a  beam,  but  its 
orientation  will  be  vectorially  parallel  to  the  associated  local  beam  system 
and  load  transfer  will  be  made  in  a  plane  normal  to  the  local  Y  axis  and 
containing  the  load  point. 

The  integrated  loads  due  to  inertial  loads  are  found  by  transferring 
each  inboard  along  the  beam  network.  A  sparsing  procedure  is  employed  in 
input  data  in  that  the  local  beam  a  group  of  mass  point  loads  first  enters 
is  defined  as  well  as  the  mass  number  range  (first  to  last)  of  the  mass  point 
group. 


=  LTLAMYJ 


The  aerodynamic  loads  are  similarly  defined  by  entering  first  to  last 
box  numbers  and  first  to  last  slender  body  element  numbers  and  the  beams  they 
are  first  associated  with.  Box  number  sequencing  and  slender  body  sequencing 
definitions  are  found  in  Volume  II  describing  the  aerodynamic  module  details. 

The  process  of  transferring  a  group  of  equilibrium  loads  to  an  integrated 
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load  station  consists  of:  1)  converting  the  equilibrium  loads  in  the  AAS 
system  to  shears  and  moments  at  the  integrated  load  point  in  the  AAS  system; 
2)  rotating  the  integrated  loads  into  the  local  beam  coordinate  system;  3) 
rotating  the  resulting  integrated  loads  properly  to  reflect  the  orientation 
of  the  equilibrium  loads. 

The  load  transfer  is  given  by 
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system  and 
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where  the  transfer  distances  in  the  AAS  system  from  the  loads  to  integrated 
loads  stations  are 

\XB  =  XRL  -  XRB>  \Yg  *  Yrl  -  Yrb>  AZb  =  ZRL  -  ZRB 

If  sectional  data  is  not  input  to  the  program,  then  MXR  =  0. 

The  integrated  load  rotation  from  AAS  to  local  beam  system  is  found  from 


PXB 
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TLAMM 


I 


Where  now  the  P  and  M  are  shears  and  moments  in  the  local  beam  system 
coordinate  directions  (pounds  and  inch  pounds). 


The  equilibrium  inertial  loads  come  from  rotating  the  inertial  loads  from 
the  orientations  of  the  nodal  deflection  vectors  used  to  the  MS  system. 
Let 


where  and  are  inertia  loads  and  moments  in  whatever  orientation 
the  idealization  has  assumed.  (Note  that  often  TLAMV  will  be  the  transpose 
of  TLAMM  for  corresponding  sets  of  mass  points  and  beam  these  first  enter). 
Then 
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and  TLAMV  is  input  by  the  user  for  each  mass  group,  and  is  seen  to  be  the 
direction  cosine  matrix  giving  the  components  of  the  inertial  loads  in  the 
MS  system. 

The  panel  aerodynamic  loads  are  defined  positive  (for  a  panel  loaded 
in  the  usual  fashion)  up  for  horizontal  surfaces  and  along  the  -YR  axis 
for  vertical  type  surfaces.  For  an  aerodynamic  panel  load  F,  let 


Panel  aerodynamic  loads  are  assumed  to  act  at  the  local  mid-span  point  of 

each  box  and  on  the  local  quarter  chord,  and  y  y  Y  here  are  the  direction 

x’  y’  z 

cosines  of  the  force  in  the  AAS  system. 

Then 


where 


mD 


+  Az 


■  51 


S‘ 


2  2 
Ay  +  AZ 


and  Ay  and  az  are  distances  from  inner  to  outer  ends  of  the  local  aerodynamic 
box  in  the  AAS  system  and  are  obtained  from  the  aerodynamic  module  data. 


Slender  body  aero  forces  are  aligned  in  the  +  ZR  and  +  YR  directions 
only,  and  act  at  the  geometric  center  of  each  body. 

The  spatial  orientation  of  an  integrated  load  is  seen  to  be  fixed  by 
the  rotation  process  the  local  beam  has  undergone  to  arrive  at  its  spatial 
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position  as  defined  previously.  Figure  11  shows  the  definition  of  the 
potential  loads  at  a  load  station,  and  Figure  12  gives  typical  orientations 
for  components  of  the  vehicle.  The  load  station  itself  is  defined  by 
specifying  its  AAS  system  coordinates  and  which  local  beam  it  is  associated 
wi  th. 

Integrated  loads  on  the  centerline  of  the  vehicle  are  checked  to 
establish  their  symmetric  and  antisymmetric  components  and  suitably  modified 
for  the  one-half  aircraft  analysis. 

Routine  GEOMMB  and  MMBMLT  form  the  matrices  PINTH,  PINTL,  PINTS 
relating  the  equilibrium  inertial  loads  to  the  integrated  loads  based  on  the 
preceding  transformation.  If  sectional  data  is  input  to  the  program,  these 
routines  also  form  the  matrices  PINTT,  PINTA,  and  PINTP  for  the  sectional 
slope  degrees  of  freedom.  Routine  MSSPHI  forms  the  matrix  PIQ  relating 
generalized  response  to  integrated  loads. 

I 

hi  ■ 

where 

[EMPHl]  =  [M] 

PIQ  is  saved  for  use  in  the  unit  gust  loads  module. 

Routines  GEOMAB,  ABAMLT,  GEOMSB ,  and  SBAMLT  form  the  matrices  PINTP, 
PINTY,  and  PINTZ  which  relate  unit  aerodynamic  forces  at  all  box  and  slender 
body  load  points  to  the  integrated  loads  using  the  transformations  discussed 
above.  Routines  AROLOD  and  AROPHA  form  the  matrices  PAQS  and  PAQA  which 
give  the  integrated  loads  for  symmetric  and  antisymmetric  motion  dependent 
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NO.  =  LOAD  ORIENTATION  CODE 


Typical  Wing  or  Horizontal  Stabilizer  Loads 
RHS  -  Aircraft  Plan  View 
Component  Code  1 


Component  Code  2 

Figure  12.  Sample  Integrated  Load  Orientations 
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Typical  Vertical  Stabilizer  Loads 
RHS  -  Aircraft  Side  View 
Component  Code  3 


Typical  Wing  Pylon  Loads 
RHS  -  Aircraft  Side  View 
Component  Code  4 


Figure  12  (cont).  Sample  Integrated  Load  Orientations 
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INNER 

5 


OUTER 


Typical  Fuselage  Side  Pylon  Pod  Loads 
RHS  -  Aircraft  Plan  View 
Component  Code  5 


OUTER 


INNER 


Typical  Centerline  Vertical  Stabilizer  Pylon  Loads 
RHS  -  Aircraft  Side  View 
Component  Code  6 


Figure  12  (cont).  Sample  Integrated  Load  Orientations 
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aerodynamic  forces,  respectively,  due  to  unit  generalized  amplitudes. 


PAQS(k) 


PINTP 

PINTY 

PINTZ 

PINTP 

PINTY 

PINTZ 


Where  and  are  respectively  the  motion  dependent  aerodynamic  box  and 
body  forces  due  to  all  symmetric  and  antisymmetric  input  modes.  The  matrices 
PAQS,  PAQA,  PINTP,  PINTY,  and  PINTZ  are  saved  for  use  in  the  unit  gust  loads 


and  trim  modules . 


A  stress  matrix  may  be  defined  by  the  user  which  relates  local  stresses 
to  the  defined  integrated  loads: 

{0}  =  [STRESS]  {P} 

Routine  THRUST  calculates  the  integrated  loads  and  generalized  forces 
due  to  unit  thrust  loads  for  all  engines  specified.  The  components  of  the 
thrust  loads  are  found  from  direction  cosines  yz)  determined  by  the 

line  of  action  of  the  thrust  which  is  found  from  the  two  mass  points  per 
engine  thrust  specified  in  the  input.  If  only  one  mass  station  per  engine 
is  used  with  the  sectional  input  data,  then  a  dummy  station  must  be  specified 
to  define  the  line  of  action  of  the  thrust.  The  resulting  forces  are  taken 
at  the  first  input  thrust  mass:  ^ 
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This  aligns  the  forces  with  the  inertial  forces  at  the  thrust  mass  point. 
The  integrated  loads  due  to  the  ith  unit  thrust  then  are 

F 

F 

F 

The  generalized  forces  in  the  symmetric  modes  are  given  by  premultiplying 
Fxt  Fyt  Fzt  by  the  symmetric  modal  row  corresponding  to  the  thrust  mass: 

PHIX. 

i 

PHI  Y.. 

PHIZi 

'  '  i 

THRL0D  and  THRGNF  are  saved  for  later  use  in  the  trim  module. 


{THRGNF.j 


5.  ACTIVE  CONTROL  MODULE 


The  active  control  module  forms  the  data  and  arrays  necessary  to 
augment  the  equations  of  motion  for  the  feedback  effects.  Input  data  required 
for  this  module  consist  of;  1)  definition  of  the  kinematics  of  the  control 
system  and  2)  'block'  data  defining  the  individual  transfer  functions.  Data 
for  symmetric  and  antisymmetric  control  are  entered  separately.  The  kinematics 
are  described  by  the  definition  (for  each  symmetric  and  antisymmetric  transfer 
function)  of  the  mass  point  number  and  orientation  (x,  y,  z)  of  the  sensed 
motion,  a  scalar  for  scaling  the  transfer  function  and,  the  mode  number  of 
the  input  modal  column  which  provides  the  control  force.  The  ‘block1 
input  data  consist  of  the  definition  of  numerator  and  denominator  polynomials 
which,  when  multiplied  together,  give  the  total  output/input  polynomial  in  s, 
where  s  is  the  LaPlace  operator.  The  maximum  order  of  any  single  block 
polynomial  in  the  entire  group  is  input  as  MX0BLK,  and  the  maximum  number  of 
blocks  in  the  largest  transfer  function  as  MXBLK.  NTFS  and  NTFA  define  the 
total  number  of  symmetric  and  antisymmetric  transfer  functions.  Figure  13 
details  the  major  routines  in  this  module. 

The  incremental  commanded  motion  of  a  set  of  degrees  of  freedom  by  a 
sensor  may  be  related  to  the  sensed  quantity  by  a  transfer  function  T(s). 

For  the  frequency  response  solutions,  s  may  be  replaced  by  i<*\  For  illustrative 
purposes,  let  the  degrees  of  freedom  be  those  deflections  associated  with  the 
rotation  of  a  control  surface. 

=  He^  6 

where  are  the  h^  deflections  due  to  a  unit  rotation  6.  The  commanded 
rotation  by  an  active  systen  then  may  be 


r  ■  iTH< ,  jtiib 
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Dimensions  arrays  and  calls  ACSMOD. 


Controls  routine  operation  and  reads  in 
control  system  definition  data. 


Forms  special  modal  arrays  from  input 
inertial  mode  shapes. 


TFBLIN  ) 


Reads  in  block  transfer  function  data 
and  checks  it. 


TFBLOK 


Forms  transfer  function  polynomials  in 
frequency  form  for  subsequent  use. 


Checks  polynomials  for  one  rad/sec 
frequency. 


Figure  1  3  Active  Control  Module  Routines 


52 


6C  =  T  ( ico)hs 

where  h$  is  the  amplitude  of  an  inertial  degree  of  freedom  somewhere  on  the 
vehicle.  In  the  modal  solution,  h$  is  composed  of  motion  due  to  aircraft 

t 

rigid  body  perturbations  plus  elastic  motion,  and  in  terms  of  the  generalized 
amplitudes  q: 

hs  =  USJ  (q) 

The  surface  commanded  deflection  at  any  frequency  w  then  must  be: 

{AhE}  =  T(iuj)  {<j>5E)  L<!>SJ  fq> 

The  total  motion  hE^.  of  the  surface  then  is  given  by  the  sum  of  the  elastic, 
rigid  body  and  commanded  motion: 

{hE^  =  +  (<()6E^  Iql 

Since  there  may  be  many  transfer  functions  and  several  control  surfaces  and 
many  sensors,  the  incremental  deflections  may  be  written  as 


mi-  (  . 

{Ah}  =  l  t  . (iw)  J  }  L<t>6j J  tq> 


Note  that  the  modal  amplitude  commanded  by  the  control  system  and  the 
transfer  function  must  be  consistent.  That  is,  if  the  transfer  function  repre¬ 
sentation  is  given  in  terms  of  degrees/unit  sensed  and  the  control  surface 
input  modal  amplitude  is  normalized  to  radian,  then  the  transfer  function 
must  be  scaled  down  by  57.3  degrees/radian. 

In  the  active  system  module  the  data  necessary  to  form  the  augmented 
motion  is  required  in  terms  of  specification,  for  each  transfer  function,  of 
the  data  input  number  of  the  mode  driven,  the  mass  number  corresponding  to 
the  sensed  inertial  degree  of  freedom,  a  code  specifying  which  orientation 
of  that  degree  of  freedom  (x,  y  or  z)  and  a  scalar  multiplier  for  the 
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transfer  function.  This  scalar  multiplier  may  be  used  as  a  direction  cosine 
for  sensed  orientations  different  from  that  of  the  mass  point  or  for  inter¬ 
polation  or  for  general  scaling.  The  transfer  functions  are  calculated  by 
polynomial  multiplication  of  'block'  polynomials  describing  the  active  elements 
of  each  path.  Loops  which  describe  symmetric  control  surface  motion  and 
antisymmetric  motion  are  separately  entered.  A  transfer  function  then  has 
the  form: 


T(1«) 


NN 

I 

0=1 
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l 

n=l 
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The  generalized  mass  and  aerodynamic  matrices  are  augmented  by  the 
active  system  with  the  following  matrices: 


NTF 

[am]  -  I  Tt(1.)  (m,)  1*,^  J 


NTF 

[A-£>(  1u)]  =  l  Tt(iu)  {<&z(ia>)}  jjJs  J 
1  & 


The  integrated  load  matrices  are  augmented  in  a  similar  fashion: 


[aPIQ]  =  l  T^iw)  {PIQ^} 

Jew.)]  - 

_PAQA(iw)J  *=1  1  (  PAQA£(ito)  )  L  l5tJ 


where  m  ,  J)  ,  PIQ  ,  PAQS0,  PAQA  are  the  columns  of  the  matrices  associated 

X.  )t  a>  X  X 

with  the  control  surface  mode  associated  with  the  nth  transfer  function  and 
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tr 


<j>r  is  the  modal  row  giving  the  sensed  degree  of  freedom  motion  associated 
with  the  Uh  transfer  function. 


As  an  example  of  transfer  function  generation,  the  loop  shown  in 
Figure  14  commands  a  control  surface  deflection  due  to  displacement  and 
acceleration  where 

K-j  =  constant  =  0.75 

,  .  («»  *  3s_+  20) 

(s  +  2s  +  5) 

ACT  =  -1500_ 

(S  +  50) 

h  s  +  1 

K  =  (s+18)  (s2+21s+50) 

4  (s+10)  (s2+3s+25) 

u  .  s  +  12 

5  s  +  5 


In  terms  of  6c/h  this  loop  may  be  represented  by  three  simple  transfer 
functions,  each  with  five  'blocks'  as  shown  in  Figure  1  5  where 


k4-1 


s  +  18 
s  +  10 


k4-2 


s2+21s+50 

s2+3s+25 


The  transfer  functions  then  are  (shown  by  block  input) 
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Figure  14.  Sample  High  Gain  Control  Loop 
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TRANSFER  FUNCTION  NO.  1 
ACCELERATION  LOOP 


TRANSFER  FUNCTION  NO.  2 
DISPLACEMENT  LOOP  NO.  1 


v  K  -  K.  -  K  1.0  j - f  ACT 

b  I  .  Z  ■  ! 

i 

TRANSFER  FUNCTION  NO.  3 
DISPLACEMENT  LOOP  NO.  2 


Figure  15.  Sample  Control  System  Transfer  Function  Blocks 


t3(s)  =  (1.0)  (n  7K)/ sjjt  35  +  2o\  / 2500_ 

V  s  +  2s  +  5  /  \s2+1C 


it 


(s^+lOOs+2500/ 

R0UtfW  AFCS  ,0rms  “e  *ls  for  symmetric  and  antisymmetric 

sensors.  TF8LIN  reads  in  the  'block'  transfer  function  data  and  TFBL0K 

converts  it  to  polynomials  in  powers  of  iw.  TRFVAL  prints  not  the  transfer 
faction  polynomials  evaluated  at  one  rad/sec  for  checking. 
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6.  FREQUENCY  RESPONSE  MODULE 


The  frequency  response  module  solves  the  generalized  equations  of  motion 
for  the  generalized  response  to  a  unit  travelling  gust  in  the  frequency  domain. 
The  equations  solved  for  each  required  frequency  and  blast  orientation  are 

{-J-  [m]  -q  O^'(iw)]  +  (1  +  ig)  T  *j)  {q(ioj)}  =  jp  {<9g(iu>)} 

where: 

in  is  the  generalized  mass  matrix 

y,{ iu)  is  the  generalized  aerodynamic  force  matrix  due  to  motion 
*is  the  generalized  stiffness  matrix 
g  is  the  structural  damping  (only  in  the  el  astic  modes) 
q  is  the  dynamic  pressure  (psi) 
q(iw)  is  the  generalized  response 

^^(ic)  is  the  generalized  gust  vector  for  a  specific  orientation 

2  4. 

p  is  the  ambient  density  (lb-sec  /ft  ) 

V  is  the  aircraft  velocity  (fps ) 

In  the  module  the  above  equations  are  solved  for  all  blast  orientations 
available  simultaneously,  passing  automatically  through  the  routines  twice, 
once  for  symmetric  solutions  and  once  for  antisymmetric  solutions. 

Input  data  required  consist  of  the  flight  altitude  and  velocity,  the 
structural  damping  in  each  mode,  an  aerodynamic  tape,  and  the  definition  of 
frequencies  for  solution.  The  input  frequencies  for  solution  should  cover 
the  entire  range  of  given  elastic  modes  of  the  free-free  aircraft  and  out  to 
the  highest  frequency  obtainable  from  the  available  aerodynamic  solutions. 

The  major  routines  in  this  module  are  shown  in  Figure  16. 
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CFREQR 


Dimensions  arrays  and  calls  FREQRS 


I 

FREQRS 


Calculates  the  generalized  response  to 
a  unit  gust 


Calculates  incremental  matrices  for  the 
active  system  contribution  and  adds 
them  to  the  equations  of  motion 


Interpolator  for  motion  dependent 
generalized  aero  from  'hard'  point  aero 
dynamics 


Interpolates  for  unit  generalized  gust 
forces  from  'hard'  point  gust  forces 


Figure  16.  Frequency  Response  Module  Routines 


&0 


Routine  FREQRS  forms  and  solves  the  equations  one  frequency  at  a  time 

and  saves  the  solutions.  Routines  TFMATM,  TFMATF  and  TRFCMX  calculate  the 

incremental  generalized  mass  and  generalized  motion  dependent  aerodynamic 

forces  due  to  the  active  system  feedback,  where  the  incremental  generalized 

mass  matrix  is  given  by: 

NTF 

[Am(ioj)]  =  l  T  (iw)  {m,}  UT  J 
2=1  1  J  1S2 

and  the  incremental  generalized  aerodynamic  matrix  by 

NTF 

[ A  c/j ( i 10 ) ]  =  l  T  ( ia>)  {^(ico)}  L<i>TS  J 
2=1  1  J  1S2 

where 

NTF  =  number  of  transfer  functions  for  symmetric  or  antisymmetric  motion 

T^(iu))  =  the  Jtth  transfer  function  evaluated  at  u>  (in  routine  TRFCMX) 

m.  =  the  j th  column  of  the  generalized  mass  matrix  corresponding 
J 

to  the  jth  mode,  the  driving  mode  (a  control  surface  rotation  mode) 
^j(ia)  =  the  similar  column  of  the  aerodynamic  matrix 

<j>.  =  the  modal  row  giving  the  2th  sensed  degree  of  freedom  motion 

iS2 

in  terms  of  generalized  response 

Routines  TDTINT  and  INTERP  form  the  motion  dependent  generalized 
aerodynamic  forces  at  specific  frequencies  by  spline  interpolation  from  the 
aerodynamic  matrices  at  the  'hard'  points,  where 
NK 

[y. ( i w ) ]  =  l  C.(oj)  O'.(iwj)] 

j*l  J  J 

Routine  FORMC  and  COEFF  form  the  coefficients  C-(u>)  as  function  of 

\J 

the  'hard'  point  reduced  frequencies  and  the  desired  frequency  for  solution. 
Routine  GST INI  interpolates  for  the  generalized  gust  forces,  where 


SPLHP* 

i 

V1"1) 

(  l  U)  ) 

SPLHZ\ 

— 

< 

Fgz(,“>| 

_SPLHY_ 

1 

Fgy<1»)J 

and  the  SFLHP,  SPLH^,  and  SPLHY  are  obtained  from  the  aero  tape  and  are  the 
generalized  forces  due  to  local  box  gust,  body  Z  gust  forces  and  body  Y 
gust  forces  respectively.  Fgp,  F^  and  Fgy  are  the  aerodynamic  element  gust 
forces  which  are  obtained  from  the  'hard'  point  aerodynamic  forces  by 
interpolation: 

♦  r  NK 

Fg  =  e  l  Cj (iw)  RG(1Mj) 
j  =  l 

NK 

4>g  =  l  Cj(iiu)  0G(io>) 

J-l 

where  and  <f>G  are  the  moduli  and  phase  angle  of  the  local  gust  force. 

The  local  gust  forces  for  symmetric  and  antisymmetric  solutions  are 
saved  for  use  in  the  unit  gust  load  module  along  with  the  generalized 
response  solutions  and  the  interpolation  coefficients. 
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7.  UNIT  GUST  LOAD  MODULE 


The  unit  gust  load  module  calculates  the  integrated  loads  and  accelera¬ 
tions  due  to  unit  travelling  gusts  at  specific  orientations.  Input  data 
consist  of  a  unit  load  tape,  an  aerodynamic  tape  and  a  frequency  response 
tape.  The  routines  are  as  shown  in  Figure  17. 


The  routine  CFRL0D  dimensions  the  necessary  arrays  and  calls  FRL0AD. 
FRL0AD  forms  the  symmetric  and  antisymmetric  oscillatory  integrated  loads 
and  accelerations  in  the  frequency  domain  due  to  a  unit  gust,  orientation 
by  orientation.  The  oscillatory  loads  at  a  specific  frequency  and  orien¬ 
tation  are  the  sum  of  the  inertial  loads  and  motion  dependent  aerodynamic 
and  gust  (blast)  aerodynamic  loads: 


Ps(iu.) 

pn(iw) 


[PIQ] 


qs(iu>) 

qfl(iw) 


+  q 


PAQ  S(iw) 
PAQA(Iui) 
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qA(iw) 


+  p_V  [PINTP] 
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Fd  (iw) 


Fp  (iu) 
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,p_v  [PINTY] 
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I  Fy  (iw), 
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F7  (iw) 
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The  accelerations  at  a  specific  frequency  are  given  by: 


*is  0 

a$(iw) 

M 

-°  "  ll\ 

qA(iw) 

Dimensions  arrays  and  calls  FRL0AD 


\ 

i 


Calculates  integrated  loads  due  to 
inertial  and  aerodynamic  motion  dependent 
equilibrium  forces  (inerpolates  on  aero 
forces) 


{  FGMULT 

V _ 


Forms  integrated  loads  due  to  aerodynamic 
gust  forces  for  unit  gust 


Figure  17.  Unit  Gust  Load  Module  Routines 
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The  generalized  response  solutions  qs  and  qA  are  available  from  the 
frequency  response  module.  PIQ,  PINTP,  PINTY,  and  PINTZ  are  frequency 
independent  and  available  from  the  unit  loads  module. 

PAQS  and  PAQA  are  frequency  dependent  and  available  from  the  unit 
loads  module  for  only  specific  frequencies  (the  'hard'  points).  Spline 
interpolation  is  used  to  obtain  these  arrays  at  all  other  frequencies, 
using  the  same  coefficients  which  were  generated  in  the  frequency  response 
module.  That  is 


’PAQS(iw)" 

NK 

=  £ 

C>) 

PAQS  (  i  a;  .  ) 

J 

PAQA(iu)) 

j=l 

j 

PAQA(iu.) 

L  J  J 

The  symmetric  and  antisynmetric  local  gust  forces  Fp(.,  FpA>  Fzs, 

FZA’  FYS’  and  FYA  are  obta^ned  ^rom  *be  frequency  response  module  output. 
Routine  FGMULT  carries  out  the  premultiplication  by  the  appropriate  PINT 
arrays . 

The  active  system  contributions  to  the  loads  are  calculated  by 
augmenting  the  PIQ,  PAQS  and  PAQA  arrays  in  the  same  fashion  as  that  found 
in  the  frequency  response  module. 
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8.  TRIM  MODULE 


The  trim  module  trims  the  aircraft  for  a  specified  flight  condition.  A 
flight  condition  is  defined  as  flight  at  specified  altitude,  speed  (hence  Mach 
number),  maneuver  and  thrust.  The  principal  routines  used  are  CTRIM  and  TRIM, 
where  CTRIM  dimensions  the  arrays  and  calls  the  TRIM  routine  which  carries 
out  the  trim  solution.  CTRIM  is  called  from  the  driving  routine  of  the  blast 
and  time  response  modules.  Input  data  for  this  module  consist  of  an  aero¬ 
dynamic  tape,  a  unit  loads  tape  and  load  factor,  maneuver,  rate  of  climb 
and  engine  thrust: 

n  =  total  load  factor  for  flight  (1.0  for  trimmed  level  flight) 

KMAN  =  maneuver  constant  (0  level  flight,  1  for  pullout,  2  for  turns) 

Z  =  EFAS  rate  of  climb  (or  descent)  in  a  turn  (ft/sec) 

Flight  altitude  and  velocity  for  the  condition  are  obtained  from  the  unit 
gust  load  tape.  Figure  18  details  the  major  routines  found  In  this 
module. 

The  linearized  force  equations  of  motion  for  the  solution  of  trimmed 
flight  loads  are: 

IF)  =  [KFF]  (h)  =  {Fj}  +  (Fa)  +  {Ft> 

where 

Kpp  =  free-free  stiffness  matrix 
Fj  =  inertial  forces 
F^  =  aerodynamic  forces 
Fj  =  thrust  forces 

These  solutions  assume  that  drag  equals  the  forward  component  of  thrust  and 
that  the  aircraft  is  either  in  a  steady  climbing  (or  descending)  turn  or  in 


Dimensions  arrays  and  calls  TRIM 


Calculates  trim  parameters  and  trimmed 
flight  integrated  loads 


Figure  18.  Trim  Module  Routines 


level  flight  or  a  symmetric  pullout. 


The  inertial  force  distribution  for  the  above  restrictions  are  given  by 

'  FI  |  =  n  9  CM] 

where 

n  =  load  factor  (ratio  of  lift  to  weight)  and  the  local  accelerations 
relative  to  rigid  body  oscillations  are  small 
M  is  the  distributed  mass  matrix 

4>j^  is  the  rigid  body  plunge  mode  for  the  inertial  degrees  of  freedom 

The  thrust  forces  are  given  by 
I  FT  \  '  ^XTH^  th 

where:  is  a  matrix  of  direction  cosines  giving  the  components  of  the 

thrust  vectors  aligned  with  the  appropriate  inertial  degrees  of 
freedom  and  was  generated  for  unit  thrust  loads  in  the  unit  load 
module 

Th  is  the  engine  thrusts  for  the  flight  condition  (symmetric  only) 

The  aerodynamic  forces  are  given  by: 

{  fa  }  -  *  5[6]{hJ 

where 

q  is  the  flight  dynamic  pressure 
h^  are  the  aerodynamic  point  deflections 
are  the  aerodynamic  point  velocities 
D,  D  are  the  aerodynamic  AIC's  relating  aerodynamic  force  to  point 
displacement  and  point  velocity. 

The  aerodynamic  point  deflections  are  caused  by  rigid  body  motion, 
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aircraft  jig  shape  (twist,  camber  and  relative  surface  offset  from  the  body 
axes,  but  not  dihedral),  control  surface  position,  and  elastic  deformations. 
The  aerodynamic  point  velocities  are  due  to  the  constant  rotational  velocity 
of  the  vehicle  in  the  prescribed  maneuvers.  Then: 


where 


*AJS,A  = 
*Aa 
% 

'•‘Ae 

^A6e 

♦a«a 

*AES’  %A 


AX,  AZ  = 


jig  modes,  symmetric  and  antisymmetric 
rigid  body  pitch  mode 
rigid  body  yaw  mode 
rigid  body  roll  mode 

symmetric  trim  mode  (elevator,  horizontal  stabilizer,  canard) 
yaw  control  surface  mode  (usually  rudder) 
roll  control  surface  mode  (usually  aileron) 

=  symmetric  and  antisymmetric  aerodynamic  mode  shapes  from 
the  elastic  modes 

distances  from  the  reference  point  for  the  rigid  body  modes 
and  the  c.g.,  and  are  obtained  from  the  appropriate  elements 
of  the  generalized  mass  matrix 


i 


i 
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Antisymmetric  set: 


The  elements  of  the  above  equations  are  obtained  from  the  available  general¬ 
ized  matrices  by  extraction. 

The  generalized  aerodynamics  for  velocity  are  obtained  from  the  lowest 
non-zero  reduced  frequency  (k^)  complex  matrices: 


where  b^  is  the  reference  semichord  from  the  aerodynamic  module. 

The  value  of  n,  P,  Q  and  R  are  obtained  from  input  data  specifying  the 
trim  condition  desired.  For  level  flight  or  a  level  pullout  at  specified 
load  factor: 


P  =  R  =  0 
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where  the  load  factor  n  is  input. 


For  climbing  (or  level  or  descending)  turn: 

Q  =  jp-  0  ~  »  <f>g  =  cos 

n 

ec 

P  =  -  y—  g  tan  <^>B 

R  =  3.  sin 

where  <j>„  is  the  bank  angle  calculated  from  the  input  load  factor  and  e  is 

D  C 

the  climb  angle.  ec  is  calculated  from  input  rate  of  climb  (Z)  and  R  is 
calculated  from  the  bank  angle: 

ec  =  sin_1 

V2 

D  =  1 _ 

g  tan  <j>B 

The  symmetric  and  antisymmetric  integrated  loads  for  the  trimmed  flight 
solution  are  given  by 

J  PSTRIM  J  =  "  9'P'0<*h>!  *5[PAQS(k.0)]|<,s| 

+  5  (  CB_  )  Imag  [PAOSlkJ]  (qs  J  +  [THRWD]  {  T  H  ] 

JPATR1m|  ■  R  [PAQA(k  -  0)]  )qfl[  *q(^)  Imag  [PAQA(k(>]  j  qfl  J 

Additional  aerodynamic  data  are  output  for  the  trim  condition.  These 
are  the  load  distributions  in  varying  degrees  of  detail.  The  aerodynamic 


parameters  and  geometric  data  printed  include  the  following: 

1.  Lifting  pressures  and  their  center  locations  in  each  box  on  all 
lifting  surfaces; 

2.  Lift  coefficient,  moment  coefficient  (about  1/4-chord),  and  center 
of  pressure  on  each  lifting  surface  strip  as  a  function  of  spanwise 
location; 

3.  Spanwise  distribution  of  loading  (in  the  form  c^c/c  on  each  lifting 
surface  as  a  function  of  span; 

4.  Vertical  and  lateral  running  loads  on  slender  bodies  (force  per 
unit  length)  are  output  as  a  function  of  longitudinal  location 
(body  element  centers)  on  all  bodies; 

5.  Total  force  and  moment  coefficients  for  the  entire  aircraft  and 
for  the  lifting  surfaces  separately. 

The  equations  that  are  used  in  the  calculation  of  the  various 

aerodynamic  loadina  coefficients  are  given  in  Section  VI ,  Part  1  of  Vol . 

1 1  of  thi s  report. 
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9.  BLAST  AND  TIME  RESPONSE  MODULES 


These  modules  establish  the  blast  gust  time  and  density  history  at  the 
moving  AAS  origin  of  the  aircraft  and  the  Fourier  transforms  of  the  product 
of  velocity  and  the  local  density,  and  calculate  the  frequency  response  of 
the  aircraft  due  to  the  blast  from  the  unit  gust  frequency  response  and 
acceleration  solutions.  (The  trim  module  is  called  from  the  driving  routine 
CGUSTR).  The  program  then  transforms  these  gust  and  acceleration  reponses 
back  to  the  time  domain.  Summation  of  the  trim  loads  and  the  time  histories 
of  the  perturbation  loads  gives  the  total  loads  on  the  vehicle  as  a  function 
of  time.  The  maximum  positive  and  neqative  loads  are  compared  to  the  input 
allowables  and  the  solution  iterated  for  critical  range  (if  desired). 

Input  data  required  for  this  module  consist  of: 

A  unit  gust  load  tape,  a  unit  load  tape  and: 

N0RMAX  =  the  number  of  orientations  to  be  analyzed 

N0R  =  the  orientation  definition  number  for  each  orientation  analyzed 
REST  =  the  initial  slant  range  estimate  for  each  orientation  analyzed 
EFR  =  burst  yield  (kilotons) 

TIMEMX  =  maximum  time  in  seconds  of  the  response  time  history 
KLPT  =  iteration  control  for  establishing  critical  range: 

0 ,  for  no  i teration 
1 ,  for  i teration 

KGRD  -  ground  reflection  control: 

0,  for  no  reflection 
1,  for  including  reflection 
HGRD  '  height  above  sea  level  of  the  ground  (ft) 

KL0AD  *  allowable  load  modification  constant: 
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0,  use  maximum  allowable  loads  from  unit  loads  tape, 

1,  input  and  modify  allowables 

STALDS(I,J)  =  maximum  positive  (j  =  7)  and  negative  (J  =  8)  allowable 
loads  for  control  range  (entered  only  if  NEWMAX  =  1) 

NCRITS  =  control  constant  for  determining  allowable  range  based  on 
critical  stresses  (is  NSTRSS  not  zero) 

0,  no  critical  stress 

1,  input  and  use  critical  stress  for  range  determination 

ALL0WS(I,J)  =  maximum  positive  (J  =  1 )  and  negative  (J  =  2)  allowable 
stresses  (entered  if  NCRITS  non  zero) 

Gust  orientations  for  the  blast  intercept  and  their  code  numbers  are 
shown  in  Figure  19.  A  side  directed  burst  orientation  may  be  developed  from 
its  symmetric  counterpart  on  the  opposite  side,  hence  the  need  for  fewer 
asymmetric  orientations  in  the  aero  and  response  modules.  A  tabular  corres¬ 
pondence  between  the  orientation  numbers  of  the  aero  module  and  the  blast 
module  is  used  to  take  advantage  of  this.  Table  3  gives  the  standard  direction 
cosines  for  the  standard  orientations. 

The  substantial  difference  in  this  analysis  procedure  and  that  of  the 
VIBRA-4  program  is  seen  to  be  the  use  of  Fourier  transforms  to  and  from  the 
frequency  domain  to  establish  the  time  histories  of  loads  rather  than  piece- 
wise  integration  in  the  time  domain.  This  procedure  is  necessary  because  the 
complex  three  dimensional  unsteady  aerodynamic  solutions  are  known  only  for 
constant  amplitude  oscillatory  motion.  The  Fourier  transform  process  properly 
accounts  for  the  time  lags  of  the  blast  induced  gust  impingement  across  the 
aircraft  spatially  in  that,  in  the  limit,  every  possible  phase  lag  across 
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Figure  19.  Blast  Orientations  for  Mast  Module 
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TABLE  3 


DIRECTION  COSINES  FOR  STANDARD  ORIENTATIONS 


ORIENTATION 

NUMBER 

lx 

Lc 

Iz 

1 

0.8660254 

0.0 

0.50 

2 

0.50 

c.o 

0.8660254 

3 

0 

0.0 

1.0 

4 

-0.7071068 

0.0 

0.7071068 

5 

-0.7071068 

0.0 

-0.7071068 

6 

0 

0.0 

-1.0 

7 

0.50 

0.0 

-0.8660254 

8 

0.8660254 

0.0 

-0.50 

9 

0.7071068 

-0.7071068 

0.0 

10 

0 

-1.0 

0.0 

11 

-0.7071068 

-0.7071068 

0.0 

12 

-0.7071068 

0.7071068 

0.0 

13 

0 

1.0 

0.0 

14 

0.7071068 

0.7071068 

0.0 

15 

0.0 

0.7071068 

0.7071068 

16 

0.0 

-0.7071068 

0.7071068 

17 

0.0 

-0.7071068 

-0.7071068 

18 

0.0 

0.7071068 

-0.7071068 
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k  =  reduced  frequency 
M  =  free  stream  Mach  number 
b  =  reference  semi  chord 

yx  =  direction  cosine  of  the  blast  wave  front  with  respect  to  the 

XAAS  axis 
t  =  shifted  time 


i  =  y  x  +  y  y  +  Y  z 

where 

x  =  X^s  coordinate  of  the  point 

y  =  coordinate  of  the  point 

z  =  coordinate  of  the  point 

y  ,  yz  =  direction  cosines  of  the  blast  wave  with  respect  to  the 

YAAS  and  ZAAS  axes>  respectively. 

The  shifted  time  t  is  used  since  the  gust  downwash  distribution  over  the 
vehicle  is  based  on  distance  from  the  origin  of  the  AAS  system.  Since  t, 
true  time  must  be  zero  at  time  of  wave  interception  with  the  aircraft,  shifted 
time  is  given  by 
t  =  t  +  T 

where  t  is  true  time  and  T  is  the  time  for  the  wave  to  move  from  interception 
point  to  the  origin  of  the  AAS  system,  and  is  given  by 
Yx  X,  ♦  Yv  V,  *  Y;  Z, 

? - ?  ~x 


where 

X j ,  Yj,  Zj  =  AAS  coordinates  of  the  point  of  shock  interceiption 
V$s  =  speed  of  sound. 
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The  inverse  Fourier  transform  of  the  frequency  response  of  a  stable  aircraft 
to  gust  forces  resulting  from  the  above  downwash  distributions  may  be  shown 
to  be  the  time  history  of  the  aircraft  subjected  to  a  unit  impulsive  gust. 

The  Fourier  transform  of  the  gust  time  history  convoluted  in  the  frequency 
domain  with  the  aircraft  frequency  response  function  and  subsequently 
inverse  transformed  will  be  the  time  response  of  the  vehicle  to  the  actual 
gust  time  history.  The  Fourier  transformation  requires  that  gust  function, 
G(t),  be  absolutely  integrable;  and  such  is  the  case.  Examination  of  typical 
material  velocity  time  histories  shows  that  their  Fourier  transforms  are 
effectively  zero  in  the  range  of  40  to  60  Hz.  This  implies  that  the  inverse 
Fourier  procedure  will  predict  the  time  history  of  the  aircraft  well  if  all 
roots  of  the  aeroelastic  system  out  to  the  maximum  frequency  necessary  to 
transform  the  gust  function  are  available  in  the  aircraft  frequency  response 
function.  The  implication  of  this  is  that  either  the  frequency  domain  aero¬ 
dynamics  must  be  known  out  to  this  frequency  range  or  the  modal  set  describing 
the  aircraft  vibration  characteristics  must  be  truncated  at  sufficiently  low 
frequency  such  that  no  roots  exist  between  the  last  known  frequency  response 
solution  frequency  and  the  frequency  of  zero  moduli  (uiMAX). 

The  driving  routine  GUSTDR  calls  the  time  module  for  the  specified 
maneuver,  forming  the  symmetric  and  antisymmetric  trimmed  flight  loads  P$TRIM 
and  pATRjM  and  then  loops  on  specified  orientations  for  the  resulting  load 
time  histories.  Subroutine  GSTHST  calculates  the  burst  scaling  factors, 
gets  the  overpressure  ratio  from  PRESS  and  the  gust  velocity  from  TPEVAL  and 
GUSHR0  as  a  function  of  time  and  range.  Shock  time  of  arrival  is  calculated 
in  TAR.  Range  is  updated  each  time  step  in  routine  FLTPOS  assuming  the  mean 
flight  path  reasonably  follows  the  trajectory  dictated  by  the  maneuver  speci- 
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fication  (climbing  turn,  pullup  or  level  flight).  Routines  TPINT  and  INTRBL 
are  used  to  include  the  gound  reflected  and  incident  shock  waves.  ROUTINES 
TPEVAL,  GUSRH0,  TAR,  PRESS,  INTRPB  and  TPINT  are  taken  from  Reference  1.  The 
EFAS  burst  coordinates  at  time  of  burst  are  given  by 


where 

SQ  is  the  slant  range  at  time  zero  (intercept  time) 

e  is  the  climb  angle 
c 

<1>B  is  the  bank  angle 

y  y  y  are  direction  cosines  of  the  burst  orientation  with  respect 
Tx*  'y  z 

to  the  AAS  system 

RT  and  hQ  are  the  turn  radius  and  altitude  at  time  zero. 


The  aircraft  position  in  the  EFAS  at  time  t  is  given  by 


i;ei 

iej 


/  -Rj  sin  Rt 

;  RT  cos  Rt 
T 

v  h  +  (Vsine  )t 
0  ^ 


climbing  turn 


where 

hQ  =  flight  altitude 

o  - 

v 

R  -  ^ 

KT  0 

V  =  aircraft  flight  velocity  (fps) 

t  =  time  (secs) 


1 
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(X  \  (  -R„  sin  Q. 

EM  P„  ‘ 

I  E  )  0 

*  Z-  '  I  k  +d  n.c 


where 


ho+Rp[l-cos(Qt)j 


6  -  (2il)g 
4  v 

v2 

rp  “  X^Dg 

n  =  load  factor 


symmetric  pullout 


)Xe)  rvti 

iir  U 


for  level  flight 


The  slant  range  S  at  time  t  is  given  by 

s(t)  =  y  (xE  -  xB)2  +  (ye  -  yb)2  +  (XE  -  zB)2 

The  time  intervals  are  taken  as  0.05  seconds  until  the  gust  velocity  goes 
negative,  then  at  0.50  second  intervals  till  the  gust  velocity  is  less  than 
2  fps  (time  T^^).  An  extrapolation  is  the  i  made  to  find  T^^,  the  time 
at  which  the  gust  velocity  has  returned  to  zero. 


Routine  TIMHST  takes  the  output  time  history  of  the  gust  velocity 
and  local  density  p  and  with  routine  TRFFT  Fourier  transforms  the  product 
into  the  frequency  domain: 

G(t)  =  p(t)  •  Vg(t)/paiJbjent 


I 

and  g(iw)  =  J 


G(t)e'1'ut  dt 
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It  is  convenient  to  use  an  analytic  function  for  G(t)  for  which  the  exact 
Fourier  transform  is  known.  The  approximation 
G(t)  =  GQ(2e‘a(t+T/  -  e'0(t+T)) 

where  GQ  is  the  gust  velocity  at  time  t  =  -t  .  fits  the  G  function  very 
well  for  numerous  time  histories  of  gust  velocity  with  respect  to  the  moving 
airframe.  This  is  generally  the  case  when  the  vehicle  is  moving  towards 
the  burst.  The  Fourier  transform  of  G(t)  then  is  given  by: 


g(iw) 


+i 


)(¥ T 

2  2 

COSwT  +  uj  / 

(|_\a  +u 

e  +w 

l 

ITT  "  Xl'l  sinwT  '  w  Vrf  ' 

[a  0  +oj  J  \a  +w 


—x — x-  si  nun 

6+oj/ 


COSlDT 


) 

) 


Routine  TRFFT  checks  the  fit  of  the  function  G(t)  in  the  region  of 
time  when  the  gust  velocity  is  negative  (which  is  the  area  of  poorest  fit) 
and  if  the  fitting  function  predictions  exceed  five  percent  of  the  actual 
time  histories  the  Fourier  transform  is  taken  as  the  sum  of  a  series  of 
impulses  over  the  time  history  of  the  burst. 


The  gust  and  density  function  is  then  convoluted  in  the  frequency 
domain  with  the  syinnetrlc  (P^)  and  antisymmetric  (P^)  load  frequency 
response  to  the  unit  gust: 

Pg$(iw)  =  PS(ia))  •  9(iw) 

pgA(iu)  =  PA(io))  ’  g(iw) 
to  form  the  perturbation  loads  for  the  blast. 
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The  time  history  of  the  perturbation  loads  are  calculated  in  I  FT : 


V*)  =  -i- 


Imag  [PgS(  iw)  sino)t]dw 


=  "tt 


Imag  [Pg^(ioj)  sinwtjdcu 


A  discussion  of  the  method  of  Fourier  transforms  used  herein  may  be 
found  in  References  7,  8  and  9. 


Routine  LOADCH  forms  the  right  side  and  left  side  time  history  loads 
generated  from  the  half  aircraft  analysis: 

PRHS(t)  =  PSTRIM+  PATRIM+  ps(t)  +  P,(t) 

%S(t)  =  PSTRIM  "  PATR I M  +PS ( ^  ~  PA(t) 

Maximum  positive  and  negative  loads  are  saved  and  compared  to  allowable 
positive  and  negative  loads.  If  a  stress  matrix  has  been  defined  in  the 
limit  loads  module,  the  stress  time  histories  are  found  and,  if  so  flagged, 
the  maximum  allowable  gust  velocity  established  on  the  basis  of  input 
allowable  stresses.  The  maximum  allowable  gust  velocity  for  the  burst 
conditions  are  estimated  from 


7.  Wylie,  C.R.  Jr.,  Advanced  Engd.nceAd.ng  Mathematics ,  McGraw-Hill  Book 
Co.,  Inc.,  New  York,  1951. 

8.  Solodovnikov,  V.  V.,  IntAoduction  to  the  Statistical  Dynamics  of, 
Au-toma-ted  Contact  Sybtemi,  Dover  Publications,  Inc.,  New  York,  1960. 

9.  Hurty,  W.C.  and  Rubinstein,  M.  F.,  Dynamic, 5  oft  StAucTuAci, Prentice- 
Hall,  Inc.,  New  Jersey,  1964 
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*  MAX(t  *+0>  ■  W‘>  '  V  ■t0> 

9  _ 2_ 

^A1 low 

and  the  maximum  allowable  overpressure  from: 


where  is  the  ambient  SDeed  of  sound  (FPS)  and  dq  is  the  ambient  pres¬ 
sure  (psi). 

A  revised  slant  range  is  calculated  by  routine  RANGE,  which  is  an 
inversion  of  routine  PRESS.  If  iteration  for  range  has  been  specified, 
the  time  history  load  calculation  process  is  repeated  until  convergence 
on  an  allowable  gust  velocity.  At  the  conclusion  of  a  response  solution, 
the  distance  from  burst  to  intercept  (in  the  EFAS)  is  displayed. 

Figure  20  details  the  major  routines  in  this  module. 


85 


CGUST 


Dimensions  arrays  and  calls  GUSTDR 
and  CTRIM 


GUSTDR 


y 


Driving  program  for  Trim,  Blast  and 
Time  History  modules 


CTRIM 

TRIM 


Calculates  trimmed  flight  conditions 
and  steady  loads  (Trim  module) 


''  GSTHST 


Calculates  time  history  of  gust  and 
density  as  a  function  of  spatial  orient¬ 
ation  with  respect  to  blast  point 
(Blast  module) 


(TPEVAL 


J 


Sets  up  calculation  for  gust  velocity 
calculations,  including  ground  reflection 


"PRESS  '  ^Calculates  overpressure,  given  range 

TPINT  ^Calculates  height  and  range  of  triple 
point  path  interception 

GUSRHO  Calculates  gust  velocity  and  local 
/  density  from  given  overpressure 


z  '  IKTRPB  Interpolates  as  required 


TAR 


,y  Calculates  shock  arrival  time 


/  FLTPOS 


/timhst 


Calculates  range  from  burst  to  aircraft 
AAS  origin  at  any  time  t 


Performs  Fourier  transform  of  gust  time 
density,  forms  frequency  domain  solution 
for  blast  and  calculates  inverse  transfonn 
time  histories  of  symmetric  and  anti- 
symnetric  solutions.  (Time  history  module) 


CTRFFT  \  Fourier  transform  routine 
IFT  Inverse  Fourier  transform 


* 

Figure  20.  Blast  and  Time  Response  Module  Routines 
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1 

AD-/1106  480  DOUGLAS 

f  NUCLEAR 

AUG  81 

UNCLASSIFIED 

AIRCRAFT  CO  LON*  (EACH  CA 

BLAST  RESPONSE  COMPUTER  PROMAN.  VOLUME 
J  A  MCBREVr  J  P  SIESINSi  T  P  KALMAN 

AFML-TR-B1-3*- 

F/G  18/3 

I*  PROGRAM  DESC— ETCClJI  fl 

DNA001-73-C-0216  1 

VOL-1  NL  1 

tm 

j 

■■ 

WI^M 

■■ 

■■■ 

^■1 

■ 

jj 

■ _ 

i  LOADCH 

V _ > 


Calculates  left  and  right  side  pertur¬ 
bation  loads  as  a  function  of  specific 
orientation  and  sums  in  trim  loads. 
Maximum  and  minimum  loads  are  found  and 
maximum  overpressure  for  allowables 
calculated  for  next  iteration,  if  any. 


STRSCH 


Calculates  left  and  right  side  total 
stresses  versus  time. 


Calculates  new  allowable  range  for  air¬ 
craft  for  specified  overpressure  and 
gust  velocity. 

Finds  distance  from  aircraft  to  burst 
at  burst  time 


Figure  20  (contd).  Blast  and  Time  Response  Module  Routines 
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10.  RIGID  MODULE 


The  RIGID  module  has  been  added  to  permit  correlation  studies  with 
experimental  measurements  of  the  blast  loading  of  (effectively)  rigid  models. 
The  experimental  data  on  time  histories  of  the  moving  reference  point 
overpressure,  density,  and  material  velocity  are  input  by  the  user  and  are 
utilized  in  place  of  the  present  internal  calculations.  The  correlation 
pressure  points  are  also  specified  by  the  user  in  terms  of  the  locations  of 
the  aerodynamic  boxes  and  slender  body  elements  so  that  the  pressure  time 
histories  for  a  steady  angle  of  attack  of  the  rigid  vehicle  can  be  determined. 

The  module  first  calculates  the  trim  forces  for  the  specified  aerodynamic 
points  by  multiplying  the  forces  from  the  AERO  file  by  the  dynamic  pressure 
and  the  input  values  of  generalized  coordinates  for  any  modes  input  to  the 
AERO  module.  Next,  the  module  performs  a  Fourier  transform  on  the  product 
of  the  input  material  velocity  and  density  in  order  to  put  it  in  the  frequency 
domain.  After  multiplying  the  result  by  the  dynamic  pressure  and  the  gust 
forces  from  the  AERO  file,  an  inverse  Fourier  transform  is  performed  to 
obtain  the  forces  in  the  time  domain.  Adding  these  to  the  trim  forces  and 
dividing  by  box  areas  or  body  lengths  leads  to  the  desired  pressures  and 
running  loads. 

The  output  from  this  module  is  then  the  time  histories  of  the  pressures 
and  running  loads  at  the  specified  aerodynamic  boxes  and  slender  body 
elements.  These  may  be  compared  to  corresponding  experimental  results. 

Figure  21  details  the  major  routines  in  this  module. 
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Dimensions  arrays  and  calls  RIGID 


Calculates  rigid  loads 


Calculates  blast  intercept 


Interpolates  gust  forces  from  'hard' 
point  gust  forces 


Fourier  transform  routine 


Inverse  Fourier  transform  routine 


Figure  21.  Rigid  Module  Routines 


11.  MERGE  MODULE 


The  MERG  module  has  been  added  to  permit  the  merging  onto  one  file, 
the  aerodynamic  matrices  for  different  reduced  frequencies  from  two 
separate  files.  This  module  adds  a  capability  that  will  save  computing  time 
in  that  a  set  of  aerodynamic  matrices  for  a  new  set  of  reduced  frequencies 
can  be  obtained  without  rerunning  any  reduced  frequencies  that  have  already 
been  previously  obtained.  This  feature  also  saves  computing  when  an  AERO 
job  runs  out  of  time  before  computing  all  of  the  reduced  frequencies  requested. 
The  user  simply  runs  the  AERO  module  for  the  frequencies  not  completed  and 
merges  with  the  previous  file  to  obtain  the  complete  set  of  aerodynamic 
matrices. 

Onlv  two  files  can  be  input  for  each  run,  but  the  output  file  can  always 
be  input  to  a  subsequent  merge  run.  The  filenames  of  the  input  files  are 
TAPE  17  and  TAPE  18,  and  the  file  name  of  the  output  file  is  TAPE  19. 
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SECTION  IV 


PFOGFAM  INPUT 

1.  GENERAL  DESCRIPTION  OF  INPUT  DATA 

The  program  input  is  divided  into  two  (2)  card  decks  called 
the  Fixed  Data  Deck  and  the  Fun  Data  Deck  respectively.  The  Fixed 
Data  Deck  is  divided  into  four  (4)  major  q roups,  according  to  the 
program  module  which  uses  the  data,  and  contains  all  of  the 
information  which  define  the  aircraft  geometry,  mass,  mode  shapes, 
frequencies,  loads,  etc.  The  Run  Data  Deck  contains  the 
information  necessary  to  define  the  conditions  describing  the 
specific  run;  for  example,  the  blast  orientations  to  be  considered 
and  the  associated  trial  ranges. 

Each  input  card  consists  of  either  all  integer  numbers  or  all 
real  numbers  input  in  fields  of  12  columns.  Integer  numbers  are 
input  right  justified  in  the  field;  that  is,  the  number  is  punched 
in  the  card  so  that  it  ends  in  the  right-most  column  of  the  field. 
Real  numbers  are  generally  left- justified,  which  means  the  number 
is  punched  in  the  card  so  that  it  starts  in  the  left-mcst  column 
of  the  field.  In  the  input  descriptions  that  follow,  the  ending 
card  column  for  integer  numbers  or  the  beginning  card  column  for 
real  numbers  are  indicated  above  the  dashed  lines  and  the  input 
variables  are  indicated  below  the  dashed  line.  Each  item  number 
begins  on  a  new  card. 

Units  of  data  describing  geometry  must  be  consistent  and  are 
expected  in  inches,  thouqh  an  input  parameter  SIZFCT  is  provided 
in  the  Run  Data  Deck  to  adjust  for  geometric  data  in  units  ether 
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than  inches.  Units  for  the  modal  deflections  are  expected  to  be  inches  for 
linear  deflections  and  radians  for  angular  deflections.  These  linear  or 
angular  deflections  are  given  for  a  unit  modal  amplitude.  Jig  mode  input  for 
calculating  the  basic  lift  and  moment  forces  on  the  undeformed  vehicle  must 
be  in  the  same  units  used  for  the  geometry.  Modes  which  describe  the  aero¬ 
dynamic  and  inertial  deflections  in  terms  of  rotations,  such  as  rigid  body 
pitch  or  control  surface  rotation  modes,  should  be  for  one  radian  of  rotation. 
An  input  parameter,  RBRADF,  is  provided  in  the  trim  module  inout  in  the  event 
some  other  rotational  base  value  has  been  used.  All  such  rotation  modes  must 
have  the  same  reference  rotational  value. 

All  other  unit  requirements  are  as  described  in  the  detailed  input 
data  section. 
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2. 


FILE  CONTROL  STATEMENTS 


All  input  and/or  output  data  files  used  by  VIBRA-6  are 
sequential  files  arid  can  be  stored  on  tape  or  disk  directly  at  the 
time  of  creation  using  standard  CDC  control  statements.  The 
control  statements  for  saving  a  file  are  the  REQUEST  statement  for 
saving  on  tape  and  the  CATALOG  statement  for  saving  on  disk.  The 
control  statements  for  using  a  previously  saved  file  are  the 
REQUEST  statement  for  a  file  saved  on  tape  and  the  ATTACH 
statement  for  a  file  saved  on  disk.  Files  may  be  copied  from  tape 
to  disk  or  disk  to  tape  by  the  use  of  the  COPYBF  control  statement 
without  affectinq  the  operation  of  the  program.  For  the  proper 
coding  of  control  statements,  the  appropriate  operating  system 
manual  should  be  consulted. 


The  following  correspondence  must  be  maintained  between  the 
data  files,  the  fortran  unit  designator,  and  the  CDC  filename. 


Fo£trip._UciT 


Run  Lata 

5 

Fixed  Data 

31 

Aerodynanics 

19 

Unit  Loads 

34 

Frequency  Response 

35 

Unit  Gust  ^oads 

36 

TAPE5 

TAPE31 

TAPE19 

TAPE34 

TAPE35 

TAPE36 
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3 


FIXED  DATA  DECK  INPOT  DESCFIPTION 


The  four  major  groups  in  the  Fixed  Data  Deck  are  the 
sectional  input  data,  the  aerodynamic  module  input,  the  inertial 
module  input,  and  the  load  module  input.  Each  group  is  described 
in  more  detail  preceeding  the  specific  input  descriptions  that 
follow.  If  a  qiven  module  is  not  to  be  executed  in  a  run,  the 
input  data  for  that  module  does  not  need  to  be  removed  from  the 
Fixed  Data  Deck.  The  program  can  find  the  data  necessary  to  run 
any  specific  module  via  the  data  designator  card  which  is  the 
first  card  in  each  major  group.  The  Fixed  Data  Deck  is  input  usinq 
fortran  unit  31  (TAPE31). 


a.  Sectional  Input  Data  (SECT) 


The  input  data  for  this  module  consist  of  flags  for 
dimensioning  and  module  control,  locations  of  and  deflections  at 
the  inertial  nodal  points,  and  modal  frequencies  and  structural 
dairpinq.  This  data  is  optional,  but  if  it  is  input,  it  must  be 
the  first  lata  in  the  Fixed  Data  Deck.  Either  the  SECT  or  the 
IMCT  data  must  be  input  and  the  program  will  use  the  first  one 
that  appears  in  the  Fixed  Data  Deck. 


ITEM  NJ. _ 1  Data  designator  card 

cc  1  2  3  4  6 

iic? 


Jfariablg  Description 

SECT  =SECT ,  designates  that  the  data  following  are  the 

inertial  module  input  data 


ITEM_NJ^2 


Dimensioning  data 


cc 

1 

2 

3 

4 

6 

7 

1 

_ 2__ 

_ 4 _ 

_ 2 _ 

_ Q _ 

_ i 

NMS 

NDCF 

NSYM 

NASYM 

NENGS 

KPRLDS 

Variable 


Inscription 


NMS 


NDOF 
NSYM 
NAS  YM 
NENGS 
KPF IDS 


number  of  'right  hand  side'  mass  points,  including 
centerline  masses 

number  of  degrees  of  freedom  per  mass  station 
number  cf  symmetric  modes 
number  of  antisymmetric  modes 

number  cf  engines  on  right  hand  side  and  centerline 
*1,  print  unit  load  matrices  generated  (if  any) 

=0,  do  not  print 


1 


I 
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ITEM  NJ, 

.2 

Mass  station  locat 
stations,  1=1, NMS) 

ions. 

(Pepeated 

for  all  irass 

cc 

1 

2 

3 

4 

6 

1 

3 

5 

7 

9 

1 

EIXIO 

ELYIO  ELZIO 

Variable 

2es£riBticn 

SLXIO 

XAAS  coordinate 

of 

mass. 

in. 

FLYIO 

YAAS  coordinate 

of 

mass. 

in. 

SL2IO 

ZAAS  coordinate 

of 

mass. 

in. 

ITEM  NJ, 

4 

Inertial  modal 

f requenci 

es,  FREQ,  and  mode  shapes 

PHIX ,  PHI Y ,  and 

PHIZ.  (Pepeated  for 

all  nodes. 

N=1,NCMCE,  whore  NCMOD=NSYM-fNASYM) 

cc 

1 

2 

3 

4 

6 

1 

3 

5 

7 

9 

1 

FPEQ 

FREQ  FRED 

FRED 

FPEQ 

FPEQ 

Variable  Description 

FRED  modal  frequency,  Hz 


ITEM  !i),  5 

Modeshaf  es 

(Pepeated  for  all  mass  stations 

,  1  =  1 

NMS  and  all 

modes,  N=1  to  NSYM+NASYM) 

FIRST  CAPD 

cc 

1 

2 

3 

4 

6 

1 

3 

5 

7 

9 

1 

F 

T_, 

H 

THETA 

ALPHA 

PSI 

Variable 

Description 

F 

linear  deflection  in  the  x- 

direction 

*4 

linear  deflection  in  the  y- 

direction 

H 

linear  deflection  in  the  z- 

direction 

THETA 

rotational 

deflection  about 

the  x-axis 

ALPHA 

rotational 

deflection  about 

the  y-axis 

PSI 

rotational 

deflection  about 

the  z-axis 
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SECOND  CARD  (Input  only  if  NDOF=7  or  8) 

cc  1  2  3  4  6 

1 _ 3 _ 5 _ 7 _ 9 _ 1 

Seta  delta 


Variable 


Description 


PETA 


DELTA 


rotational  deflection  abcut  the 
surface  hinge-line 
rotational  deflection  about  the 
surface  hinge-line 


primary  control 
secondary  control 


ITLM_N  3, _ 6  Mass  property  table  (Repeated  for  all  mass  stations, 

1=1, NMS),  Note  that  this  is  not  the  mass  matrix  that 
is  used  by  the  program,  tut  a  table  of  mass 
properties  that  is  used  to  form  the  mass  matrix. 

£I£ST_CAFJD 


cc 

1  2 

3 

4 

6 

1 

3  5 

7 

9 

_  1 

M 

MDELX  MDELY 

MDELZ 

Variable 

Description 

M 

mass 

MPELX 

mass  unbalance. 

x-direction 

MDELY 

mass  unbalance. 

y-direction 

MDELZ 

mass  unbalance. 

z-direction 

SE£OND_cafd 

cc  1 

2 

3  4 

6 

1-  —  3 

5 

7  9 

1 

IXX  IYY 

IZZ 

I X  Y  I YZ 

IZX 

T  X  X 
I  Y  Y 
IZZ 
I  X  Y 
I  YZ 
IZX 


moment  cf 
moment  cf 
moment  of 
product  of 
product  of 
product  of 


inertia  about  x-axis 
inertia  about  y-axis 
inertia  about  z-axis 
inertia  about  x-  6 
inertia  about  y-  6 
inertia  about  z-  £ 


y-axis 

z-axis 

x-axis 


! 


i 
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lillEQ-SiiEfi  (Input  only  if  NDOF=7  or  8) 


cc  1  2  3  4  6 

1 _ 3 _ 5 _ 2 _ £ _ 1 _ 

SHE  PAB  IBB  PTE  SFB  ~  PPB 


V  ariable 


D§SSEiEii£JD 


SHE 
PAB 
IBE 
PTr 
S  Ft 
PPF 


mass  unbalance  for  primary  control  surface 
product  of  inertia  for  primary  control  surface 
moment  cf  inertia  for  primary  control  surface 
product  of  inertia  relating  theta  and  beta 
mass  untalance  relating  f  and  beta 
product  of  inertia  relating  psi  and  beta 


FOUFTK_CAFD  (Input  only  if  ND0F=8) 

cc  1  2  3  4  6 

1 _ 3 _ 5 _ 2 _ 9 _ 1 

SHE  ~  PAD  ~  PBD  ~  ICE 


Variable 


Description 


SHD 

p;,r 


PED 

rrr 


mass  unt.alance  for  secondary  control  surface 
product  of  inertia  for  secondary  control  surface 
relating  alpha  and  delta 

product  of  inertia  for  secondary  control  surface 
relating  teta  and  delta 

moment  cf  inertia  for  secondary  control  surface 


EIITlLSSiSB 


(Input  only  if  ND0F=8) 


.  f  1  2  3  4  6 

1.  _ 3 _ 5 _ 2 _ 9 _ 1 

PT.  SFC  PPD 

V 

Vdri,.  Description 


PTC 

SFD 

PPD 


product  of  inertia  for  secondary  control  surface 
relating  theta  and  delta 

mass  untalance  for  secondary  control  surface  relating 
'  f  and  delta 

\  product  of  inertia  for  secondary  control  surface 
'^relating  psi  and  delta 

V 


v 
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IIEM_NQj._2  Modal  definitions,  AMODNO. 
FI£ST_CAFD 


1 

2 

3 

4 

6 

2 

4 

6 

8 

_ £ 

M0DE1 

MODE  2 

M0DE3 

MODE4 

MODE  5 

Variable  pescription 

MOT  El  mode  number  of  rigid  body  plunge  mode 

MCDE2  mode  number  of  rigid  body  pitch  mode 

MODE 3  mode  number  of  rigid  body  fore  and  aft  mode 

M0EE4  mode  number  of  first  symmetric  elastic  mode 

M0CE5  mode  number  of  last  symmetric  elastic  mode 


SECOND  CAPD 

cc  1  2 

1 _ 2 _ 4 

MOCE6  "*  M0CE7 


3  4  6  7 

_ 6 _ £ _ £ _ 2 

MODE8  MODE9 


Variable  Description 

M0LE6  mode  nunber  of  pitch  trim  mode 

MODE7  mode  number  of  symmetric  jig  mode 

MOD 38  mode  number  of  first  symmetric  mode  to  be  deleted  in 

this  analysis 

MODE9  mode  number  of  last  symmetric  mode  to  be  deleted  in 

this  analysis 


IHIE£_c£FD 


cc 

1 

2 

3 

4 

6 

1 

2 

4 

_ § _ 

_  _  £  . 

£. 

M0DE11 

MODE  12 

M0DE13 

MO  DEI  4 

MODE 1 5 

idtidble  Description 

M0DE1 1 
MOD El  2 
MOT El  3 
MOD314 
MOrE15 


mode  number  of  rigid  body  roll  mode 

mode  number  of  rigid  body  yaw  mode 

mode  number  of  rigid  body  lateral  mode 

mode  number  of  first  antisymmetric  elastic  mode 

mode  number  of  last  antisymmetric  elastic  mode 
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F0UFTH_£AE2 


CC  1  2  3  4  6 

kodeT6  mod  e 17  mode! I  mcdItI  MODE  2  0 


Variable 


Description 


MODE16 
MODE17 
MOTEIS 
MOD  519 

MODE20 


Note : 


mode  number  of  roll  trim  mode 

mode  number  of  yaw  trim  mode 

mode  number  of  antisymmetric  jig  mode 

mode  number  of  first  antisymmetric  mode  to  be  deleted 

in  this  analysis 

mode  number  of  last  antisymmetric  mode  to  be  deleted 
in  this  analysis 

Trim  modes  may  also  be  defined  as  elastic  modes.  For 
example,  the  yaw  trim  mode  (usually  rudder  rotation) 
might  also  be  used  in  the  yaw  damper  system  and  hence 
be  defined  as  an  elastic  mode,  a  trim  mode  and  in  the 
ACS  definition. 


ITEM  NO. _ 8  Modal  structural  damping,  CPXDPG. 

(Fepeated  until  the  damping  for  all  modes  has  been 
read  in  for  N=1,NOMOD) 


cc 

1 

2 

3 

4 

6 

1 

3 

_ 5 _ 

7 

9 

_ 1 _ 

CPXDPG 

CPXDPG 

CPXDPG 

CPXDPG 

CPXDPG 

CPXDPG 

Variable  Description 

CPXDPG  modal  structural  dampinq 

The  structural  damping  should  be  input  zero 
in  all  modes  except  the  elastic  modes. 


Note : 


b.  Aerodynamic  Module  Input  Data  (AERO) 


The  input  data  for  this  module  consist  of  the  following  five 
groups  of  data:  general  data,  panel  data,  body  data,  modal  spline 
interpolation  data,  and  qust  data.  The  general  data,  panel  data, 
and  irodal  spline  interpolation  data  must  always  be  input,  but  the 
body  data  and  qust  data  are  optional  as  specified  by  flags  in  the 
general  data. 

The  general  data  consist  of  flags  for  dimensioning  and 
program  control,  locations  of  and  deflections  at  the  aerodynamic 
nodal  points,  and  constants.  The  panel  and  body  data  consist  of 
the  input  necessary  to  define  completely  each  panel  and  body.  The 
modal  spline  interpolation  data  consist  of  the  data  to  relate  the 
aerodynamic  nodal  points  to  the  panel  and  body  data.  The  gust 
data  allow  nonstandard  blast  orientations  to  be  input  for  the  gust 
calcu lations. 

Units  of  data  describing  geometry  must  be  consistent  with  the 
units  for  the  other  major  qroups. 

GFCUP  NO.  1  General  Data 

ITEM  NO, _ 2  Data  designator  card 

cc  1  2  3  U  6 

1 _ 3 _ 5 _ 7 _ 9 _ ] _ 

AEFO 

Variable  Descrj.ptj.on 

AEFO  =AEFO ,  designates  that  the  data  following  are  the 

aerodynamic  module  input  data 


101 


ITEM  NO. _ g  Dimensioning  data 


cc 

1 

2 

3 

4 

6 

7 

1 

2 _ 

4 

6 

8 

_  -0  .  __ 

_  .2 

NODES 

NSYM 

NASYM 

mfixT 

MFIX2 

ld£i£bl£ 

Description 

NODES 

number  of  aerodynamic  nodal  points 

NS  YM 

number  of  symmetric  modes 

NASYM 

number  of  antisymmetric  modes 

MFIX1 

mode  number  of  first  mode 

to  monitor 

MFIX2 

mode  number  of  second  mode 

to  monitor 

Note: 

Gusts  can  not  be  monitored 

• 

ITEM  NO. 

3 

Dimensioning 

data 

cc 

1 

2 

3 

4 

6 

7 

1 

2 

4 

6 

8 

0 

_2 

NP 

MSTRIP 

NSMAX 

NCMAX  NBOXES 

Variable 

Description 

NP 

Total  number 

of  panels  on 

all  lifting 

surfaces 

MSTPIP 

Total  number 

of  strips  for 

all  panels 

NSMAX 

Maximum  number  of  strips  per  panel 

NCMAX 

Maximum  number  of  chordwise  boxes  per 

panel 

NPOXES 

Total  number 

of  lifting  surface  boxes 

ITEM  NO. _ 4 

Dimensioning 

data 

cc  1 

2 

3 

4  6  7 

1 _ 

4 

_ 

6 

8  0  2 

N6 

MS  BE 

MBE 

VdXidJble 

NB 

Total  number 

of 

bodies 

MSEE 

Total  number 

of 

slender  body 

elements  for  all  bodies 

MEF 

Total  number 

of 

interference 

body  elements  for  all 

bodies 


ITEM-NOi. 


cc 

1 _ 

ELXIA 

Variable 

ELXIA 
ELY  I A 
2LZIA 


ITEM  NO. 


cc 

1 _ 

PHINA 

Variable 

PHINA 

PHIZA 

PHI  YA 


Coordinates  of  aerodynamic  nodal  points.  Item 
numbers  5  and  6  are  omitted  if  SECT  data  has  been 
input  since  these  data  will  be  obtained  by  the 
program  from  the  SECT  data. 

(Repeated  for  all  aerodynamic  points,  1=1, NODES) 

1  2  3  4  6 

3 _ 5 _ 7 _ 9 _ 1 _ 

ELYIA  ”  ELZ IA~ 

Description 

XAA.S  coordinate  of  aerodynamic  nodal  point 
YAAS  coordinate  of  aerodynamic  nodal  point 
ZAAS  coordinate  of  aerodynamic  nodal  point 


Modeshapes  (Repeated  first  for  all  aerodynamic  points 
1=1, NODES,  and  then  for  all  modes  J= 1 , NS YM+NAS YM) 

1  2  3  4  6 

3 _ 5 _ 2 _ 9 _ 1 _ 

PHIZA  PHI  YA 

Description 

Mode  shape  in  normal  direction  at  panel  aerodynamic 
nodal  point 

Mode  shape  in  z-direction  at  body  aerodynamic  nodal 
point 

Mode  shape  in  y-direction  at  body  aerodynamic  nodal 
point 


6 

_ 2 _ 


7 

2 


IIIM_no. _ 7 

cc  1  2  2 

1 _ z _ i _ 6 

rPFl  IPF2  IPF3 


4 


NGUST 


Variable  Description 

TPK1  Print  flag  for  the  r.ormalwash  and  gust  boundary 

conditions 

=1,  print  all  boundary  conditions 
=0,  no  print 

I PF 2  Print  flag  for  point  forces  that  are  saved  on  the 

AEFO  file 
=1,  print  forces 
=0,  no  print 

I PF  3  Print  flag  for  pressures,  body  loadings,  and  downwash 

factors 

=1,  print  pressures  only 

=2,  print  pressures  and  body  loadings  only 
=3,  print  pressures,  body  loadings,  and  the  down- 
wash  factor  matrix,  DT  and  DTA 
=0,  do  not  print  any  of  the  above 

NGUST  Gust  direction  cosine  override  flag 

=0,  use  standard  direction  cosine  matrix  for  13  gust 
conditions 

=  n,  override  some  (or  all)  elements  of  the  direction 
cosine  matrix  for  the  gust  boundary  conditions. 
n=number  of  total  gust  conditions  (max=20) 

See  input  under  Group  No.  5 


lULtU^Qi. _ 8 

cc  1  2  3  4  6  7 

1 _ l _ 4 _ 6 _ g _ 0 _ l 

NKD  NKP  MK1  MK2 

Variable  Description 

NKD  Number  cf  reduced  frequencies  for  doublet  lattice 

calculations 

NKP  Number  of  reduced  frequencies  for  piston  theory 

calculations 

MK1  sequence  number  of  first  box  on  first  panel  repre¬ 

senting  a  body  surface,  whenever  this  body  is  at  zero 
incidence;  otherwise  MK1=0 

MK2  Sequence  number  of  last  box  on  last  panel  represent¬ 

ing  a  body  surface,  whenever  this  body  is  at  zero 
incidence;  otherwise  MK2=0 

Ngte;  Panels  on  body  surfaces  need  not  be  input  last. 

Maximum  number  of  reduced  frequencies  (NKD+NKP) =100 
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Variable  Description 

FMACH  Mach  number,  usual  definition 

FEFA  Reference  area,  usually  total  area  of  both  wings 

REFS  Reference  semispan 

PEFC  Reference  chord,  usually  average  chord  of  wing 

XM  Moment  axis 

SCALER  Nondime ns ional  image  radius 


ITEM  NO.  10  (Repeated  until  all  FREQ (I)  are  input  for  1=1, NK, 
where  NK  =  NKD+NKP) 

cc  1  2  3  4  6 

1 _ 3 _ 5 _ 7 _ 9 _ 1 _ 

FRFQ  FPEQ  FREQ  ~  FREQ  FREQ  FFEQ 

Variable  Description 

Ffeq  Reduced  frequency 
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GR2UIL.N0, 


ITEM  NO. 


Panel  data  (Items  1  thru  5  repeated  for  all  panels 
N=  1 , NP) 


Variable 


2ei££i£tign 

inboard  leading  edge  x-coordinate 
inboard  trailing  edge  x-coordinate 
outboard  leading  edge  x-coordinate 
outboard  trailing  edge  x-coordinate 


ITEM  N0t _ 2 


Variable  Description 

Y1  inboard  edge  y-coordinate 

V 2  outboard  edge  y-coordinate 

Z 1  inboard  edge  z-coordinate 

Z2  outboard  edge  z-coordinate 


IUE«_NOi. _ 3 


3 

_ 6. 

"iGR  UP 


7 

_ 2 


Variable  Description 

Nc  Number  of  chordwise  boxes  for  panel 

NS  Number  cf  spanwise  strips  for  panel 

IGF  Up  Group  number  of  panel.  Usually,  panels  that  are 

butted  up  trailing  edge  to  leading  edge  are  members 
of  the  same  group,  and  then  the  monitored  span  loads 
will  Le  calculated  for  the  ’combined  strips'  of  such 
panels.  For  panels  in  the  same  IGPUP,  strips  with 
the  3atre  y-coordinate  (for  dihedral  of  4?  degrees  or 
less)  or  the  same  z-coordinate  (for  dih«  i::al  of 
qreater  than  45  degrees)  will  be  added  into  one 
combined  strip  for  the  span  load  calculation  only. 
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ITEM  NO. _ 4  (Fepeated  until  all  TH(I)  are  input  for  I=1,NC+1) 


TH  TH  TH  TH  TH  TH 


Il£i§ble  Description 

TH  Fractional  chordvise  divisions  for  panel.  Usually 

varies  from  0.0  at  the  leading  edge  to  1.0  at  the 
trailing  edge 


ITEM  NO. _ 5  (Repeated  until  all  TAU(I)  are  input  for  I=1#NS+1) 

cc  1  2  3  4  6 

1 _ 2 _ 5 _ 7 _ 9 _ 1 _ 

TAU  TAU  TAU  TAU  TAU  TAU 

V§li§bie  Description 

Fractional  spanvise  divisions  for  panel.  Usually 
varies  from  0.0  at  the  inboard  edge  to  1.0  at  the 
outboard  edge 


TAU 


<2EQ!IE_fc!Q»— 1  Body  data  (Items  1  thru  7  repeated  for  all  bodies, 
N=1,NB.  If  NB=0,  skip  to  Group  No.  4) 


I2EM_N0s. _ 1 


cc 

1 

2 

3  4 

6 

1 

3 

5 

....  7  9 

1 

ZC 

YC  AO 

AF 

Variable 

inscription 

ZC 

z-coordinate 

of  body  axis 

YC 

y- coordinate 

of  body  axis 

AO 

average  characteristic  semi-width  of 

interference  body 

AF 

cross-sectional  aspect  ratio  of  body 

(height/width) 

ITEM  NO, 

_ 2 

cc 

1 

2 

3  4 

6 

7 

1 

2 

4 

6  _  a 

0 

NBE 

NSBE 

NFI  NFS 

NT1 

Variable 

fiescription 

NBE 

number  of  interference  body  elements 

NSBE 

number  c£  slender  body  elements 

NRI 

interference 

•radius*  flag 

=1,  RI  array 

is  input  below 

=0,  FI  ( I)  =  AO,  for  all  I=1,NBEt1 

NFS 

slender  body 

•radius*  flag 

=1,  PS  array 

is  input  below 

=  0,  PS  ( I)  *  AO,  for  all  I=1,NSBE+1 

NT1 

number  of  elements  in  the  TH1  array  below  (max=20) 

iTEM_NO, _ 3  (Pepeated  until  all  XII(I)  are  input  for  I=1,NBE+1) 

cc  1  2  3  4  6 

1 _ 3 _ 5 _ 7 _ 9 _ _ 1 _ 

XII  XII  XII  XII  XII  “  XII 

Variable 


XII 


Description 

x-coordinates  of  interference  body  element  endpoints 


ITEM_N J. _ 4  (Repeated  until  all  RI(I)  are  input  for  I = 1 , N P E+ 1 . 

Omit  if  NRI  =  0) 

cc  1  234  6 

1 _ 3 _ 5 _ 2 _ 9 _ 3 _ 

PI  RI  '  PI  RI  PI  PI 

Variable  Description 

PI  semi-widths  or  radii  of  interference  body  element 

endpoints 

U£M_N3i _ 5  (Repeated  until  all  TH1(I)  are  input  for  I=1,NT1) 

cc  1  2  3  4  6 

1 _ 3 _ 5 _ 7 _ 2 _ 1 _ 

TH1  TH1  TH1  TH1  THl  TH1 

Variable  Description 

THl  anqular  orientation  of  the  points  on  interference 

body  surfaces,  degrees 

Note:  These  points  must  not  lie  near  any  of  the  lifting- 

surface-body  intersections.  Also  for  elliptic  bodies 
THl  is  not  the  angular  orientation  of  a  point  on  the 
body  surface  but  an  angular  orientation  parameter. 

The  actual  location  of  the  point  is: 

Z=AP*A*SIN(TH1) ,  Y=A*COS (THl ) 


ITEM  NO. _ 6  (Repeated  until  all  XIS(Il  are  input  for  I=1,NSBEf1) 


cc 

1  2 

3 

4 

6 

1 

3  5 

7 

9 

1 

XIS 

XIS  XIS 

XIS 

XIS 

XIS 

Hi  iafele 

Description 

XIS 

x- coordinates  of 

slender  body 

element 

endpoints 

ITEM  NO. 

.2  (Repeated  until 

all  R S  ( I )  are 

input  for  I=1,NSBE+1. 

Omit  if  NRS=0) 

cc 

1  2 

3 

4 

6 

1 

3  5 

7 

9 

1 

PS 

RS  PS 

RS 

RS 

RS 

vatiafele 

£e§C£i£tiC2 

PS 

semi-widths  or  radii  of  slender  body 

endpoints 
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GFOUP  NO.  4  Modal  Spline  Interpolation  Data 
ITEM  NO. _ 1_ 

cc  1  2  3  4  6  7 

1 _ 2 _ 4 _ 6 _ s _ Q _ 2 

NSB  NSP  NMAX  KPRINT 

Variable  DgsgliEtign 

NSB  Number  cf  super bodies 

NSP  Number  cf  superpanels  (=0,  if  SECT  input) 

NMAX  Maximum  number  of  aerodynamic  nodal  points  in  any 

one  of  the  superbodies/superpanels 
KPRINT  Print  flag 

*1,  print  H  and  DHDX  matrices 
=0,  do  not  print 


ITEM  NO.  2  Item  numbers  2  thru  4  are  repeated  first  for  all  NSB 
supertodies,  then  for  all  NSP  superpanels. 


cc 

I _ 

1 

_ 2_ 

2 

4 

3 

6 

4 

6 

0 

IF1 

NXQ 

NS  OP 

IXCON 

NACELL 

Variable  Description 

I FI  Flag  to  select  either  surface  spline  or  linear  spline 

calculations  (Not  used  with  SECT  input  data) 

=0,  surface  spline  (used  for  superpanels) 

=1,  linear  spline  (used  for  superbodies) 

NXQ  Number  of  aerodynamic  nodal  points  in  the  current 

superbody /super panel 

NSUP  Number  of  bodies/panels  in  the  current  superbody/ 

superpanel 

IXCON  Flag  to  indicate  x  or  y=constant  for  superbody 

(not  used  for  superpanel) 

*1,  x=constant 

=0,  y=ccnstant  (usual  case) 

NACELL  flag  to  specify  a  body  as  a  nacelle  (Used  only  with 

SECT  input  data) 

*1,  body  is  a  nacelle 
=0,  body  is  not  a  nacelle 


no 


*-  * 


IIJr^_NQ1 _ 3  (Repeated  until  all  INSUP  (I)  are  input  for  I=1,NSUP) 


cc 

1 

1 

2 

2 

4 

3 

6 

4 

8 

6 

0 

7 

2 

INSUP 

INSUP 

INSUP  INSUP 

INSUP 

INSUP 

Variable 

Description 

INSUP 

Array  of  body  numbers/panel 

numbers  in  the 

current 

super body /super panel 

ITEM  NO. 

_ 4 

(Repeated  until 

all  NODE  (I) 

are 

input  for  I 

=  1  ,NXQ) 

cc 

1 

2 

3 

4 

6 

7 

1 

2 

4 

6 

8 

0 

2 

NODE 

NODE 

NODE 

NODE 

NODE 

NODE 

Variable 

££S££i£ticn 

NODE 

Array  of  aerodynamic  nodal 

point 

numbers  in 

the 

current  superbody/superpanel 

ITEM  NO. 

5 

Item  numbers  5 

thru  12  are 

input 

only  if  SECT  data 

has  beer  input. 

cc 

1 

2 

3 

4 

6 

7 

1 

2 

4 _  _ 

6 

8 

0 

2 

NSUPTO  MAXSEC  MAXPF 


Variable 


Description 


NSUPTO 

MAXSEC 

MAXPP 


total  number  of  superpanels 

maximum  number  of  sections  in  any  one  superpanel 
maximum  value  of  IPAIR  (Item  11)  for  any  section 
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V 


ITEM  NO.  6  Item  numbers  6  thru  13  are  Repeated  for  all  NSUPTO 
superpanels.  \ 


2 

4 

'nodes 


3 

_ 6_ 

NEAASS 


NODASS 


Variable 


NODES 


NEAASS 


NODASS 


Description 


Note: 


number  of  sections  in  current  super,  Viel  (see  text  of 
Volume  II  for  definition  of  a  sectio*  \ 
number  of  nodes  in  current  superpanel.  \NODES=0  for 
control  surface  without  nodes  on  its  h^  \geline 
Identification  number  of  the  superpanel  >hich 
contains  the  nodes  which  are  associated  *  ’th  the 
current  control  surface  superpanel.  \ 

The  number  of  nodes  lying  on  the  associates  'super- 
panel  that  are  to  be  used  for  the  current  co.  \tol 
surface  superpanel.  NODASS  must  be  less  than  yr 
equal  tc  the  maximum  number  of  nodes  on  the  \ 
associated  superpanel.  \ 

NEAASS  and  NODASS  pertain  to  the  case  where  a  \ 

superpanel  is  a  control  surface  and  does  not  have  l  ' 
own  nodes.  If  such  is  not  the  case  set  them  tc  zert 


ITEM  NO.  7  Degree  cf  freedom  flags  for  current  superpanel 


cc  1 

1 _ 2. 

H 


Variable  Description 

tl  =H,  use  the  H  degree  of  freedom  for  all  specified 

nodes  for  current  superpanel 
A  =A,  use  the  ALPHA  degree  of  freedom 

T  =T,  use  the  THETA  degree  of  freedom 

P  -B,  use  the  BETA  degree  cf  freedom 

D  =D,  use  the  DELTA  degree  of  freedom 


IIEN  NO.  8  Elastic  axis  endpoint  coordinates. 

cc  1  2  3  a  fc 

1 _ 3 _ 5 _ 7 _ 9 _ 1 _ 

XA  YA  ZA  XP  YE  Zf 

Var iacie  Description 

XA  XAAS  coordinate  of  inboard  end  of  elastic  axis 

YA  YAAS  coordinate  of  inboard  end  of  elastic  axis 

ZA  ZAAS  coordinate  of  inboard  end  of  elastic  axis 

XE  XAAS  coordinate  of  outboard  end  of  elastic  axis 

YL  YAAS  coordinate  of  outboard  end  of  elastic  axis 

21  ZAAS  coordinate  of  outboard  end  of  elastic  axis 

Note:  Endpoint  coordinates  rrust  not  be  the  sarre  as  any 

nodal  point  coordinates.  The  elastic  axis  fcr  a 
control  surface  superpanel  is  the  hingeline  of  the 
control  surface. 


(Pepeated  until  all  NODE  (I)  are  input  fcr  1  =  1, NODES. 
Oirit  if  NCDES=0) 


cc 

1 

1 

2 

2 

4 

3 

6 

4 

8 

6 

0 

7 

2 

NODE 

NO  C  E 

NODE 

NODE 

NODE 

NODE 

Variable 

£escripticn 

NODE 

node  nunbers  in  the  current  superpanel 

IT IX  NO. 

12 

(Eepeated  until  all  NASS(I) 

i  are  input  for  1=1, 

,NODASS. 

Crrit  if  NCDAS5=0) 

cc 

1 

2 

3 

4 

6 

7 

1 

2 

4 

6 

8 

0 

2 

NA 

ss 

NASS 

NASS 

NASS 

NASS 

NASS 

Var lable 

Description 

NAS  S 

Identification  ruirbers  of  nodes  lying  on  an 

associated  superpanel  that 

are  to  be  used  cn  the 

current  control  surface  superpanel 

which  has  no  nodes 

of  i ts  own. 

NASS  iru st  be  c 

i  subset 

of  the  nodes  of 

the  associated  superpanel. 
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IXEtS_K2._li 


Item  numbers  11  thru  13  are  repeated  for  all  NSEC 
sections  for  the  current  superpanel. 


cc  i  2  3467 

1 _ 2 _ 4 _ 6 _ 8 _ Q _ 2 

IFP  IPAIF 


Ce§££ipti£n 


IFP  =1,  section  is  cut  by  parallel  lines 

=0,  section  is  cut  by  radial  lines 
I  PAIR  number  of  pairs  of  box  numbers  to  be  input  in  item  12 


Note:  The  aerodynamic  boxes  lying  in  a  section  must  be 

identified.  These  boxes  can  be  identified  by  the 
first  and  last  box  numbers  of  a  sequence  of  boxes. 
Since  there  may  be  more  than  one  sequence  of 
consecutive  boxes  there  may  be  several  pairs  of 
first-last  box  number  combinations. 


cc  1 

2 

3 

4 

6 

1  J 

5 

7 

9 

1 

XO  YO  ZO 

Vatiable 

Description 

XO 

XAAS  coordinate  of 

reference 

for 

section 

rc 

YAAS  coordinate  of 

reference 

for 

section 

ZO 

ZAAS  coordinate  of 

reference 

for 

section 

ITEN_NOi_13 

(Repeated  for  all 

IPAIR 

pairs) 

cc  1 

2 

3 

4 

6 

7 

1  2 

4 

6 

8 

0 

2 

11 

12 

Variable 

Description 

11 

box  number  of  first  box 

cf  group 

of  boxes  for  the 

current  section 

box  numfcer  of  last  box  of  group  of  boxes  for  the 
current  section 


12 


GFOVP  NO.  5  Gust  Data  (Input  only  if  NGUST*0) 


ITEM  N3t 

cc 

1 _ 

Variable 

NN 


XCCS 

rcos 

zees 


.1  (Repeated  for  each  qust  condition  to  be  added  or 
overridden) 

11  2  3 

2  4 _ 5 _ 7 _ 

NN  XCOS  YCOS~  ZCOS~ 

Description 

A  number  from  1  to  NGUST  indicating  the  gust 
condition  to  be  input.  If  NN  is  less  than  14,  the 
standard  gust  corresponding  to  NN  will  be  overridden 
with  the  following  input,  if  NN  is  qreater  or  equal 
to  14,  the  qust  whose  direction  cosines  are  given 
below  will  be  added  to  the  existing  13  standard  gust 
orientations.  (the  maximum  number  of  standard  plus 
additional  gusts  is  20) 

NN*- 1  indicates  end  of  optional  gust  input. 

cos (alpha),  direction  cosine  from  the  x-axis  to  blast 

center 

cos (beta),  direction  cosine  from  the  y-axis  to  blast 
center 

cos (gamma) ,  direction  cosine  from  the  z-axis  to  blast 
center 
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c.  inertial  Module  Input  Data  (IMOD) 


The  input  data  for  this  module  consists  of  flags  for 
dimensioning  and  module  control,  locations  of  and  deflections  at 
the  inertial  nodal  points,  and  modal  frequencies  and  structural 
damping. 


ITEM  NO,  1  Data  designator  card 

cc  1  2  3  4  6 

1 _ 3 _ 5 _ 2 _ 2 _ 1 _ 

I  MOD 

Variable  Description 

IMOD  =IMOD,  designates  that  the  data  following  are  the 

inertial  module  input  data 


ITEM  NO. 2  Dimensioning  data 

cc  1  2  3  4  6  7 

1 _ 2 _ 4 _ 6 _ 8 _ 0 _ l 

NMS  NENGS  KPFLDS 

Variable  Description 

NMS  number  cf  'right  hand  side'  mass  points,  including 

centerline  masses 

NENGS  number  of  engines  on  right  hand  side  and  centerline 

KPFLDS  =1,  print  unit  load  matrices  generated  (if  any) 

=0,  do  not  print 


ITEM  N'O,  3  Mass  and  location  data.  (Repeated  for  all  masses 
N=  1 , NMS) 

cc  1  2  3  4  6 

I _ 3 _ 5 _ 7 _ 9 _ 1 _ 

EM  ELXIO  ELYIO  ELZIO 

Variable  Description 

EM  mass.  It,  sec.  sec/in 

EI.XIO  XAAS  coordinate  of  mass,  in. 

ELYIO  YAAS  coordinate  of  mass,  in. 

ELZIO  ZAAS  coordinate  of  mass,  in. 
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ITEM  NO.  4  Inertial  modal  frequencies,  FREQ,  and  mode  shapes 
PHI x ,  PHIY,  and  PHIZ.  (Repeated  for  all  modes, 

N= 1 , NOMCD ,  where  NOMOD=NSYM+NASYM.  NSYM  and  NASYM 
are  input  in  the  AERO  data) 

FIRST  CARD 

cc  1  2  3  4  6 

1 _ : _ 5 _ 7 _ 9 _ 1 

F  R  EC? 

Variable  Description 

FREC?  modal  frequency,  Hz 


(Repeated  for  all  mass  points,  1=1, NMS.  Order  mus 
correspond  to  that  of  ITEM  NO.  3) 

cc  1  2  3  4  6 

1 _ 3 _ 5 _ 7 _ 9 _ 1 _ 

PHIX  '  PHIY  PHIZ 


Variable  Description 


PHIX 

mode 

shape 

in 

X 

direction 

at 

mass 

station 

PHIY 

mode 

shape 

in 

y 

direction 

at 

mass 

station 

PHIZ 

mode 

shape 

in 

z 

direction 

at 

mass 

station 

I  TEM_N0j__5 

Modal  definitions. 

AMOCNO. 

FIRST_CAFD 

cc  1 

2 

3 

4 

6 

1  _  _  2 

4 

6 

8 

MODE! 

MODE 2  MODE 3 

M0DE4 

_ 

MODE5 

Pliable  Dgs^ription 

MODE1  mode  number  of  rigid  body  plunge  mode 

MCDE2  mode  number  of  rigid  body  pitch  mode 

MOLF mode  number  of  rigid  body  fore  and  aft  mode 
MODE4  mode  number  of  first  symmetric  elastic  mode 

MODE5  mode  number  of  last  symmetric  elastic  mode 
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SECOND  CngD 


cc 

1 

2 

3  4 

6 

7 

1 

2 

4 

6  8 

.  .  0_  ..  . 

_  l 

M0DE6 

MODE7 

M0DE8  M0DE9 

Variable 

Description 

M0DE6 

mode  number  of 

pitch  trim  mode 

M0DE7 

MODE  8 

mode  number  of 
mode  number  of 

symmetric  -jig  mode 
first  symmetric  mode 

to  be  deleted 

in 

M0DE9 

this  analysis 
mode  number  of 
this  analysis 

last  symmetric  mode 

to  be  deleted 

in 

THUD  CAPP 

cc  1 

2 

3  4 

6 

1  2 

4 

6  8 

_  0  _  . .  .. 

M0DE1 1 

MODE 12 

MODE 13  MODE 14 

MODE 15 

Variably 

Bg§££iEiion 

MODE11 

mode 

number 

of 

rigid  body  roll  mode 

MODE12 

mode 

number 

of 

rigid  bcdy  yaw  mode 

MODS13 

mode 

number 

of 

rigid  body  lateral 

mode 

M0DE14 

mode 

number 

of 

first  antisymmetric 

elastic  mode 

MOD  21 5 

mode 

number 

of 

last  antisymmetric 

elastic  mode 

FOUFTH  CAFD 

cc  1 

2 

3  « 

6 

1 _ 2_ 

4 

6  8 

0 

MODE 16 

MODE  17 

MODE18  MO DEI 9 

MODE 20 

liable 

Bgscj 

riPtion 

MODE16 

mode 

number 

of 

roll  trim  mode 

M0DE17 

mode 

number 

of 

yaw  trim  mode 

MODE 18 

mode 

number 

of 

antisymmetric  jig  mode 

MODE1 9 

mode 

number 

of 

first  antisymmetric 

mode  to  be  deleted 

in  this  analysis 

MOL  220 

mode 

number 

of 

last  antisymmetric 

mode  to  be  deleted 

in  this  analysis 

Note: 

Trim 

modes  may 

also  be  defined  as 

elastic  modes.  For 

example,  the  yaw  trim  mode  (usually  rudder  rotation) 
might  also  be  used  in  the  yaw  damper  system  and  hence 
be  defined  as  an  elastic  mode,  a  trim  mode  and  in  the 
ACS  definition. 
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I1£M_N0s._6 


Modal  structural  dampinq,  CPXDPG. 

(Fepeated  until  the  damping  for  all  modes  has  been 
read  in  fcr  N=1,NOMOD) 


cc 

1 

2 

3 

4 

6 

1 _ 

_ 3 _ 

_ 5 _ 

_ 7 _ 

9 

1 

CPXDPG 

CPXDPG 

“CPXDPG 

CPXDPG 

CPXDPG 

CPXDPG 

liable 

CPXDPG 

Note: 


Description 

modal  structural  dampinq 

The  structural  dampinq  should  be  input  zero 
in  all  modes  except  the  elastic  modes. 
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d.  Unit  Load  Module  Input  Data  (LOAD) 


The  input  data  for  this  module  consists  of  flags  for 
dimensioning  and  module  control  and  all  the  data  necessary  to 
define  the  desired  loads. 


ITEM  NO. _ 1  Data  designator  card 

cc  1  2  3  4  6 

1 _ 3 _ 5 _ _ 7 _ _9 _ 1 

LOAD 


Variable  Description 

LOAD  =  LOAD,  designates  that  the  data  following  are  the 

unit  lead  module  input  data 


ITEM. 

NQ.  2 

Dimensioning 

data 

cc 

1 

2 

3 

4 

6 

7 

1 

2 

4 

_ 

8 

0 

JL 

NBEAMS 

NINTLD 

NSTPSS 

NMGRP 

NABGFP 

NSBGRP 

Variable  Description 


NBEAMS 

NINTLD 

NSTFSS 

NMGPP 

NABGPP 

NSPGFP 


number  of  beams  for  integrated  loads 
number  cf  integrated  loads 
number  of  stresses 

number  cf  groups  of  masses  associated  with  beams 
number  of  groups  of  aerodynamic  boxes  associated 
with  beams 

number  of  groups  of  aerodynamic  slender  body 
elements  associated  with  beams 


ITEM  N 3 

Beam  geometry,  EEAMGM 

(Pepeated  for  all 

beams. 

N*1. 

NBEAMS) 

fIJ:§i_c&PC 

cc 

i 

2 

3 

4 

1 

3 

5 

7 

9 

XI 

y  i 

'll 

Variable 

Description 

XI 

XAA.S  coordinate  of 

inner 

end 

of  beam. 

in, 

ri 

YAAS  coordinate  of 

inner 

end 

of  beam. 

in 

m 

ZAAS  coordinate  of 

inner 

end 

of  beam. 

in, 
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SgCOKD_C££D 


2 


3 

7 _ 

CODE 


a 


6 


t 


cc  1 

1 _ 3_ 

XO  ~  YO 


5 _ 

ZO 


9  1 


Variable 


inscription 


xo 

YO 

ZO 

COLE 


XAAS  coordinate  of  outer  end  of  beam,  in. 

YAA  S  coordinate  of  outer  end  of  beam,  in. 

ZAAS  coordinate  of  outer  end  of  beam,  in. 

Code  defining  the  component  to  which  the  beam  belongs 

=1.0,  wing  or  horizontal  tail 

=2.0,  fuselage 

=3.0,  vertical  stabilizer 

=4.0,  winq  pod 

=5.0,  fuselage  pod 

=6.0,  vertical  stabilizer  pod 


ITEM  NO.  4  Fequired  integrated  loads,  STALDS 

(Fepeated  for  all  loads,  N= 1 , NINTLD) 

FlrST  CAPD 


cc 

1 _ 

PFA.MNO 


1 


CODEL 


2  3  4 

_ 5 _ 7 _ i 

CODEC 


6 

1 


Variable  Description 


BSA.NUO 


COL  EL 


CODEC 


beam  number  of  the  beam  with  which  the  integrated 
load  is  associated, 
load  code 

=1.0,  integrated  load  is  moment  Mx 
=2.0,  integrated  load  is  moment  My 
=3.0,  integrated  load  is  moment  Mz 
=4.0,  integrated  load  is  shear  Pz 
=5.0,  integrated  load  is  shear  Py 
=6.0,  integrated  load  is  shear  Px 

code  defining  the  component  to  which  the  load  belongs 

=1.0,  wing  or  horizontal  tail 

=2.0,  fuselage 

=3.0,  vertical  stabilizer 

=4.0,  wing  pod 

=5.0,  fuselage  pod 

=6.0,  vertical  stabilizer  pod 


S££CND_Ca£D 


2 

5 _ 

ZL 


6 


cc 


IdllaMe 

XL 

YL 

ZL 

PMF.X 

PMIN 


ITEM  NO.  £ 


CC 

1 _ 

STF  ESS 

Variable 

STRESS 


II£M_NOa.j6 


cc 

1 _ 

NFM 

Variable 


1 

3_ 

TL 


3 

7 _ 

PMAX 


4 

9 _ 

PMIN 


1 


Description 

XAAS  coordinate  of  location  of  integrated  load,  in. 

YAAS  coordinate  of  location  of  integrated  load,  in. 

ZAAS  coordinate  of  location  of  integrated  load,  in. 

maximum  allowable  positive  load 

maximum  allowable  negative  load  {a  negative  number) 


Matrix  defining  local  stresses  in  terms  of  the 
integrated  loads,  STRESS. 

(Repeated  until  all  of  the  STRESS  matrix  has 
been  read  in  by  rows.  Omit  if  NSTRSS=0) 

1  2  3  4  6 

3 _ 5 _ 7 _ 9 _ 1 _ 

STRESS  STRESS  STRESS  STRESS  STRESS 

Description 

element  of  STRESS  matrix  relating  local  stresses  to 
the  integrated  loads 


Definition  of  groups  of  masses  associated  with 
beams,  NFNLMB.  (Repeated  for  N=1,NMGFP) 

1  2  3  4  6 

3 _ 5 _ 7 _ 9 _ 1 

NLM  EEAMNO 

Description 

number  of  first  mass  in  group 

number  of  last  mass  in  group 

beam  number  to  which  group  is  associated 


NFM 

NLM 

BEAMNO 


IT£M  NO,  7  Direction  cosine  matrix  relating  mass  point  deflec¬ 
tion  to  the  AAS  coordinate  deflections  (x,y,z, 
ordered).  (Repeated  until  three  rows  for  each  mass 
group  have  been  read  in,  1=1,3,  N=1,NMGFP) 


cc  1  2  3  4  6 

1 _ 3 _ 5 _ 7 _ 9 _ 1. 


TLAMV1  ~~~ 

TLAMV2  TLAMV  3 

Variable 

Description 

TLAMV 1 

first  element  of 

the  ith  row  of  TLAMV 

matrix 

TLAMV 2 

second  element  of  the  ith  row  of  TLAMV  matrix 

TLAMV 3 

third  element  of 

the  ith  row  of  TLAMV 

matrix 

ITEM  NJ.  8 

Definition  of  groups  of  aero  boxes  associated  with 

beams,  NFNLAB. 

(Repeated  for  N=1,NABGRP) 

cc 

1 

2 

3  4 

6 

1 

3 

5 

7  9 

1 

NF/-3 

NLAB  EEAMNO 

V^Ll^ble 

Description 

NFA3 

number  of  first 

aero  box  in  group 

NLAB 

number  of  last  aero  box  in  group 

BEAMNO 

beam  number  with 

which  group  is  associated 

ITEM  NO.  9 

Definition  of  groups  of  slender  body 

elements 

associated  with 

beams,  NFNLSB.  (Repeated  for 

N=  1 ,  NSEGF  P) 

cc 

1 

2 

3  4 

6 

1 

5 

_Z_  2_ 

1 

NFSP  NLSB  EEAMNO 


Ijr Besgrjptjon 


NFSB 

NLSB 

BEAMNO 


number  of  first  slender  body  element  in  group 
number  of  last  slender  body  element  in  group 
beam  number  with  which  group  is  associated 


4 
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ITEM  N3. _ 12  Definition  of  engine  mass  points,  NENGM. 

(Fepeated  for  N=1,NENGS.  Omit  if  NENGS=0 ) 

cc  1  2  3  U  6 

1 _ 3 _  _ 5  2 _ 9 _ _1 _ 

ENG Ml  ENGM2 

Variable  BS§gription 

ENGM1  number  of  first  mass  defining  the  engine 

ENGM2  number  of  second  mass  defining  the  engine 

Note:  Force  due  to  thrust  acts  at  ENGMl  along  a  vector 

pointed  from  ENGM2  to  ENGMl 
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fUN  DATA  DECK  INPUT  DESCRIPTION 


4  , 


The  Fun  Data  Deck  must  always  be  available  and  is  the 
standard  fortran  input  unit  5  (TAPE5)  . 


2£OUP_NOJ_1^  Identification  of  configuration,  print  control, 
control  of  execution,  and  case  constants 


ITEM  N3. _ X  Identification  of  configuration  and  print  control 

cc  1  23467 

1 _ 2 _ 4 _ _  _6 _ 8 _ _ y _ _  2. 

I  DENT 

Variable 

I  DENT 
KPPCXQ 
KPFCHO 

KPf CXL 

KPf CHL 


KPRCXO  KPRCHO  KPPCXL  KPRCHL 

Description 

identification  number  of  configuration 
=n,  print  generalized  response  for  orientation  n 
=  n,  prir.t  check  matrices  in  frequency  response 
module  for  orientation  n 

=  n ,  print  integrated  load  modulii  due  to  unit  gust 
for  cnentation  n 

=n,  print  matrix  data  for  unit  qust  calculations  for 
orientation  n 
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ITEM  NO.  2  Module  execution  control  card 


cc 

1 _ 4. 

5  PCD 


1 

.9 _ 2. 

OPCD 


1  2 

J_ _ 0. 

OPCD 


2  2 

_5 _ 8_ 

OPCD 


3  3 

.3 _ 6_ 

OPCD 


4  4 

1 _ 4. 

OPCD 


Variable 


OPCD 


Note: 


Description 

defines  nodules  to  be  executed 
=AEPO,  execute  aerodynamic  module 
=  (JNI T,  execute  unit  loads  module 
=ACSM,  execute  active  control  system  module 
=FRSP,  execute  frequency  response  module 
=GUST,  execute  unit  gust  loads  module 
=BLST,  execute  trim  loads,  blast,  and  time  history 
nodules 

=RIGE,  execute  rigid  loads  module 
=MEFG,  execute  aero  file  merge  module 

The  user  must  specify  which  modules  are  to  be 
executed  in  any  given  pass  through  the  program.  If  a 
module  is  specified  to  be  executed  and  sufficient 
data  is  net  available  (from  any  required  source)  the 
program  will  halt  with  an  explanatory  statement.  The 
active  control  system  module  need  only  be  executed 
when  running  with  FRSP  and/or  GOST  specified,  but  the 
control  system  definition  must  be  the  same  if  these 
two  modules  are  not  run  in  the  same  pass. 


ITEM  NO,  3  Case  constants  (If  MEPG  or  AEPO  has  been  specified  on 
the  module  control  card,  omit  this  item) 


VKEAS 


2 

5 _ 

SIZFCT 


Vaii^ble  Description 

ALT  altitude,  ft 

VKEAS  velocity,  equivalent  airspeed  in  knots 

SIZFCT  size  factor  defining  the  units  of  the  input  geometry 

data  relative  to  inches.  default=1.0 
PLQ  =n,  plot  generalized  response  for  nth  orientation 

=0.0,  do  not  plot 

PLL  =n,  plot  nth  load  due  to  unit  gust  for  all 

orientations 
=0.0,  dc  not  plot 
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GftOUP  NO.  2  Active  control  system.  (If  ACSM  has  not  been 

specified  on  the  module  control  card,  omit  this 
qroup) 

Input  modes  must  contain  a  mode  which  describes  the  motion  of 
the  surface  drivinq  the  system  as  commanded  by  the  active  system. 
The  input  modal  amplitude  of  this  mcde  must  be  consistent  with  the 
commanded/ sensed  movement  as  defined  by  the  transfer  function 
input  data.  The  scalar  multiplier  may  be  used  for  this  purpose. 
Note  that  several  transfer  functions  may  be  used  to  drive  any 
sir.qle  mode,  but  any  single  transfer  function  may  drive  only  one 
mode.  Thus  if  twc  surfaces  are  being  used  and  driven  by  the  same 
mathematical  transfer  function  and  two  separate  modes  are  used  to 
define  these  control  surface  motions,  the  transfer  function  must 
be  entered  twice  and  defined  by  two  different  transfer  function 
numbers. 


ITEM  NO.  1 

Control  and 

dimensioning 

data 

cc  1 

2 

3 

4 

6 

7 

1  2 

U 

6 

8 

0 

_  l 

NTFS 

NTFA 

MXBLK 

MXOBLK 

Variable 


Description 


NTFS 

NTFA 

MXELK 

MX05LK 


number  cf  symmetric  transfer  functions 
number  of  antisymmetric  transfer  functions 
maximum  number  of  blocks  input  in  any  single  transfer 
function  with  items  3  or  5 

maximum  order  of  any  block  polynomial  in  input  trans¬ 
fer  function  blocks  with  items  3  or  5 
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ITEM  NO.  2  Data  describing  kinematics  of  active  control  system 

for  symmetric  control,  TMSS.  (Repeated  for  1=1# NTFS. 
If  NTFS=C ,  skip  to  item  no.  4) 


cc 

1 _ 

TM6  31 


1 

3 _ 

TMSS  2 


5 _ 

TMSS  3 


3 

7 _ 

TMSS  4 


U 

9 


6 

1 


v&Ei^ble 


Descripticn 


IMF  SI 
TMSS2 


IMS  S3 
TMES4 


mass  point  number  of  sensed  mass  point 
degree  of  freedom  number  of  sensed  degree  of  freedom 
* 1 . c  # ph I X ;  =2.0, PHI  X ;  =3.0, PHIZ;  or  if  SECT  input: 
=1.0, F;  =2. 0 , L;  =3.0, H;  =4.0, THETA;  =5.0, ALPHA; 

=6 . 0 , PSI ;  =7.0, BETA;  =8. 0, DELTA 

scalar  multiplier  on  sensed  degree  of  freedom  motion 
mode  number  of  the  driving  mode  (should  be  a  control 
surface  mode) . 


XjiM  Nd.._3  Data  defining  elements  of  symmetric  transfer 
f  unction. 


cc 

NTF 


1  2  3  4  6 

3 _ 5 _ 2 _ 2 _ 3 _ 

NBLK  “  NOFD  AN  AD  TLAST 


Variable  Description 


njxF  transfer  function  number 

=-1.0,  no  more  input  data  for  this  item 
NPX ff  block  number  of  these  elements 

;;OF-D  power  of  LaPlace  operator  (s)  ,  associated  with  these 

element  s 

an  numerator  coefficient 

,*.D  denominator  coefficient 

TLAST  last  transfer  function  associated  with  this  data 

=0.0,  entry  is  for  transfer  function  NTF  only 
*0.0,  entry  is  for  transfer  function  NTF  thru  TLAST 
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■s 


.n 


ITEM  NO.  4  Data  describing  kinematics  of  active  control  system 
for  antisymmetric  control,  TMSA. 

(Fepeated  for  I=1,NTFA.  If  NTFA=0,  skip  to  group  3) 


cc 

1 

2 

3 

4 

6 

1 

3 

5 

7 

9 

1 

T  MSA  1 

TMSA2 

TMSA3 

TMSA4 

Variable  Description 

TMSA 1  masr  point  number  of  sensed  mass  point 

TMSA 2  degree  of  freedom  number  of  sensed  degree  of  freedom 

=1.0 ,  PHI  X ;  =2.0,PHIY;  =3.0, PHIZ;  or  if  SECT  input: 

=  1.0, f;  =2.  0,L;  =3.0, H;  =4.0, THETA;  =5.0, ALPHA; 

=6.0, PSI ;  =7.0, BETA;  =8.0, DELTA 

TMSA 3  scalar  multiplier  on  sensed  degree  of  freedom  motion 

TMSAU  mode  number  of  the  driving  mode 


_ 5  Data  defining  elements  of  antisymmetric  transfer 

f unctior . 

cc  1  2  3  u  6 

I _ 3 _ 5 _ 2 _ 9 _ 1 _ 

NTF  NBLK  NOFD  AN  AD  TLAST 

Variable  Description 

NTF  transfer  function  number 

=-1.0,  no  more  input  data  for  this  item 

NBLK  block  number  of  these  elements 

NOFD  power  of  LaPlace  operator  (s) ,  associated  with  these 

element  s 

AN  numerator  coefficient 

AD  denominator  coefficient 

TLAST  last  transfer  function  associated  with  these  data 

=0.0,  entry  is  for  transfer  function  NTF  only 
*0.0,  entry  is  for  transfer  function  NTF  thru  TLAST 

Nfitex  Inputing  AN,  AD  and  ALAST,  blank  fields  are  not  the  same 
as  0.0;  therefore  0,0's  must  be  input  explicitly. 
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GPOUP  NO.  3  Frequency  response,  (if  FRSP  has  not  been  specified 
on  the  irodule  control  card,  omit  this  group) 

ITEM  NO.  1  Dimensioning  and  control  data 

cc  1  2  3  4  6  7 

1 _ 2 _ 4 _ 6 _ £ _ 0 _ 2 

NFRGR 

Variable  De scription 

NFRGF  number  of  frequency  groups  to  be  input  in  item  no.  2 


ITEM  NO.  2  Frequency  groups  (Repeated  for  N=1,NFFGF) 

cc  1  2  3  4  6 

1 _ 3 _  5 _ 7 _ _ 9 _ 1 

FI  F2  DF 

Variable  Description 

FI  starting  frequency  (hz)  for  this  group 

F2  ending  frequency  (hz)  for  this  group 

DF  incremental  frequency  for  this  group 


Note:  input  frequencies  roust  be  entered  in  an  ascending 

order  and  groups  must  not  overlap,  nor  should  any 
frequency  be  repeated. 
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GROUP  NO.  u  Unit  gust  loads.  (If  GUST  has  not  been  specified  on 
the  module  control  card,  omit  this  group) 


ITFM  NO.  1  Dimensioning  and  control  data 

cc  1  23467 

i _ i _ a _ 6 _ a _ a _ 2 

NACC 

Variable  De  script],  on 

NACC  number  of  mass  points  for  acceleration  calculations 


ITEM  NO.  ?  Acceleration  indexes  (Repeated  until  all  mass  point 
numbers  have  been  input  for  N=1,NACC.  Omit  if 
NACC=0) 

cc  1  2  3  4  6  7 

1 _ 2 _ 4 _ 6 _ a _ Q _ 2 

MASSNO  I DOF  MASSNO  IDOF  MASSNO  IDOF 

Variable  Description 

mass  point  number  for  which  acceleration  is  desired 
degree  cf  freedom  for  which  acceleration  is  desired 
=1,  x-direction 
=2,  y-direction 
=3,  z-direction 


MASSNO 

IDOF 


GSOUP_Ng^ 


ITEM  NO. 
cc 

1 _ 

AN 

Variable 

AN 

KMA.N 


ZDOT 

KPFTFM 

FPF ADF 


IIESLL'Qi.. 

cc 

1 _ 

TIMFMX 

Variable 

TI MEMX 

DELT 

EFF 

hoft 

KFFT.VI1 

KPLTMH 


Ngt^: 


£  Trim  and  blast  data.  (If  BLST  has  not  been  specified 
on  the  module  control  card,  omit  this  group) 

Trim  corstants 


1 

2 

3 

4 

6 

3 

5 

7 

9 

1 

KMAN 

ZDOT 

KPRTRM 

PBPADF 

Description 

load  factor 
manuever  constant 
=0.0,  no  manuever 

=1.0,  symmetric  pull-up  or  push-over 
=  2.0,  turn 

climb  rate,  fps  (input  only  if  KMAN=2.0) 

=1.0,  print  all  matrices  in  trim  solution 
=0.0,  dc  net  print 

scale  factor  for  rotation  modes  input  at  ether  than 
for  ere  radian  rotation,  =  1  rad/value  used  in 
radians  (def ault= 1 . 0) 


Time  history  constant  and  blast  characteristics 


1 

2 

3 

4 

6 

3 _ 

_ 5_ 

_ 7 

_ 9 _ 

J _ 

DELT 

EFP 

HGPD 

KPPTMH 

KPLTMH 

Description 

maximum  time  for  time  history,  sec 
delta  time  for  time  histories 
weapon  yield,  KT 

heiqht  of  ground  above  sea  level,  ft 
=1.0,  print  total  load  time  histories 
=0.0,  do  not  print 

=1.0,  plot  load  time  histories  for  FHS  of  aircraft 
=2.0,  plot  load  time  histories  for  LHS  of  aircraft 
=3.0,  plot  both  LHS  and  RHS 
=0.0,  dc  not  plot 

DELT  is  increased  by  a  factor  of  5  after  0.25  seconds 
of  response  and  again  by  a  factor  of  2  after  1.0 
seconds . 


132 


I2EM_N0.5__3  Control  flags 


cc 

1 

2 

3 

4 

6 

7 

1 

2 

4 

6 

_ 8 _ 

£  _ 

2 

NORMAX 

KGRD 

KLPT 

KLOAD 

NCFITS 

KPRBLS 

Variable  Description 

NORMAX  maximum  number  of  orientations  to  be  considered 

KGRE  control  constant  for  ground  reflection 

=0,  no  ground  reflection 
=1,  include  ground  reflection 
KLPT  =  1,  iteration  for  critical  range  is  desired 

=0,  nc  iteration 

KLCAI  =  1,  new  maximum  allowable  loads  input  in  item  no.  7 

NCP ITS  =1,  maximum  stresses  input  in  item  no.  8 

KPRBLS  =1,  print  blast  matrices  for  checkout 

=  n,  do  rot  print 


ITEM  MO. _ 4  Thrust  input  data,  THRUST.  (Repeat  until  all  thrust 

(I)  have  been  input  for  1-1 ,NENGS.  If  NERGS=0 , 
omit  this  item) 


cc 

1 

2 

3 

4 

6 

1 

3 

_ 5 _ 

7 

_ 9  _ 

1 

THF  UST 

THRUST 

THRUST 

THRUST 

THRUST 

THRUST 

Variable  Description 

THRUST  engine  thrust,  lb. 


ITEM_NO, _ 5  orientation  numbers,  NOFS.  (Repeat  until  all  NORS 

(I)  have  been  input  for  1=1,  NORMAX) 

cc  1  2  3  4  6  7 

1 _ 2 _ « _ 6 _ S _ 0 _ l 

NCR  NCR  NOR  ~n5p  NOR  NOR 

V i^ble  Description 

NOR  orientation  numoer  of  blast  for  desired  solution 

Note:  NORMAX  and  NOF  must  be  consistent  with  the  set  of 

orientations  established  during  execution  of  the  FRSP 
module  and  available  on  the  AERO  file. 
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ITEM  NO.  6 


Ranqe  estimates,  PEST.  (repeat  until  all  PEST(I) 
have  been  input  for  l  =  1,NOPMAX) 


cc 

1 

1 

3 

2 

5 

3 

7 

4 

9 

6 

1 

PEST 

REST 

FEST 

REST 

REST 

REST 

Variable 

Description 

PEST 

estimated  ranqe  at 

shock  arrival  time,  ft 

ITEM  NO,. 

_7 

Maximum 

allowable 

loads,  STALDS.  (Repeat 

for 

N= 1 , NINTLD. 

If  KLOAD=0 ,  omit  this  item) 

cc 

1 

2 

3 

4 

6 

1 

3 

5 

7 

9 

1 

PM7.X 

PMIN 

Variable 

Description 

PMAX 

maximum 

allowable 

positive  load 

PMIN 

maximum 

allowable 

negative  load  (a  negative  number) 

ITEM  NO. 

_8 

Maximum 

allowable 

stresses. 

(Repeat  for 

N=1  , 

NINTLD. 

If  NCF ITS 

i=0,  omit  this  item) 

CC 

1 

2 

3 

4 

6 

1 

_3_ 

5 

7 _ 

9 

1 

3MAX  SMIN 


Variable  Description 


SMAX  maximum  allowable  positive  stress 

SMIN  maximum  allowable  neqative  stress  (a  negative  number) 


GROUP  NO.  6  Figid  loads.  (If  FIGD  has  not  been  specified  on  the 
module  ccntrol  card,  omit  this  group) 


ITEM  NO.  1  Dimensioning  and  control  data 


cc 

1 

2 

3 

4 

6 

7 

1__ 

_ 2 _ 

4 

6 

8 

_ a _ 

.  2 

j  ; 

NTMGST 

NBCX 

NSBE 

NOR 

I  PLOT 

IPREQ 

, 

Variable  Description 

NTMGST  number  of  input  time  points  for  gust 

NBOX  number  of  aerodynamic  boxes  for  which  time  histories 

of  pressures  are  desired 

NSPE  number  cf  slender  body  elements  for  which  time 

histories  of  pressures  are  desired 

NOE  gust  orientation  requested 

I PIOT  plot  flag  for  Fourier  transform  of  RHO*VG  and 

pressures 

IFFEO  print  flag  for  checkout  print 


ITEM  NO,  2  Time  history  constants 

cc  1  2  3  4  6 

tTmemx  delt  *  .  ~ 

Variable  pescjipticn 

TIMEMX  maximum  time  for  output  time  history,  sec 

DELT  delta  time  for  output  time  histories,  sec 

Note:  DELT  is  increased  by  a  factor  of  5  after  0,25  seconds 

of  response  and  again  by  a  factor  of  2  after  1.0 
seconds . 


ITEM  NO.  3  Number  cf  frequency  groups 


NFPGP 


Yat ia£l£ 

NFFGP  number  cf  frequency  groups  to  be  input  in  item  no.  4 
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ITEM  NO,  4 

Frequency  groups  (Fepeated  for  N=1,NFFGF) 

cc 

1 

2  3 

4 

6 

j 

1 

3 

5  7 

9 

1 

fi 

F2 

CF 

4 

Variable 

Description 

FI 

starting  frequency  (hz)  for  this  group 

F2 

ending  frequency  (hz) 

for  this  group 

DF 

incremental  frequency 

for  this  group 

Note: 

input  frequencies  must 

be  entered  in  an  ascending 

order  and  groups  trust  not  overlap,  nor  should  any 
frequency  be  repeated. 


ITLM  NO.  5  Aerodynamic  box  numbers  at  which  presssures  are 

desired.  (Pepeat  until  all  IBOXES(I)  have  been  input 
for  N= 1 ,  NEOX) 

cc  1  2  3  4  6  7 

i _ 2 _ 4 _ 6 _ g _ 0 _ 2 

I BOX  IBCX  I BOX  IBCX  IBOX  IBOX 

Variable  Description 

I PC X  aerodynamic  box  number  at  which  pressure  is  desired 


ITEM  NO,  6  Slender  body  element  numbers  at  which  pressures  are 

desired.  (Pepeat  until  all  ISBE(I)  have  been  read  in 
for  K= 1 ,NSEE) 

cc  1  2  3  4  6 

1 _ 2 _ 4 _ _ 6 _ 8 _ £ _ 

"iSBE  “  IDIF  ISBE  IDIF  ISBE  ~  IDI 

Variable  Description 

ISBE  slender  body  element  number  at  which  pressure  is 

desi red 

IDIF  force  direction:  =1,  z-direction;  =2,  y-direction 
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IIEM_NJJ._2  Input  generalized  coordinates.  (Fepeat  until  all 
0(1)  have  been  input  for  N=1,NOMOD) 


cc 

1 

2  3 

4 

6 

1 

.3. 

5  7 

.9 _ 

J _  _ 

0 

0 

0  0 

Q 

Q 

Variable 

Description 

0 

input  generalized  coordinate 

lTeM_NO. _8 

Input  time  histories  of  overpressure. 

density,  and 

material  velocity.  (Repeat  until  all 

time  points 

have  teen  input  for  N=1 , NTMGST) 

cc 

1 

2  3 

4 

6 

1 

3 

5  7 

_9 _ 

1 

TIKE 

DP  FHO  VG 

V  ariable 

Description 

tike 

time,  seconds 

dp 

overpressure,  psi 

PHC 

density,  lb- sec**2/f t**4 

VG 

material  velocity,  ft/sec 
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GROUP  NO.  7  Merqe  o t  aerodynamic  files  from  TAPE17  and  TAPE18 

onto  TAFE19.  (If  MERG  has  not  been  specified  on  the 
module  control  card,  omit  this  qroup) 


ITEM  NO.  1  Number  cf  input  reduced  frequencies 


cc 

1 


1 

_ 2. 

INK  1 


Variable 

INK1 
INK  2 


2 

_ 4_ 

INK2 


3 

6 


4 

8 


6 

0 


Description 

number  of  frequencies  from  TAPE17  file  to  be  merged 
number  cf  frequencies  from  TAPE18  file  to  be  merged 


iISM_N  _ 2 

cc 

1 _ 

RKIN1 

Variable 

PKIN1 


Reduced  frequencies  from  TAPE17  file,  (Repeated 
until  all  RKINI(I)  have  been  input  for  I*1,INK1) 


1 

3 _ 

RKIN1 


2 

PKINl" 


3 

7 _ 

'pkinT 


4 

9 


6 

1 


RKIN1 


RKIN1 


pescriptjcn 


value  of  reduced  frequency  from  TAPE17  file  for 
which  aero  is  to  be  merged 


ITEM_N2i,  3  Reduced  frequencies  from  TAPE18  file.  (Repeated 
until  all  RKIN2 (I )  have  been  input  for  1=1, INK 2) 


cc 

1 _ 

SKIN2 

V 51 i able 

PKIN2 


1 

.3 _ 

PKIN2 


2 

5 


4 

9 


PKIN2 


RKIN2 


RKIN2 


6 

J _ 

FKIN2 


Description 


value  cf  reduced  frequency  from  TAPE18  file  for 
which  aero  is  to  be  merged 
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(to  -4 


5 


DISCUSSION  OF  INPUT  DATR 


Fixed  data  deck  input  for  the  structural  model  is  essentially 
the  same  as  that  given  in  Reference  1  with  the  exception  of  the 
aerodynamic  model  data  and  the  unit  lead  data.  However,  rigid 
body  modes  and  trim  must  be  input  to  both  the  aerodynamic  module 
and  the  inertial  module.  Rigid  body  modes  should  consist  of  rigid 
body  plunge,  pitch,  fore  and  aft,  roll,  yaw,  and  lateral 
translation.  Trim  modes  describe  surface  motion  necessary  to  trim 
the  aircraft  and  consist  of  a  mode  each  for  pitch,  roll,  and  yaw. 
These  and  the  rigid  body  rotation  modes  should  consist  of  the 
control  point  deflections  for  one  radian  of  rotation.  In  the 
event  that  some  other  reference  rotation  is  used,  then  a 
correction  factor  must  be  input  to  the  trim  module  (RBFADF). 
Details  of  aerodynamic  data  preparation  are  given  in  Volume  II  of 
this  repert.  The  choice  of  reduced  frequencies  for  the 
aerodynamics  module  must  include  k=0  and  the  second  reduced 
frequency  should  be  small  (approximately  0.1).  The  maximum  value 
should  be  high  enouqh  to  extend  the  frequency  somewhat  beyond  the 
hiqhest  natural  free-free  frequency  of  the  elastic  model.  Reduced 
frequency  spacing  should  be  relatively  uniform  and  consist  of 
about  4  values  for  each  1.0  k  used,  though  3  will  often  suffice. 

Unit  load  data  required  is  explained  in  the  unit  load  module 
description.  Care  should  be  taken  with  this  data  to  insure  that 
the  beam  segments  used  to  define  the  load  paths  are  contiguous  or 
that  all  appropriate  load  stations  (mass  and  aero)  feed  load  to 
the  appropriate  beam  segment. 
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Active  systeir  definition  input  data  is  discussed  in  that 
mod  ule. 

The  choice  of  frequencies  for  solution  depends  on  the  vehicle 
under  analysis.  Small  increments  of  frequency  should  be  used  in 
the  aircraft  riqid  body  frequency  ranqe  and  near  the  elastic  modal 
frequencies.  Broader  spacinq  may  be  used  between  separated  modes. 
In  the  event  a  frequency  is  input  which  is  beyond  the  maximum 
definable  by  the  unsteady  aerodynamic  matrices  available,  the 
proqram  will  desrqr.ate  it  as  unusable.  The  maximum  frequency 
available  is  qiven  by 

f  =  Vk/6. 28b 

where  V  is  the  true  airspeed,  k  is  the  maximum  reduced  frequency 
available,  and  b  is  the  reference  semi-chord. 

fun  data  are  as  described  in  that  section,  except  that  CELT 
should  be  small  encuqh  to  qive  the  response  of  the  hiqhest 
freuuer.cy  mode,  DELT  qreater  than  or  equal  0.2f.  The  time  of 
solution  (TIMEMX)  should  be  at  least  5S/Vss  where  S  is  the  maximum 
lenqth  of  the  vehicle  (tip  tc  tip  or  nose  to  tail)  and  Vss  is  the 
speed  of  sound  at  fliqht  altitude. 
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SECTION  V 


PROGRAM  OUTPUT 

The  program  output  is  generally  self-explanatory  in  that  all  output 
data  are  preceded  by  headings  and  variable  definitions  which  are  for  the  most 
part  consistent  with  the  data  definitions  found  in  the  module  description 
and  input  data  sections  of  the  report. 

Input  data  from  the  Fixed  and  Run  Data  Decks  for  any  given  run  are 
printed  out  for  the  purpose  of  checking  input  data. 

Generated  data  riles,  output  by  the  several  modules,  are  described  in 
Section  VI. 

The  use  of  a  generated  data  file  is  followed  by  a  printout  of  the 
data  file  header  which  identifies  the  file  and  gives  a  brief  summary  of  the 
data  on  the  file.  A  rewind  of  any  file  initiates  a  message  to  that  effect. 

During  the  execution  of  any  module,  messages  informing  the  user  of  the 
number  of  incremental  words  (in  decimal)  required  for  the  module  are  printed 

Several  auxiliary  print  flags  are  available  for  use  throughout  the 
program.  These  are  intended  primarily  for  checkout  purposes.  The  data 
printed  consist  of  either  input  data  to  one  module  from  the  generated  file 
data  of  a  preceding  module  or  the  matrix  equations  formed  within  the  module 
in  execution.  The  equations  are  described  in  the  module  description  section 
of  the  report.  The  user  should  note  that  the  use  of  these  print  flags  will 
result  in  a  large  amount  of  print. 

Output  from  the  aerodynamic  module  is  described  in  detail  in  Volume  II 
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The  user  should  check  the  aerodynamic  output  at  a  reduced  frequency  of 
zero  to  insure  that  the  generalized  forces  and  moments  due  to  rigid  body 
pitch  and  yaw  for  one  radian  rotations  are  numerically  identical  to  the 
gust  forces  for  the  vertical  and  lateral  gust  orientations  respectively. 

Lack  of  agreement  indicates  an  error  in  modal  deflection  inputs  for  the  aero 
node  points.  The  error  will  usually  be  in  the  sign  of  the  input  for  a 
group  of  points. 

Output  from  the  inertial  module  consists  of  the  mass  and  mass  point 
geometry,  the  inertial  mode  shapes,  the  modal  definition  table,  the  structur¬ 
al  damping  coefficients  and  the  generalized  mass  and  stiffness. 

Output  from  the  unit  loads  module  consists  of  a  summary  of  all  input 
data  for  the  module.  In  addition,  the  sweep  and  dihedral  for  each  local  beam 
are  printed  and  an  array  defining  the  beams  loaded  (by  input  number)  by 
any  load  first  entering  each  defined  beam.  The  TLAMY  matrices  and  TLAMM  for 
each  beam  are  displayed.  The  printout  following  consists  of  the  AAS  location 
of  each  load  point,  inertial  and  aero,  by  groups  as  defined  in  the  input. 

If  the  print  flag  KPRLDS  (Fixed  Data  Deck)  has  been  input  non-zero,  the 
integrated  loads  due  to  each  inertial  and  aero  load  point  will  be  printed 
out  along  with  the  integrated  loads  due  to  unit  modal  amplitude  for  inertial 
motion  and  aero  motion. 

The  active  system  module  output  consists  of  the  input  definition  data, 
a  summary  of  the  values  of  all  block  polynomials  used  to  define  the  system 
and  the  transfer  function  polynomial  in  powers  of  iu>.  The  final  transfer 
functions  are  also  printed  out  for  a  frequency  of  one  rad/sec  for  check 
purposes.  The  antisymmetric  data  follow  the  symmetric  data. 
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The  frequency  response  module  printed  output  consists  only  of  the 
gust  orientation  number,  direction  cosines,  the  fliqht  condition,  and 
the  frequencies  of  solution.  If  the  print  flag  KPRCXQ  has  been  input  for 
an  orientation  there  will  be  printed  the  generalized  response  solution 
vector  q  for  each  frequency.  If  the  print  flag  KPRCHQ  is  input,  the  file 
input  generalized  motion  dependent  aero  matrices  and  the  gust  forces  on  all 
aerodynamic  elements  will  be  printed  out.  The  aero  force  integration  arrays 
will  also  be  printed.  In  addition,  for  each  frequency  of  solution,  the 
interpolated  generalized  gust  forces,  motion  dependent  aero  forces  and  the 
equations  for  the  generalized  response  will  be  printed.  If  an  active 
system  has  been  defined,  all  incremental  mass  and  aero  matrices  will  also 
be  printed.  The  symmetric  cases  print  first,  and  are  followed  by  the 
antisymmetric  (if  an  antisymmetric  solution  is  required  for  the  particular 
orientation  flagged  for  print). 

The  basic  print  output  of  the  unit  gust  loads  module  consists  only  of 
a  message  of  the  number  of  orientations  analyzed.  If  the  plot  flag  PLQ 
has  been  specified  nonzero,  the  generalized  response  for  all  modes  for 
the  gust  orientation  specified  will  be  plotted.  If  the  plot  flag  PLL  is 
input,  the  integrated  load  so  specified  will  be  plotted  versus  frequency 
for  all  orientations.  If  the  print  flag  KPRCXL  is  nonzero,  there  will  be 
displayed  the  complex  integrated  loads  per  unit  gust  for  the  specified 
orientation  and  their  moduli  for  all  frequencies.  Settinq  the  flag 
KPRCXL  greater  than  100  will  cause  a  print  of  the  integrated  load  moduli 
for  all  orientations.  The  use  of  the  flag  KPRCHL  will  invoke  a  printout 
of  all  data  used  to  generate  the  loads  due  to  unit  gust  for  the  specified 
orientation.  These  data  consist  of  the  frequencies  of  solution,  the 
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interpolation  coefficients,  the  generalized  response  vectors  for  all  fre¬ 
quencies,  the  integrated  loads  due  to  unit  inertial  modal  amplitude  and 
unit  aero  modal  amplitude  (PIQ,  PAQS,  PAQA)  and  finally  the  integrated 
loads  from  summation  of  all  forces. 

The  output  from  the  Trim,  Blast  and  Time  History  Modules  consists  of 
a  summary  of  input  control  data  for  these  modules  and  a  summary  of  the  modal 
control  data  for  the  symmetric  trim  solution  and  flight  condition  data.  The 
generalized  coordinate  solutions  for  the  symmetric  part  of  the  solution  are 
displayed  and  the  symmetric  trim  parameters  for  the  maneuver  condition 
follow.  If  the  maneuver  specified  requires  an  antisymmetric  trim  solution, 
similar  data  for  the  antisymmetric  solutions  are  printed.  These  data  are 
followed  by  a  summary  of  symmetric  (and  antisymmetric,  if  present)  trim 
loads  at  the  specified  integrated  load  stations. 

In  the  event  that  the  user  has  flagged  the  trim  matrix  print  flag 
(KPRTRM)  there  will  be  displayed  the  aerodynamic  and  inertial  matrices  used 
in  the  trim  solution  and  the  matrix  equations  solved  for  the  symmetric  and 
antisymmetric  trim  solutions. 

The  amount  of  blast  and  time  history  output  is  controlled  by  the  three 
print  flags,  KPRTMH,  KPITMH,  and  KPRBLS.  The  minimum  print,  which  should 
be  the  case  for  range  iteration  runs,  consists  of  a  summary  of  blast  conditions 
and  the  delta  time  of  intercept  with  respect  to  the  AAS  origin  (x)  of  the 
blast  wave  with  the  first  aerodynamic  load  producing  element  of  the  vehicle 
and  its  AAS  coordinates.  These  are  followed  by  the  EFAS  coordinates  of  the 
AAS  origin  at  time  of  intercept,  the  EFAS  coordinates  of  the  burst,  the  time 
required  for  the  blast  wave  to  reach  the  vehicle  and  the  estimated  highest 
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frequency  of  the  gust  function.  The  output  following  consists  of  a  summary 
of  the  maximum  positive  and  negative  loads  experienced  at  all  integrated 
load  stations  and  the  maximum  allowables.  If  a  stress  matrix  has  been 
specified,  similar  data  are  displayed  for  the  stresses. 

The  next  printed  output  is  a  summary  of  the  present  blast  conditions 
and  the  estimated  allowable  blast  conditions  which  are  based  on  the  allowables 
for  load  or  stress.  The  next  output  (for  the  current  iteration)  consists 
of  an  estimate  of  the  EFAS  location  of  the  aircraft  for  the  estimated  slant 
range,  based  on  the  allowables,  and  the  relative  position  of  the  aircraft 
with  respect  to  the  burst  at  time  of  burst  and,  if  a  converged  solution 
has  been  achieved,  a  notation  of  same. 

The  final  output  for  each  orientation  is  a  summary  of  the  integrated 
load  definitions. 

If  the  user  has  flagged  the  blast  plot  flag,  KPLTMH ,  there  will  be 
plotted  the  material  velocity  time  history  and  the  time  history  of  all  the 
integrated  loads. 

If  the  print  flag  KPRTMH  has  been  input,  complete  load  time  histories, 
will  be  printed  for  all  integrated  loads.  If  the  print  flag  KPRBLS  has  been 
input  there  will  result  a  complete  printout  of  the  gust,  density  and  range 
time  history,  the  Fourier  transform  of  the  gust  material  velocity  times 
density  ratio  and  the  symmetric  (and  antisymmetric  if  necessary  for  the 
blast  condition)  frequency  response  function  which  consists  of  the  product 
of  the  frequency  response  of  the  aircraft  to  the  unit  gust  with  the  above 
transform  function.  The  complete  symmetric  time  history  perturbation 
(and  antisymmetric  if  required)  solution  for  all  integrated  loads  will 
also  be  printed. 
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SECTION  VI 


PROGRAM  OPERATION 

VIBRA-6  has  been  coded  to  operate  on  a  CDC  6600.  The  concept  of 
dynamic  core  has  been  utilized  throughout  the  program,  thus  removing  any 
restriction  on  the  problem  size  other  than  the  maximum  core  available. 

Dynamic  core  on  the  CDC  is  achieved  by  using  blank  COMMON  for  all 
dimensioned  variables.  Since  blank  COMMON  is  loaded  in  central  memory 
following  the  last  subroutine,  the  remainder  of  core  is  available  for 
use  as  dynamically  allocated  core.  An  RFL  control  card  specifying  the 
total  central  memory  (program  +  dynamic  core)  must  be  specified  before 
execution  because  blank  COMMON  has  been  dimensioned  so  that  it  has  a 
length  of  2.  Dimensioning  data  are  available  from  input  data,  the  length 
of  each  array  is  calculated,  and  the  starting  location  in  COMMON  of  that 
array  is  referenced  to  the  first  location  of  COMMON  as  a  subscripted  array. 
The  array  dimensioning  for  each  of  the  main  program  modules  is  shown  in 
Tables  4  through  13.  The  tables  show  the  variable  length  of  each  array 
and  the  relative  location  of  that  array  with  respect  to  the  other  arrays. 
Figure  22  shows  graphically  the  relationship  between  the  primary  program 
modules  and  the  various  dimension  levels.  The  total  amount  of  core  required 
to  run  any  specific  module  can  be  calculated  by  using  Table  14. 

The  various  COMMON  blocks  used  throughout  the  program  are  described 
in  detail  in  Tables  15  through  19.  The  tables  are  self-explanatory 
except  for  Table  19.  The  COMMON  block  DDTBLS  is  essentially  10  tables 
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of  20  words  each,  where  the  second  subscript  of  DDTBL  is  used  to  specify 
a  given  table.  The  tables  are  used  by  the  data  file  access  subroutines 
to  read  and  write  the  appropriate  data  from  the  generated  data  files. 

The  organization  of  the  data  on  the  files  is  described  in  Tables  20  through 
24. 


The  Fortran  listing  of  all  subroutines  is  given  in  Volume  III  of 
this  report. 

The  program  is  modularized  such  that  each  module  can  be  executed 
independently,  thus  an  efficient  overlay  feature  is  used  to  reduce  the 
amount  of  central  memory  necessary  for  any  given  run.  This  overlay  is 
accomplished  using  the  CDC  SEGLOAD  loader.  The  SEGLOAD  directives 
required  are  listed  in  Table  25. 
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TABLE  « 


CONTROL  MODULE  (DACGOST)  ARFAY  DIMENSIONING 


DIM 

LEVEL 

FEL 

LOC 

FFOG 

ID 

LENGTH 

APPAY 

NAME 

DEFINITION 

MDL1 

1 

LI 

20 

AMODNO 

MODE  LOCATION  DEFINEP 

2 

L2 

NOMOD 

WP 

MODAL  FREQUENCIES  (Hz) 

3 

L3 

NOMOL.NOMOD 

EMBAP 

GENERALIZED  MASS 

4 

L4 

NOMOL 

EMWR2 

GENERALIZED  STIFFNESS 

MC 12 

5 

L4  1 

NOMOC 

CPXDPG 

STRUCTURAL  DAMPING 

6 

L9 

NMS -NOMOD 

PHI  X 

X“ I NEPTIAL  MODE  SHAPES 

7 

Li  0 

NMS- NOMOD 

PHI  Y 

Y- INERTIAL  MODE  SHAPES 

fc 

Li  1 

NMS. NOMOD 

PHIZ 

Z- INERTIAL  MODE  SHAPES 

9 

L51 

NTFS. 4 

TMSS 

ACTIVE  SYSTEM  KINEMATICS 
SYMMETRIC 

10 

L5  2 

NTFA-4 

TMSA 

ACTIVE  SYSTEM  KINEMATICS 
ANTISYMMETRIC 

11 

L53 

2 .NTFS* MX OPE 

TFCS 

TRANSFER  FUNCTION  POLYNOMIAL 
SYMMETRIC 

12 

L54 

2  .NTFA-MXORD 

TFCA 

TRANSFER  FUNCTION  POLYNOMIAL 
ANTISYMMETRIC 

13 

L55 

NTFS-NSYM 

PHISS 

MODAL  AMPLITUDES  SENSED  (SYM) 

1  4 

L56 

NTFA-NASYM 

PHISA 

MODAL  AMPLITUDES  SENSED  (ASM) 

MDL3 

15 

L5 

IMS. NMS 

EM 

MASS  PROPERTY  TABLE 

16 

L6 

NMS 

ELXIO 

XAAS  COORDINATES  OF  MASSES 

17 

L7 

NMS 

ELY  10 

YAAS  COORDINATES  OF  MASSES 

18 

L8 

NMS 

ELZIO 

ZAAS  COORDINATES  OF  MASSES 

14* 

LI  2 

NCF-NMS. NOMOD 

EMPHI 

MASS  TIMES  MODE  SHAPES 

20* 

LI  3 

2  •  N  E  NG  S 

NEMGM 

ENGINE  THRUST  MASSES 

21** 

LI  4 

NCF.NEf.NMS 

EMS 

SECTIONAL  MASS  MATRIX 

MDL4 

15 

L60 

LENGTH 1 

D 

SCRATCH 

15 

L61 

2  •  NTCTAP*  NG 

F 

GUST  FOPCE 

16 

L62 

LENGTH 2 

WOP  K 

SCRATCH 

16 

L6  3 

NTOTAP.NSYM 

SPLS 

AERO  FOPCE  INTEGRATION  MATRIX 
SYMMETRIC 

17 

L64 

NTOTAP.NASYM 

SPLA 

AERO  FORCE  INTEGRATION  MATFIX 
ANTISYMMETRIC 

15 

L7C 

6.NECX 

GEOMBX 

BOX  AAS  COORDINATES 

16 

L71 

6.NSPETO 

GEOMED 

BODY  AAS  COORDINATES 

TABLE  4  (CONT*  D) 


CONTROL  MODULE  (DACGUST)  ARRAY  DIMENSIONING 


DEFINITIONS: 


NOMOD 

NTFS 

NTFA 

MXORL 

MXOBLK 

MXPLK 

NSYM 

NASYM 

IMS 

NM£ 

NDF 

NENGS 

NBCX 

NSBETO 

NTOTAP 

NG 

LENGTH1 

LENGTH2 


TOTAL  NUMEEF  OF  MODES 

NUMBER  OF  SYMMETRIC  TRANSFEF  FUNCTIONS 
NUMBEF  OF  ANTISYMMETRIC  TRANSFER  FUNCTIONS 
{MXOBLK-fl )  (MXBLK) 

MAX.  ORDER  OF  LARGEST  TRANSFER  FUNCTION  BLOCK 

NUMBER  OF  ELOCKS  IN  LARGEST  TRANSFER  FUNCTION 

NUMBER  OF  SYMMETRIC  MODES 

NUMBEF  OF  ANTISYMMETRIC  MODES 

NUMBER  OF  ITEMS  IN  MA.SS  PROPERTY  TABLE 

NUMBER  OF  MASSES 

NUMBER  OF  DEGREES  OF  FREEDOM  PER  MASS  STATION 

NUMBEF  OF  ENGINES 

NUMBEF  OF  PANEL  BOXES 

NUMBEF  OF  SLENDER  EODY  ELEMENTS 

NUMBEF  OF  TOTAL  AERO  ELEMENTS  (NBOX+2 -NSBETO) 

NUMBER  OF  GUST  ORIENTATIONS 
MAX (2-NSYM- NSYM,  2- NASYM. NASYM, 

NTOTAP-NSYM  +  NTO TAP -NASYM,  2 -NTOTAP-NG) 

MAX (2-NEOX-NG,  2-NSBETO-NG,  2-NSYM-NSYM,  2 -NASYM- NAS YM) 


*  Dimensioned  only  for  unit  loads  calculations 

**  Dimensioned  only  if  sectional  data  are  input 
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TABLE  5 


AEKODY  NAMIC  MODULE  ARRAY  DIMENSIONING 


DIM 

LEVEL 

PEL 

LOC 

PFOG 

ID 

LENGTH 

ARRAY 

NA*ME 

DEFINITION 

ADL1 

1 

LI 

NODES 

ELXI 

XAAS  OF  NODAL  POINT 

2 

L2 

NODES 

ELYI 

YAAS  CF  NODAL  POINT 

3 

Li 

NODES 

ELZI 

ZAAS  OF  NODAL  POINT 

4 

L4 

NODES. MODES 

PHINA 

MODE  SHAPE  IN  NOFMAL  DIR. 

5 

L5 

NCDES-MODES 

PHIZA 

MODE  SHAPE  IN  Z-DIRECTION 

6 

L6 

NODES. MODES 

PHI  YA 

MODE  SHAPE  IN  Y-DIPECTION 

7 

LI 

NF 

NAS 

NO.  OF  ASSOCIATED  EODIES  PER 
PANEL 

8 

L8 

NP.NE 

NASB 

ASSOCIATED  BODY  NUMBERS 

9 

L9 

NF 

N BAR AY 

LAST  BOX  NO.  OF  EACH  PANEL 

10 

Li  0 

NP 

NCARAY 

NO.  OF  CHOPDWISE  BOXES  PEP 
PANEL 

11 

Li  1 

NP 

NSARAY 

NO.  OF  SPANWISE  BOXES  PER 
PANEL 

12 

LI  2 

MSTRIP 

ISSTR 

SUPEPSTRIP  NO.  OF  EACH  STRIP 

13 

LI  3 

MSTP IP 

NSSTF 

NO.  OF  STRIPS  PEP  SUPERSTFIP 

14 

LI  4 

2.NE 

I  FLA 

SEQUENCE  NOS.  OF  IBE'S 

15 

L 1  5 

2  -  NP 

NEEA 

IBE  FLAGS 

16 

L 1 6 

2.NE 

NT  12 

NO.  OF  ANGULAR  ORIENTATIONS 
PER  BODY 

17 

LI  7 

NE 

NSBEA 

NO.  OF  SBE'  S  PER  BODY 

IB 

Li  8 

NE 

YE 

YAAS  CF  BODY  CENTEPLINES 

19 

LI  9 

NE 

ZB 

ZAAS  OF  BODY  CENTERLINES 

20 

L20 

NE 

APE 

BODY  CROSS-SECTIONAL  ASPECT 
RATIO 

21 

L21 

NB 

AVF 

BODY  AVE.  CHARACTERISTIC 

HALF- WIDTH 

22 

L22 

NE 

XLE 

XAAS  OF  LEADING  EDGE  OF  BODY 

1 3 

L23 

NE 

XTE 

XAAS  OF  TRAILING  EDGE  OF  BODY 

24 

L24 

MBE 

RIA 

RADII  OF  IBE'S 

25 

L25 

NB-10 

TH1 A 

TH1  VALUES  FOP  ALL  BODIES 

26 

L27 

MSPE 

AO 

SBE  HALF- WIDTHS 

27 

L28 

MSEE 

AOP 

X-DEFIVATIVE  OF  THE  AO  VALUES 

2  B 

L29 

MSEE 

XIS1 

XAAS  OF  LEADING  EDGE  OF  SBE'S 

29 

L30 

MSEE 

XIS2 

XAAS  OF  TRAILING  EDGE  OF 

SBE'S 

30 

L31 

MSTRIP+NB 

CG 

COSINE  OF  DIHEDPAL  OF  STRIPS 

J1 

L32 

MSTP IP+NP 

cs 

CHORD  LENGTH  OF  STRIPS 

32 

L33 

MSTRIP+NB 

EE 

HALF- WIDTH  OF  STRIPS 

53 

L34 

MSTRIP+NB 

SG 

SINE  OF  DIHEDRAL  OF  STRIPS 

34 

L35 

MSTRIP+NB 

YS 

YAAS  OF  STRIP  CENTERLINE  AND 
BODY  ELEMENTS 

35 

L36 

MSTRIP+NB 

ZS 

ZAAS  OF  STRIP  CENTERLINE  AND 
BODY  ELEMENTS 
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DIN 

LEVEL 

PEL 

LOC 

PROG 

ID 

LENGTH 

ARRAY 

NAME 

DEFINITION 

ADL 1 
COLT. 

36 

L37 

MSTRIP+NE 

XIJ 

XAAS  OF  LEADING  EDGE  OF  STRIP 
CENTERLINES 

37 

L38 

MSTRIP-fNB 

YIN 

YAAS  OF  INBOARD  EDGE  OF  PANEL 

38 

L39 

MSTP  IP-fNB 

ZIN 

ZAAS  OF  INBOARD  EDGE  OF  PANEL 

39 

L40 

MSTR  IP-fNB 

COORD 

SPANK ISE  COOFDINATE  OF  STRIPS 

UO 

L41 

NTOTAL 

X 

XAAS  OF  3/4  CHORD  OF  BOXES 

AND  MIDPOINT  OF  IBE'S 

41 

L42 

NTOTAL 

XIC 

XAAS  OF  1/4  CHOFD  OF  BOXES 

42 

L43 

NTOTAL 

CELX 

AVE.  CHORD-LENGTHS  OF  BOXES 
AND  IBE'S 

43 

L44 

NTOTAL 

XLAM 

TANGENT  OF  SWEEP  OF  1/4  CHORD 
OF  BOXES 

44 

L4  5 

NTOTAL 

H 

A  COLUMN  CF  h  MATRIX 

45 

L46 

NTOTAL 

CHDX 

A  COLUMN  OF  dh/dx  MATFIX 

46 

L47 

NTOTAL 

DELA 

AREAS  OF  BOXES 

46  A 

L4  7A 

4.  NTOTAL 

XC 

X-COOPDINATES  OF  CORNERS  OF 
LIFTING  SURFACE  BOXES 

47 

L48 

2 -NTOTAL 

CT 

A  ROW  OF  THE  DT  MATFIX 

48 

L49 

2 -NTOTAL 

DTA 

A  ROW  OF  THE  DTA  MATRIX 

49 

L50 

2-  (MCDES-fNG) 

PHS 

TEMPORARY  ARRAY 

A  CL  2 

50 

LI 

NCMAX 

TH 

PANEL  CHORDI WISE  DIVISIONS 

51 

L2 

NSMAX 

TAU 

PANEL  SPANWISE  DIVISIONS 

52 

L3 

NEE 

XII 

XAAS  OF  IBE'S 

53 

L4 

MBF 

PI 

SEMI-WIDTHS  OF  IBE'S 

54 

L5 

MSDE 

XIS 

XAAS  OF  SBE'S 

55 

L6 

MSBE 

RS 

SEMI-WIDTHS  OF  SBE'S 

56 

L7 

NP 

GMA 

DIHEDRAL  ANGLES  OF  PANELS 

57 

L8 

MSTRIP-fNB 

DYS 

dy  OF  PANEL  STRIPS 

58 

L9 

MSTR  IP-fNB 

DZS 

dz  CF  PANEL  STRIPS 

59 

LI  0 

MSTRIP+NB 

GMAR 

DIHEDRAL  ANGLE  OF  STRIPS 

60 

Li  1 

NTOTAL 

XII 

XAAS  OF  INBOARD  OF  1/4  CHOFD 
LOCATION  OF  EOXES 

61 

LI  2 

NTOTAL 

XI 2 

XAAS  OF  OUTBOARD  OF  1/4  CHORD 
LOCATIONS  OF  BOXES 

62 

Li  3 

NTOTAL 

ETA1 

YAAS  OF  INBOARD  OF  EOXES  AND 
CENTERLINE  OF  SBE'S 

t>3 

L14 

NTOTAL 

ETA  2 

YAAS  OF  OUTBOARD  OF  BOXES  AND 
CENTERLINE  OF  SBE'S 

64 

LI  5 

NTOTAL 

ZETA1 

ZAAS  OF  INBOARD  CF  EOXES  AND 
CENTERLINE  OF  SEE'S 

65 

Li  6 

NTOTAL 

ZETA2 

ZAAS  OF  OUTBOARD  OF  BOXES  AND 
CENTERLINE  OF  SBE'S 

66 

LI  7 

N TOTAL 

ETA 

YAAS  CF  CENTERLINE  OF  3/4 
CHORD  LOCATION  OF  BOXES 
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TABLE  5  (CONT'D) 

AEFODYNANIC  MODULE  ARRAY  DIMENSIONING 


DIM 

LEVEL 

FEL 

LOC 

PFOG 

ID 

LENGTH 

AKPAY 

NAME 

DEFINITION 

67 

LI  8 

NTCTPL 

ZETA 

ZAAS  OF  CENTEFLINE  OF  3/4 
CHOPD  LOCATION  OF  BOXES 

68 

LI  9 

MSBE 

ETAS 

YAAS  OF  LEADING  EDGE  OF  SBE*  S 

69 

L20 

MSBE 

ZETAS 

ZAAS  OF  LEADING  EDGE  OF  SBE'S 

ALL  3 

5C 

N1 

NE 

NBCOM 

SBE  INDEXES 

51 

N2 

MAX  (NEfNP) 

INSUP 

BODY/PANEL  NUMBERS 

52 

N3 

NODES 

NODE 

NODE  NUMBEFS  IN  SUPEFEODY/ 
SUPER  PANEL 

53 

N  4 

NODES 

XIC 

XAAS  OF  NODAL  POINTS  IN 
SUPEFBODY/SUPEPPANEL 

54 

N5 

NODES 

ETAQ 

YAAS  OF  NODAL  POINTS  IN 
SUPEFBODY/SUPEPPANEL 

55 

N6 

MSTF IP4NP 

CSG 

COSINES  OF  MODIFIED  DIHEDRALS 
OF  STRIPS 

56 

N8 

NTOTAL 

X3L 

SBE  ENDPOINTS  PER  SUPEFBODY 

57 

N9 

NTCT7L 

YL 

YAAS  OF  ELEMENTS  IN  SUPEFBODY 

58 

N10 

NTOTAL 

XP 

XAAS  OF  BOXES  IN  SUPERPANEL 

59 

N1  1 

NTOTAL 

YP 

YAAS  OF  BOXES  IN  SUPEF PANEL 

60 

N12 

NMAX-NMAX 

XKD 

WORK  ARRAY  FOP  SPLINE 

61 

N13 

NMAX-NMAX 

FHS 

WOFK  ARRAY  FOR  SPLINE 

ADL4 

50 

LI 

NTMAX 

COL 

COLUMN  OF  h  OP  dh/dx 

51 

L2 

NTMAX.MODE 

WORK 

WOFK  ARRAY  USED  BY  ORGN 

A  CI.  5 

50 

LI 

2.NTCTAL 

DT 

ROW  OF  DT  MATRIX 

50 

LI 

2-NTCTAL 

CPZ 

ROW  OF  DPZ  MATRIX 

51 

L2 

2.NTCTAL 

CP  Y 

ROW  OF  DPY  MATRIX 

52 

L3 

2-NTCTAL 

CTA 

ROW  OF  DTA  MATRIX 

52 

L3 

2-NTCTAL 

CPZA 

ROW  OF  DPZA  MATRIX 

53 

L4 

2-NTCTAL 

DPYA 

ROW  CF  DPYA  MATRIX 

54 

L5 

1 12-NCMAX 

WOFK 

WORK  ARRAY  FOR  GEND 

A  Cl.  6 

50 

LI 

4 -NTOT -MSBE+ 

6-NTCT4 

10-MSEE 

WOFK 

WORK  AREA  FOF  SB 

ADL7 

50 

LI 

4- LENGTH 

BFS 

ROW  OF  SYMMETRIC  BFS  MATRIX 

51 

L2 

4.  LENGTH 

EFSA 

ROW  OF  ANTI -SYMMETRIC  BFS 
MATRIX 

ADL8 

50 

LI 

2-NTCTAL 

W 

COLUMN  OF  W  MATFIX 

51 

L2 

2- NTOTAL 

DM 

COLUMN  OF  DW  MATRIX 

52 

L3 

2-NTCTAL 

COL 

COLUMN  FOR  WORK  STORAGE 
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TABLE  5  (CONT'D) 


AERODYNAMIC  MODULE  ARFAY  DIMENSIONING 


DIM 

LEVEL 

PEL 

LOC 

PROG 

ID 

LENGTH 

ARRAY 

NAME 

DEFINITION 

A  CL  9 

50 

LI 

2-NTCTAL 

WGS 

COLUMN  OF  SYMMETRIC  GUST 
BOUNDARY  CONDITIONS 

51 

L2 

2  •  NTCTAL 

WGA 

COLUMN  OF  ANTI-SYMM2TFIC  GUST 
BOUNDARY  CONDITIONS 

ADL10 

50 

LI 

2-NTCT-NMODE 

WORK 

WORK  AREA  FOR  FHSIDE 

ADL1 1 

50 

LI 

2-MAX (3-NPM, 
NTOT-NMODE) 

WORK 

WORK  AREA  FOP  SOLVIT 

ADL12 

50 

LI 

2-LENGTH 

DT 

WORK  ARFAY  FOR  MATMUL 

51 

L2 

4-LENGTH-MSBE 

WOFK 

WORK  ARRAY  FOP  MATMUL 

A  ELI  3 

50 

LI 

2-NTMAX 

FORCE 

POINT  FORCE  COLUMN 

51 

L2 

2-NTMAX 

GF 

GENERALIZED  FORCE  COLUMN 

52 

L3 

2-NTMAX 

WORK 

WORK  ARRAY 

ADL 1  4 

50 

LI 

2-NTOTAL 

DCP 

NODAL  POINT  PRESSURES 

51 

L2 

2-NRMAX 

FZ 

BODY  ELEMENT  z-FOFCES 

52 

L3 

2-NPMAX 

FY 

BODY  ELEMENT  y-FOFCES 

53 

L4 

2-NRMAX 

CM 

STRIP  LIFT  COEFFICIENTS 

54 

L5 

2-NFMAX 

CM 

STRIP  MOMENT  COEFFICIENTS 

55 

L6 

2-NRMAX 

SPLD 

SPAN  LOADS 

56 

L7 

2-NB 

CZB 

BODY  LIFT  COEFFICIENT  IN 
Z-DIFECTION 

57 

L3 

2-NE 

CYB 

BODY  LIFT  COEFFICIENT  IN 
y-DIRECTION 

58 

L9 

2-NP 

CNB 

BODY  YAWING  MOMENT  COEFF. 

59 

L10 

2-NB 

CMB 

BODY  PITCHING  MOMENT  COEFF. 

60 

L1 1 

MSTPIP 

CPR 

STRIP  CENTER  OF  PRESSURE 
(REAL  PART) 

61 

LI  2 

MSTRIP 

CPI 

STRIP  CENTER  OF  PRESSURE 

(IMAGINARY  PART) 
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TAELE  5  (CONT* D) 

AERODYNAMIC  MODULE  ARRAY  DIMENSIONING 


DEFINITIONS: 

SEE  SLENDER  BODY  ELEMENT 

I  BE  INTERFERENCE  BODY  ELEMENT 

NODES  NUMBER  OF  AERODYNAMIC  NODAL  POINTS 

MOLES  NUMBER  OF  MODES  =  NSYM+NASYM 

MODE  MAX (NSYM,NASYM) 

NMODE  MODE+NG 

NP  TOTAL  NUMEER  OF  PANELS  ON  ALL  LIFTING  SURFACES 

NE  TOTAL  NUMEER  OF  EODIES 

MSTFIF  TOTAL  NUMEER  OF  STRIPS  FOP  ALL  EODIES 

MEL  TOTAL  NUMEEP  OF  IBE'S  FOR  ALL  BODIES  +  1 

MSEE  TOTAL  NUMEER  OF  SEE'S  FOR  ALL  BODIES  +  1 

NTCTAL  TOTAL  NUMEER  OF  LIFTING  SURFACE  BOXES  +  2 • MSBE 

NTOT  TOTAL  NUMEEF  OF  LIFTING  SURFACE  BOXES  +  2-MEE 

N’PM  NTOT+NMCDE 

NMA.X  NTOTAL+3 

NCMAX  MAX.  NC.  CF  CHORDS ISE  BOXES  PER  PANEL 

NSMA.X  MAX.  NO.  CF  STRIPS  PEF  PANEL 

NTMAX  MAX (NTCTAI , 2 -MEE) 

LENGTH  N TOTAL  +  2-MBE 

NFMAX  MAX  (MSTRI i, MSBE,MOD£S+2 -NG) 

NG  NUMEER  OF  GUST  ORIENTATIONS 
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TAELE  6 


TAELE  6  (CONT'D) 

UNIT  LOADS  MODULE  ARFAY  DIMENSIONING 


DIM 

LEVEL 

FEL 

LOC 

PFOG 

ID 

LENGTH 

ARRAY 

NAME 

DEFINITION 

24 

L22 

NSEETO 

ETAS  2 

SBE  OUTER  YAAS 

25 

L23 

NSBEFO 

ZETAS2 

SBE  OUTET  ZAAS 

LDL6 

17 

Li  5 

MXFOW -MX COL 

FPH 

BOX  OF  BODY  FORCE  DUE  TO 

MODAL  AMPLITUDE 

18 

Li  6 

NINTLD- MX COL 

FAO 

INTG*  D  LOAD  DUE  TO  MODAL 
AMPLITUDE 

DEFINITIONS: 


SPF  SLENDER  BODY  ELEMENT 


N BEAMS 

NINTLD 

NMGEP 

NABGF P 

NS EG! P 

1,  ST  R  SS 

NMASF 

NBCXES 

NBC  X 

NAER  SI- 

LENGTH 

LENGTHA 

LENGTH B 

N SB ETC 

NENGS 

NSYM 

NASYM 

NM 

MXFOW 

MXCOI 


NUMBER  OF  LOAD  BEAMS 

NUMBEF  OF  INTEGRATED  LOADS 

NUMBER  OF  MASS  GROUPS 

NUMEEP  OF  AERO  EOX  GfOUPS 

NUMBER  OF  SLENDER  BODY  ELEMENT  GROUPS 

NUMEEP  OF  STRESSES 

MAX,  NO.  OF  MASSES  IN  THE  NMGEP  MASS  GROUPS 
MAX.  NO.  CF  AERO  BOXES  IN  THE  NABGRP  BOX  GFOUP 
NUMBEF  OF  AERO  BOXES  (TOTAL) 

MAX.  NO.  SLENDEF  BODY  ELEMENTS  IN  THE  NSEGF  P  BODY  GROUP 

NINTLD.NMGFtP.NMASS 

NINTLD -NAEGFP-N BOXES 

NINTLD- NSEGF P«NAEFSE 

NUMBEF  OF  SLENDER  BODY  ELEMENTS  (TOTAL) 

NUMBEF  OF  ENGINES 
NUMBEF  OF  SYMMETRIC  MODES 
NUMBEF  OF  ANTISYMMETRIC  MODES 
TOTAL  NUMEEP  OF  MODES  (NSYM+NASYM) 

MAX (NEOX , NS  EETO) 

2-MAX (NSYM, NASYM) 


* 


Dimensioned  only  if  sectional  data  are  input 
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TABLE  7 


ACTIVE  SYSTEM  MODULE  AFRAY  DIMENSIONING 


DIM 

PEL 

PROG 

ARRAY 

LEVEL 

LOC 

ID 

LENGTH 

NAME 

DEFINITION 

CDI1 

1 

LI 

NTFM -MXOFD 

ANUM 

TEMPORARY  ARRAYS  FOF 

2 

L2 

NTFM-MXOFD 

ADEN 

FORMING  TRANS FEF 

3 

L3 

2.MXCRC 

ATEN 

FUNCTION  POLYNOMIALS 

4 

L4 

2* MXCRD 

ATFC 

TFCS  AND  TFCA 

5 

L5 

MXOFD 

ANUMT 

6 

L6 

MXOFD 

ACENT 

DEFINITIONS: 

NTFM  MAX.  NO.  TFANSFEF  FUNCTIONS  (NTFS  OR  NTFA) 

MXBLK  NO.  BLOCKS  IN  LARGEST  TRANSFER  FUNCTION 

MXOBLK  MAX  ORDER  OF  LARGEST  FUNCTION  BLOCK 

MXCRD  MXBLK- (MXCBLK+1) 
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TABLE  8 


FREQUENCY  RESPONSE  MODULE  ARRAY  DIMENSIONING 


DIM 

LEVEL 

FEL 

LOC 

PROG 

ID 

LENGTH 

ARRAY 

NAME 

DEFINITION 

FDL  1 

1 

Li  2 

2 -NMMX. NMMX 

D 

GEN.  FORCE  DUE  TO  MOTION 

2 

LI  3 

2 -NMMX 

F 

GEN.  FORCE  DUE  TO  GUST 

3 

LIU 

NVE« 

OMEGA 

FREQUENCY  (FAD/SEC) 

u 

Li  5 

2-NM 

Q 

GEN.  RESPONSE 

5 

L18 

2 -NMMX -NMMX 

TFM 

SCRATCH  FOR  ACTIVE  SYSTEM 

6 

LI  9 

NTOT^P-NS YM 

SPLS 

AERO  FORCE  I NT EG R  MAT  (SYM) 

7 

L20 

NTOT7P -NASYM 

SPLA 

AERO  FORCE  INTEGF  MAT  (ASM) 

FDL2 

8 

N 1 S 

2-NSYM2-NK 

DSALL 

SYM.  GEN.  AERO  (ALL  NK) 

9 

N2S 

2 -NTCTAP-NK 

FSALL 

SYM.  GEN.  GUST  (ALL  NK) 

10 

N3S 

NTOTAP -NK 

FS 

SYM.  MODULUS  GUST  (ALL  NK) 

11 

N4S 

NTCTAP-NK 

TANS 

SYM.  PHASE  GUST  (ALL  NK) 

12 

N5S 

2 -NTOTAP 

FS 

SYM.  GUST  FORCE 

13 

N1A 

2-NASYM2-NK 

DAALL 

ASYM.  GEN.  AERO  (ALL  NK) 

1U 

N2A 

2 -NTCTAP-NK 

FAALL 

ASYM.  GEN.  GUST  (ALL  NK) 

15 

N3A 

NTOTAP -NK 

RA 

ASYM.  MODULUS  GUST  (ALL  NK) 

ie 

N4A 

NTOTAP -NK 

TANA 

ASYM.  PHASE  GUST  (ALL  NK) 

17 

N5A 

2 -NTOTAP 

FA 

ASYM.  GUST  FORCE 

DEFINITIONS: 


NMMX  MAX (NSYM.NASYM) 

NVEH  NO.  FFEOUENCIES  FOF  SOLUTION  (NFREQ) 

NM  TOTAL  NO.  MODES  (NSYM+NASYM) 

NK  NO.  HARD  FOINT  AEFO  MATRICES 

NTOTAP  TOTAL  NO.  AERO  ELEMENTS  (NBOX+2 -NSBETO) 

NSYM  NO.  SYM  MODES 

N^SYM  NO.  ASYM  MODES 

NSYM2  NSYM-NSYM 

NASYM2  NASYM-NASYM 
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TABLE  9 


UNIT  GOST  LOAD  MODULE  ARRAY  DIMENSIONING 


DIM 

PEI 

PROG 

ARRAY 

LEVEL 

LOC 

ID 

LENGTH 

NAME 

DEFINITION 

GCL 1 

1 

L60 

NFREC 

PLOTF 

PLOTTING  ARRAY 

2 

L61 

NFREC 

PLOT  V 

PLOTTING  ARFAY 

3 

LI 

NINTLD. NFREC 

PS 

SYMMETRIC  UNIT  GUST  LOADS 

4 

L2 

NINTLD. NFREC 

PA 

ANTISYMMETRIC  UNIT  GUST  LOADS 

5 

L3 

NFREC 

OMEGA 

FREOOENCIES  (FAD/SEC) 

5  A 

L3A 

2.NACC 

INDACC 

ACCELERATION  INDICES 

GDL2 

6 

L4 

NK- NFREC 

COEF 

INTERPOLATION  COEFFICIENTS 

7 

L5 

2. NM. NFREC 

0 

GEN.  RESPONSES 

GDL3 

8 

L6 

NINTLD. NM 

PIC 

INERTIA  FORCES  DUE  TO  GEN. 
RESPONSE 

GDLU 

b 

L7 

2 .NINTLD. KSYM 

FACS 

AEPO  FORCES  DUE  TO  SYM  RESP. 

9 

L8 

2- NINTLD. NASM 

PAQA 

AEPO  FORCES  DUE  TO  ASYM  RESP. 

GDI  5 

6 

L51 

NINTLD 

SYMCOD 

CENTERLINE  LOAD  MOD  (SYM) 

7 

L52 

NINTIE 

ASMCOD 

CENTERLINE  LOAD  MOD  (ASYM) 

8 

L9 

3.NAEGPP 

NFNLAE 

SEE  LOAD 

9 

LI  0 

3  .NSBGRP 

NFNLSE 

SEE  LOAD 

10 

Li  1 

LENGTHA 

PINTP 

SEE  LOAD 

1 1 

Li  2 

LENGTHE 

PINTZ 

SEE  LOAD 

12 

LI  3 

LENGTHE 

PIKTY 

SEE  LOAD 

GET.6 

13 

L62 

2-NTCTAP 

FGPS 

SYM  GUST  FORCES 

GDL7 

13 

L6  3 

2-NTCTAP 

FGPA 

ASYM  GUST  FORCES 

DEFINITIONS: 


NFREC 

NINTLD 

NK 

NM 

NS  YM 

N/-.SM 

NABGFP 

NSFGFP 

N BOXES 

NAERSE 

LENGTHS 

LENGTHB 

NTOTAP 


NO.  FRECOENCIES  FOP  SOLUTION 
NO.  INTEGFATED  LOADS 
NO.  HARD  POINT  AEPO  MATRICES 
NO.  TOTAL  MODES  (NSYM+NASM) 

NO.  SYMMETRIC  MODES 

NO.  ANTISYMMETRIC  MODES 

NO.  AEPO  EOX  GROUPS 

NO.  SLENDEF  EDDY  ELEMENT  GROUPS 

MAX  NO.  AERO  BOXES  IN  ANY  NABGPP  BOX  GROUP 

MAX  NO.  SIENDEF  BODY  ELEMENTS  IN  ANY  NSBGRP  BODY  GROUP 

NINTLD. KAEGFP.NBOXES 

NINTLD- NSBGRP. NAERSE 

TOTAL  NO.  AEFO  ELEMENTS  (NBOX+2 -NSBETO) 
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TABLE  10 


TP IM  MODULE  AFRAY  DIMENSIONING 


DIM 

LEVEL 

PEL 

LOC 

FROG 

ID 

LENGTH 

ARFAY 

NAME 

DEFINITION 

TCL1 

1 

LI 

NM 

0 

GEN.  STATIC  FESPONSE 

2 

L2 

NM 

QDOT 

GEN.  STATIC  FATE  FESPONSE 

3 

LI  4 

NINTLD-NENGS 

THFLOD 

INTGD  LOADS  DUE  TO  THRUST 

u 

L15 

NSTM* NENGS 

THFGNF 

GEN.  FOPCE  DUE  TO  THPUST 

TEL  2 

5 

L3 

NMS- NMS 

ZS  YM 

SYM.  EQNS.  STATIC  TPIM 

e 

L  4 

NMS 

FSYM 

RHS  EONS.  STATIC  TRIM 

7 

L5 

2- NSYM- NSYM 

DS  YM 

GEN.  SYM  AEPO 

TDL3 

5 

L6 

NMA. NMA 

ZASYM 

ASYM.  EQNS.  STATIC  TFIM 

6 

L7 

NMA 

FASYM 

FHS  EQNS.  STATIC  TPIM 

7 

L8 

2* NASM. NASM 

DASYM 

GEN.  ASYM  AEFO 

TDL4 

5 

Li  1 

NINTLD.NM 

Ply 

INTGD  INERTIAL  LOADS, 

INERTIAL  RESPONSE 

TDL5 

5 

LI  2 

2 • NINTLD. NS YM 

PAQS 

INTGD  SYM  LOADS,  AERO  GEN. 
RESPONSE 

6 

Li  3 

2  .NINTLD.  NASM. 

PAQA 

INTGD  ASM  LOADS,  AEPO  GEN. 
RESPONSE 

DEFINITIONS: 


NM  NO.  TOTAL  MODES  (NSYM+NASM) 

NSYM  NO.  SYM  MODES 

NASM  NO.  ASM  MODES 

NIKTLD  NO.  INTGD  LOADS 

NENGS  NO.  ENGINES 

NMS  2+NO.  ELASTIC  MODES  (SYM) 

NMA  3+NO.  OF  ELASTIC  MODES  (ASYM) 


TABLE  11 

BLAST  AND  TIME  HISTOPY  MODULE  ARFAY  DIMENSIONING 


DIM 

LEVEL 

PEL 

LOC 

PPOG 

ID 

LENGTH 

ARRAY 

NAME 

DEFINITION 

3C1 1 

1 

LI 

NINTLD 

PSTR 

SYM.  TRIM  LOADS 

2 

L2 

NINTLD 

PATR 

ASM.  TRIM  LOADS 

3 

L21 

NSTPfS. NINTLD 

STRESS 

STRESS/INTGD  LOAD 

u 

L22 

NENGS 

THRUST 

ENGINE  THRUST 

5 

L23 

NINTLD 

SYMCOD 

CENTER  LINE  LOAD  MOD.  (SYM) 

6 

L24 

NINTLD 

ASMCOD 

CENTER  LINE  LOAD  MOD.  (ASYM) 

BDL2 

7 

L3 

NOFMAX 

NOP 

EL AST  ORIENTATION  CODES 

3 

L4 

NORMAX 

PEST 

BLAST  ESTIMATED  PANGES 

9 

L5 

8. NINTLD 

STALDS 

INTEGRATED  LOAD  DEFINITION 

10 

L6 

NFFEC 

OMEGA 

FREQUENCIES  (RAD/SEC) 

11 

L7 

LENGTH 

PS 

SYM.  GUST  LOADS  (FREQ  DOMAIN) 

11 

L7 

NSTFES-NTMRSP 

STF 

RIGHT  SIDE  STF ESSES 

12 

L8 

LENGTH 

PA 

ASM.  GUST  LOADS  (FKEQ  DOMAIN) 

12 

L8 

NSTFSS.NTMPSP 

STL 

LEFT  SIDE  STRESSES 

1? 

L9 

NTMGST 

VELG 

GUST  VELOCITY 

14 

LI  0 

NTMGST 

TIMG 

GUST  VELOCITY  TIMES 

15 

Li  1 

NTMGST 

PHO 

GUST  DENSITY 

16 

Li  2 

NFFEC 

FLOTF 

PLOT  ARRAY 

17 

LI  3 

NFFEC 

PLOTA 

PLOT  ARRAY 

Ifl 

LIU 

2 -NFFEC 

FTG 

FOURIER  TRANSFORM  OF  VELG-RHO 

19 

Li  5 

NINTLD -NTMPSP 

PST 

RIGHT  SIDE  LOADS (TIME  DOMAIN) 

20 

LI  6 

NINTLD-NTMRSP 

PAT 

LEFT  SIDE  LOADS  (TIME  DOMAIN) 

21 

Li  7 

NTMPSP 

TIME 

LOAD  TIMES 

22 

Li  8 

MAXI 

PKL 

PEAK  LOADS  (OP  STRESSES) 

23 

Li  9 

NINTLD- 2 

ALLOWS 

MAX.  ALLOWABLE  STFESSES 

11 

L7 

6-MBCX 

GEOMBX 

BOX  AAS  COOFDINATES 

12 

L8 

6-NSEETO 

GEOMBD 

BODY  AAS  COOFDINATES 

DEFINITIONS: 

NINTLD 

NO. 

INTEGRATED  LOADS 

NSTFSS 

NO. 

STRESSES 

NENGS 

NO. 

ENGINES 

NOF MAX 

NO. 

OP IENTATIONS  TO  ANALYZE 

NFPEQ 

NO. 

FREQUENCY  SOLUTIONS 

NTMFSP 

NO. 

SOLUTION  TIME  POINTS 

NTMGST 

NO. 

GUST  HISTOPY  TIME  POINTS 

NBOX 

NO. 

PANEL  BOXES 

NSEETO 

NO. 

SLENDER  BODY  ELEMENTS 

LENGTH 

2- NINTLD- NFFEQ 

MAXI 

MAX (NINTLC- 8,NSTRSS-8) 
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TABLE  12 


RIGID  MODULE  ARRAY  DIMENSIONING 


DIM 

LEVEL 

FEL 

LOG 

FFOG 

ID 

LENGTH 

AFPAY 

NAME 

DEFINITION 

FDL 1 

1 

LI 

NFFEC 

OMEGA 

FFEQUENCIES 

2 

L2 

MAXI 

TIME 

TIME 

3 

L  3 

NFOFC 

GEOM 

FOFCE  GECMETFY 

4 

L4 

2.NFCPC 

CPT 

FOFCES  DUE  TO  TFIM 

5 

L5 

2-NFOFC.NFREQ 

CP 

FOFCES  AND  PFESSUFES  VS  FREQ 

KDL2 

6 

L6 

NBOX 

I BOXES 

AEFO  BOX  NUMBERS 

6A1 

L7 

2.NSEE 

ISBES 

SLENDER  BODY  ELEMENTS 

6A2 

L8 

NSXM 

Q 

GENERALIZED  COOF DI NATES 

6A3 

L20 

2.NFCFC.NK 

FI 

INPUT  AEFO  GUST  FOFCES 

6AU 

L21 

NK 

VOBUS 

PEDUCED  VELOCITIES  FFOM  AEFO 

6A5 

L22 

NFOFC- NK 

FS 

MODULUS  GUST 

6A6 

L23 

NFOFC* NK 

TANS 

PHASE  GUST 

6A5 

L9 

MEOX-5 

GEOP 

PANEL  GECMETFY 

6A5 

LI  0 

MSBE.6 

GEOS 

BODY  GEOMETF Y 

6A5 

L10A 

3.NOFMAX 

CP 

GUST  DIRECTION  COSINES 

6A5 

L 1 1 

2.MECX.MAX2 

C 

BOX  AERODYNAMICS 

6A5 

Li  2 

2.MSEE.MAX2 

EB 

SLENDER  BODY  AEFO 

FDL  3 

6 

LI  3 

NTMGST 

DP 

CVEPPRSSSURE 

6B1 

Li  4 

NTMGST 

FHO 

GUST  DENSITY 

6E2 

LI  5 

NTMGST 

VG 

MATERIAL  VELOCITY 

6E3 

LI  6 

NFFEC 

PLOTF 

PLOT  ARRAY 

6E4 

Li  7 

MAX  3 

PLOTA 

PLOT  ARRAY 

6E5 

L18 

2-NFFEQ 

FTG 

FOURIER  TRANSFORM  OF  VG-FHO 

6B6 

LI  9 

NFOFC. NTMRSP 

P 

FOFCES  6  PRESSURES  VS  TIME 

DEFINITIONS: 

NFFEQ 

NO. 

FREQUENCY  SOLUTIONS 

NTMGST 

NO. 

GUST  BISTOPY  TIME  POINTS 

NTMFSP 

NO. 

SOLUTION  TIME  POINTS 

NFOFC 

NO. 

OF 

FOFCES 

NPOX 

NO. 

OF 

AEFO  BOXES 

NSBE 

NO. 

OF 

SLENDER  BODY  POINTS 

NSYM 

NO. 

OF 

SYF METRIC  MODES 

NK 

NO. 

OF 

FCEUCED  FREQUENCIES  ON  AERO  FILE 

MBOX 

NO. 

OF 

AEFO  BOXES  ON  AERO  FILE 

MSBE 

NO. 

OF 

SLENDER  BODY  ELEMENTS  ON  AEFO  FILE 

NORMA X 

NO. 

OF 

ORIENTATIONS  TO  BE  ANALYZED 

NGUST 

NO. 

OF 

GUST  ORIENTATIONS  ON  AERO  FILE 

MAXI 

MAX (NTMGST, NTMFSP) 

MAX2 

MAX (NSYM, NGUST) 

MAX3 

MAX ( NF  REQ , NTMRS  P) 
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TABLE  13 


MERGE  MODULE  AFFAY  DIMENSIONING 


DIE. 

pel 

FROG 

AFFAY 

LEVEL 

1,0  c 

ID 

LENGTH 

NAME 

DEFINITION 

X  CL  1 

1 

LI 

I NK 1 

FKIN1 

FREQUENCIES  TO  EE  MERGED  FROM 
TAPE1 7  FILE 

2 

L2 

INK  2 

FKIN2 

FREQUENCIES  TO  BE  MERGED  FROM 
TAP E 1 8  FILE 

3 

L  3 

INK1+INK2 

PK 

DESIRED  FREQS  ON  TAPE19  FILE 

u 

L4 

NK1 

RK 1 

REDUCED  FREQS  ON  TAPE17  FILE 

5 

L5 

NK2 

FK2 

REDUCED  FREQS  ON  TAPE18  FILE 

b 

L6 

LENGTH 

WORK 

WORK  ARRAY 

DEFINITIONS 

;- 

INK  1 

NO. 

OF  FREQUENCIES 

TO 

BE  MERGED  FROM 

TAPE 17 

FILE 

INK  2 

NO. 

CF  FREQUENCIES 

TO 

EE  MERGED  FROM 

TAPE1  8 

FILE 

NKl 

NO. 

OF  FREQUENCIES 

ON 

TAPE17  FILE 

NK  2 

NO. 

OF  FREQUENCIES 

ON 

TAPE 18  FILE 

LENGTH 

6. MAX  (NEOX.NSPE) 

NBOX 

NO. 

OF  AERO  BOXES 

Nsr : 

NO. 

OF  SLENDER  EODY 

ELEMENTS 
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Core  Overlay  Structure 
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TABLE  14 

CENTRAL  MEMORY  REQUIREMENTS 


MODULE 

CENTRAL  MEMORY  (DECIMAL  WORDS) 

Program 

40,000 

AERO 

ADL1  +  MAX  (ADL6,  ADL  11,  ADL  12) 

UNIT 

MDLl  +  MDL2  +  MDL3  +  LDL1  +  LDL3  +  LDL6 

ACSM 

Other  modules  are  always  longer 

FRSP 

MDLl  +  MDL2  +  FDLl  +  FDL2 

GUST 

MDLl  +  MDL2  +  GDL3  +  GDL5  +  GDL6 

BLST 

MDLl  +  BDL1  +  BDL2 

RIGD 

RDLl  +  MAX  (RDL2,  RDL3) 

MERG 

XDL1 

Note:  If  more  than  one  module  is  to  be  executed  then  the 
amount  of  central  memory  required  is  that  necessary  for 
the  longest  module. 
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TAELE  15 


XTF  COMMON  FLOCK  DESCRIPTION 


oc 

Variable 

Description 

1 

A  MACH 

Mach  number 

- 

GAM  X 

direction  cosine  in  x-direction  for  current 
qust  orientation 

VSS 

speed  of  sound 

u 

P0 

atmospheric  pressure  at  altitude 

C 

NMS 

number  of  'right  hand  side*  mass  points, 
including  centerline  masses 

6 

NSY.M 

number  of  symmetric  modes 

7 

NAS  VM 

number  of  antisymmetric  modes 

b 

NOMOL 

total  number  of  modes,  NSYM  +  KASYM 

q 

NFREQ 

number  of  frequencies  in  frequency  response 

1C 

NENGS 

number  of  enqines  on  right  hand  side  and 
centerline 

11 

I  DENT 

identification  number  of  configuration 

12 

B  ZERO 

reference  semi-chord 

13 

RHOO 

density  at  sea  level 

14 

SIGMA 

density  ratio 

15 

EOUVAS 

equivalent  airspeed  factor 

16 

UNITS 

airspeed  units  flag;  =  2,  knots 

17 

VKEAS 

velocity,  equivalent  airspeed  in  knots 

18 

VTFPS 

velocity,  true  airspeed  in  ft/sec 

19 

ALT 

altitude,  ft 

20 

KPFCHK 

=nm,  print  all  matrices  in  frequency  response 
and  unit  qust  load  calculations  for 
orientations  n  and  m 

21 

KPRCXI 

=nm,  print  generalized  response  (q)  and 
integrated  loads  (P)  due  to  unit  gust  for 
orientations  n  and  m 

2 2 

KPPELS 

=1,  print  load  time  history  matrices 

23 

KPFTPM 

=1,  print  all  matrices  in  trim  solution 

24 

KPRLDS 

=1 ,  print  unit  load  matrices 

25 

KPFTMli 

=1,  print  total  load  time  history 

26 

IPLO 

=n,  plot  generalized  response  (q)  for  nth 
orientation 

27 

IPLL 

=n,  plot  nth  load  due  to  unit  gust  for  all 
orientations 

28 

IPL5L 

=1,  plot  load  time  histories  for  RHS  of  aircraft 
=2,  plot  lead  time  histories  for  LHS  of  aircraft 
=3,  plot  both  RHP  and  LHS 

29 

NACC 

number  of  accelerations 

31 

NEEAMS 

number  of  beams  for  integrated  loads 

32 

NINTLD 

number  of  integrated  loads 

33 

NSTRSS 

number  of  stresses 

34 

NMGRP 

number  of  groups  of  masses  associated  with 
beams 
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1 

_  ij 

TABLE  15  (CONT'D) 

X1F  COMMON  ELOCK  DESCRIPTION 

if 

1 

1  ■ 

1 

Loc 

Variable 

Description 

1 

35 

NA3GFP 

number  of  qroups  of  aerodynamic  boxes 
associated  with  beams 

i 

36 

NSBGKP 

number  of  groups  of  aerodynamic  slender  body 
elements  associated  with  beams 

37 

NBOXES 

naximum  number  of  aerodynamic  boxes  per  group 

38 

NAESSB 

maximum  number  of  aerodynamic  slender  body 
elements  per  group 

41 

NK 

number  of  reduced  frequencies 

42 

KG 

number  of  gust  orientations 

43 

NBOX 

number  of  aerodynamic  lifting  surface  boxes 

44 

NSBETO 

number  of  aerodynamic  slender  body  elements 

• 

45 

NE 

number  of  aerodynamic  bodies 

46 

NTGTAP 

nunber  of  total  aerodynamic  points 

48 

IDIKUL 

dimension  constant  for  unit  loads 

51 

AN 

load  factor 

i 

52 

7  DOT 

cLimb  rate 

53 

BTUPN 

turn  radius 

S  4 

KMAN 

manuever  constant 
=0,  no  manuever 

=1,  symmetric  pull-up  or  push-over 
=  2,  turn 

55 

A£ 

tank  angle 

56 

AC 

climb  angle 

57 

IKDSYM 

symmetry  inuicator 

61 

NOFMAX 

maximum,  number  of  orientations  to  be  considered 

62 

TTMEMX 

maximum  time  for  time  history,  sec. 

63 

EFB 

weapon  yield,  KT 

64 

KGFD 

control  constant  for  ground  reflection 

=0,  ro  ground  reflection 

=1,  include  ground  reflection 

f  s 

K  LPT 

=1,  iteration  for  critical  ranqe  desired 
-0,  no  iteration 

t 

f.  f 

HG6D 

height  of  ground  above  sea  level,  ft 

*  * 

K  LOAD 

=  1  ,  new  max  allowable  loads  input 

r 

NCHTS 

=1,  max  stresses  input 
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TABLE  15  (CONT'D) 

XTF  COMMON  BLOCK  DESCRIPTION 


Loc  Variable  Description 


70 

MXOFD 

maximum  potential  order  of  any  transfer  function 

71 

MXCFSN 

maxirrum  order  of  numerator  of  symmetric 
transfer  function 

72 

MXORSD 

maximum  order  of  denominator  of  symmetric 
transfer  function 

73 

MXOPAN 

maximum  order  of  numerator  of  antisymmetric 
transfer  function 

7  U 

MXORAD 

maximum  order  of  denominator  of  antisymmetric 
transfer  function 

75 

NTFS 

number  of  symmetric  transfer  functions 

76 

NTFA 

number  of  antisymmetric  transfer  functions 

80 

ISECT 

*0,  sectional  input  data  available 

81 

NDOF 

number  of  degrees  of  freedom  per  mass  station 

82 

IMS 

number  of  mass  property  items 
=1,  for  mass  point  data 
=  10,  for  NDOF=  6 
=  16,  for  NCOF=7 
=23,  for  NDOF=8 

91 

FINFPQ 

final  estimated  frequency  for  zero  forcing 
function  modulus 

92 

NFG 

number  of  first  gust 

93 

NLG 

number  of  last  gust 

96 

IPFNTM 

print  flag  for  inertial  module 

97 

CELT 

initial  delta  time  for  solutions 

99 

PBPADF 

scale  factor  for  rotation  modes 

100 

SIZFCT 

size  factor  defining  units  of  the  input  geometry 
data  relative  to  inches 
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TABLE  16 


AFOCOM  COMMON  ELOCK  DESCRIPTION 


Loc 

Variable 

Description 

1 

NTI 

unit  designator  for  fixed  data  deck  (=31) 

2 

MODES 

number  of  modes 

3 

NP 

total  number  of  panels  on  all  lifting  surfaces 

4 

MSTFIP 

maximum  number  of  strips  for  all  panels 
(for  dimensioning) 

5 

NSMAX 

max  number  of  strips  per  panel 

6 

NCMAX 

max  number  of  chordwise  boxes  per  panel 

7 

NTCTAL 

total  number  of  lifting  surface  boxes  4-  2MSBE 

8 

NB 

total  number  of  bodies 

9 

MSBE 

maximum  number  of  interference  body  elements 
(fcr  dimensioning) 

1C 

MBE 

maximum  number  of  interference  body  elements 
(for  dimensioning) 

11 

ND 

not  used 

12 

NE 

qrcund  effect  lag 

13 

NB  Y 

number  of  y-oriented  bodies 

14 

NBZ 

number  of  z-oriented  bodies 

15 

NTO 

=  NTP-fNT Y+NTZ  (see  below) 

16 

NIP 

total  number  of  lifting  surface  boxes 

17 

NTY 

number  of  y-oriented  interference  body  elements 

16 

NTZ 

number  of  z-oriented  interference  body  elements 

19 

NTYS 

number  of  y-oriented  slender  body  elements 

20 

NTZ  S 

number  of  z-oriented  slender  body  elements 

21 

MAXGF 

number  of  components  (groups  of  panels) 

22 

MAXSTF 

number  of  superstrips  on  all  panels 

23 

NSEETO 

total  number  of  slender  body  elements 

24 

NSTF.IP 

number  of  lifting  surface  strips 

25 

KF 

reduced  frequency 

26 

XM 

moment  axis 

27 

FEFA 

reference  area 

28 

FEFC 

reference  chord 

29 

REFS 

reference  seni-span 

30 

FMACH 

Mach  number 

31 

LINES 

naximum  number  of  print  lines  per  page 

TABLE  17 


ZZZ  COMMON  BLOCK  DESCRIPTION 


TABLE  18 


DISK2  COMMON  BLOCK  DESCRIPTION 


Lee  Variable  Description 


1 

ND2 

2-842 

ITBL2 

843 

NFECSA 

844 

NKD 

845 

NG 

846 

I  BUMP 

847 

VOBWS 

fortran  unit  number  for  mass  storage  data,  =2 
table  of  record  numbers  of  start  of  data  items 
total  number  of  records 
number  of  reduced  frequencies 
number  of  gust  orientations 
2  ( 1+NG) 

reduced  velocities 


TABLE  19 


DDTELS  COMMON  BLOCK  DESCRIPTION 


COK  MON  /DDTELS/DDTEL  (20,10) 


Loc 

V  ariafcle 

Description 

AEFOCYNAMIC  FILE  (TAELE  1) 

1,1 

NT 

unit  designator,  =19 

2,1 

TYPE 

file  type,  =4HAEPO 

3,1 

INOUT 

=1,  input;  =2,  input  and  output;  =3,  output 

4,1 

STAT 

current  status;  =1,  read;  =0,  write 

5,1 

LPN 

last  record  number  processed 

10,  1 

NK 

number  of  reduced  frequencies 

11,1 

NSYM 

number  symmetric  modes 

12,1 

NASYM 

number  of  antisymmetric  modes 

13,1 

NGUST 

number  of  gust  orientations 

14,1 

NE 

number  of  bodies 

15,1 

NEOX 

nunber  of  lifting  surface  boxes 

16,1 

N SEE TO 

number  of  slender  body  elements 

17,1 

NSTRIP 

number  of  strips 

18,  1 

MAXSTR 

number  of  superstrips 

19,1 

NP 

number  of  panels 

20,1 

KEF 

=0,  no  errors  processing  data 
=1,  error  occurred  reading 
=2,  error  occurred  writing 

UNIT 

LOAD  FILE 

(TAELE  4) 

1,“ 

NT 

unit  designator,  =34 

2,4 

TYPE 

file  type,  =4HUNIT 

3,4 

INOOT 

=1,  input;  =2,  input  and  output;  =3,  output 

4,4 

STAT 

current  status:  =1,  read;  =0,  write 

5.4 

LRN 

last  record  number  processed 

10,4 

N INTLD 

number  of  integrated  loads 

11,4 

NSTPSS 

number  of  stresses 

12,4 

NENGS 

number  of  engines  for  thrust  calculations 

13.4 

NSYM 

number  of  symmetric  modes 

14,4 

NASYM 

number  of  antisymmetric  modes 

15.4 

NK 

number  of  reduced  frequencies 

16,4 

NA3GFP 

number  of  aero  box  groups 

17,4 

NSBGRP 

number  of  slender  body  groups 

18,4 

NBOXES 

maximum  number  of  boxes  per  aero  group 

19,4 

NAERS3 

naximum  number  of  slender  body  elements  per 
slender  body  group 

20,4 

NER 

=0,  no  errors  processing  data 
=1,  error  occurred  reading 
=2,  error  occurred  writing 

175 


TABLE  19  (CONT'D) 

DFTBLS  COMMON  BLOCK  DESCRIPTION 


Loc 

V  ariable 

Description 

FREQUENCY  RESPONSE  FILE  (TABLE  5) 

1,5 

NT 

unit  designator,  =35 

2,5 

TYPE 

file  type,  =4HFRSP 

3,5 

INOUT 

=1,  input;  =2,  input  and  output;  =3, 

output 

4,5 

STAT 

current  status:  =1,  read;  =0, 

write 

5,5 

LRN 

last  record  number  processed 

11,5 

NGUST 

number  of  gust  orientations 

12,5 

NFREQ 

number  of  frequencies 

13,5 

NSYM 

number  symmetric  modes 

14,5 

NASYM 

number  of  antisymmetric  modes 

15,5 

NK 

number  of  reduced  frequencies 

used  in 

interpolations 

16,5 

NTOTAP 

number  of  total  aero  points 

17,5 

NFG 

number  of  first  gust 

20,5 

NEF 

=0,  no  errors  processing  data 

=1,  error  occurred  reading 

=2,  error  occurred  writing 

UNIT 

GUST  LOADS 

FILE  (TABLE  6) 

1.6 

NT 

unit  designator,  =36 

2.6 

TYPE 

file  type,  =4HLOAD 

3.6 

INOUT 

=1,  input;  =2,  input  and  output;  =3, 

output 

4,6 

STAT 

current  status:  =1,  read;  =0, 

write 

5.6 

LRN 

last  record  number  processed 

11,6 

NGUST 

number  of  gust  orientations 

12,6 

NFREQ 

number  of  frequencies 

13,6 

nintld 

number  of  integrated  loads 

14,6 

NACC 

number  of  accelerations 

20,6 

NEP 

=0,  no  errors  processing  data 

=1,  error  occurred  reading 
=  2,  error  occurred  writing 
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TABLE  20 


AERODYNAMIC  (AEFO)  FILE  DESCRIPTION 


Record 
vype  Word 

Item 

Description 

1  1 

TYPE 

data  set  ■ 

type  (=4HAERO) 

2 

NK 

number 

of 

reduced  frequencies 

3 

NS  YM 

number 

of 

symmetric  modes 

4 

NASYM 

number 

of 

antisymmetric  modes 

5 

NGUST 

number 

of 

gust  orientations 

6 

NB 

number 

of 

bodi  s 

7 

NEOX 

number 

of 

lift  ng  surface  boxes 

6 

NS  BF 

number 

of 

slen  or  body  elements 

9 

NSTFIP 

number 

of 

lifting  surface  strips 

10 

MAXSTF 

number 

of 

superstrips  on  all  panels 

1  1 

NP 

number 

of 

pane  s  on  all  lifting  surfaces 

12-20 

- 

not  used 

21 

FMACH 

mach  number 

22 

PEFA 

reference 

area 

23 

FEFS 

reference 

seni-span 

24 

FEFC 

reference 

chord 

25 

XM 

moment 

axis 

26-50 

- 

not  used 

2 


1-N 


DELA  box  areas  (DELA(I) ,1=1 ,  nbox) 

DELX  box  chords  (DELX (I) ,1=1 ,NBOX) 
xic  XAAS  of  centerline  of  1/4  chord  of  boxes 
(XIC(I)  ,1=1  ,NBOX) 

CG  cosine  of  dihedral  anqle  of  strips 

(CG  (I)  , I=1,NSTPIP) 

CS  chord  length  of  strips  (CS (I) , 1=1 ,NSTFIP) 

EE  half  width  of  scrips  (EE (I) ,1=1 ,NSTFIP) 

SG  sine  of  dihedral  anqle  of  strips 

(SG  (I)  , 1  =  1 , NSTRI P) 

YS  YAAS  of  centerline  of  strips 

( YS (I) ,1  =  1 ,NSTRIP) 

ZS  ZAAS  of  centerl  e  of  strips 

(ZS (I) ,1=1 ,NSTF  ) 

XIJ  XAAS  of  leading  dge  of  strip  centerline 

(XIJ(I)  ,1  =  1  , NST  P) 

ISSTF  superstrip  number  of  each  strip 
(ISSTF  (I)  ,1=1  ,NS1RIP) 

COOPD  spanwise  coordinate  of  strips 
(COOFD(I)  ,1  =  1  ,MAaSTP) 

NBAFAY  last  box  number  of  each  panel 
(NBAFAY (I) ,1=1 , N  ) 

NCARAY  number  of  chordwise  boxes  per  panel 
(NCAFAY (I) ,1=1 ,N  ) 

NSBEA  number  of  slender  body  elements  per  body 
(NSBEA(I) ,1=1 ,NB) 


TABLE  20  (CONT*  D) 


AERODYNAMIC  (AERO)  FILE  DESCRIPTION 


Record 

type 

word 

Iteir 

Description 

YE 

YAAS  of  body  centerlines  (YB (I)  ,1= 1 , KB) 

ZB 

ZAAS  of  body  centerlines  (ZB  (I)  ,  1=1  ,NB) 

3 

1-N 

XII 

XAAS  of  inboard  edge  of  1/4  chord  of  boxes 
(XII (I) , 1=1, NBO 

ETA1 

YAAS  of  inboard  edge  of  1/4  chord  of  boxes 
(ETAl  (I)  ,1=1  ,NEOX) 

ZETA1 

ZAAS  of  inboard  edge  of  1/4  chord  of  boxes 
(ZETA1  (I)  ,1=1,  OX) 

XI 2 

XAAS  of  outboard  edge  of  1/4  chord  of  boxes 
(XI 2  ( I)  ,  1=  1 , NBOX) 

ETA2 

YAAS  of  outboard  edge  of  1/4  chord  of  boxes 
(ETA2  (I)  ,1=1, NBOX) 

ZETA2 

ZAAS  of  outboard  edge  of  1/4  chcrd  of  boxes 
(ZETA2 (I) ,1=1, NBOX) 

4 

1-N 

X 

XAAS  of  centerline  of  3/4  chord  of  boxes 
(X(I) ,1=1, NBOX) 

ETA 

YAAS  of  centerline  of  3/4  chord  of  boxes 
(ETA (I) ,I=1,NBCX) 

ZETA 

ZAAS  of  centerline  of  3/4  chord  of  boxes 
(ZETA(I)  ,  1=1, NBOX) 

Record 

type  5 

is  omitted  if  NB=0 

5 

1-N 

XIS1 

XAAS  of  leading  edge  of  slender  body 
elements  (XIS1  (I)  ,1=1 ,NSBE) 

ETAS1 

YAAS  of  leading  edge  of  slender  body 
elements  (ETAS1 (I ) ,1=1 ,NSBE) 

ZETAS1 

ZAAS  of  leading  edge  of  slender  body 
elements  (ZETASl (I) ,1=1 ,NSBE) 

XIS2 

XAAS  of  trailing  edge  of  slender  body 
elements  (XIS2 (I) ,1=1 ,NSBE) 

ETAS  2 

YAAS  of  trailing  edge  of  slender  body 
elements  (ETAS2 (I) ,1=1 ,NSBE) 

ZETAS2 

ZAAS  of  trailing  -rdge  of  slender  body 
elements  (ZETAS2 (  ),I=1,NSBE) 

6 

1-N 

CP 

direction  cosines  of  gust  orientations 
((CP(I,J)  .1=1,3)  ,  =1,NG0ST) 

7 

1-NK 

RK 

reduced  frequencies 

178 


TABLE  20  (CONT*  D) 

AERODYNAMIC  (AEFO)  FILE  DESCRIPTION 


Record 

type 

Word 

Item 

Description 

6 

1-N30X 

HPS 

deflection  at  1/4 

chord 

of 

boxes 

due  to  symmetric 

modes 

Fepeated  for 

J  =  1  ,NS YM 

Record 

type  9 

is  emitted  if  NAS YM=0 

9 

1-NSOX 

HPA 

deflection  at  1/4 

chord 

of 

boxes 

due  to  antisymmetric  modes 
Fepeated  for  J=1#NASYM 


Record  types  10  thru  13  are  omitted  if  NB=0 


10  1-NSBE  HZ S  deflection  in  z-direction  at  midpoint  of 

slender  body  elements  due  to 
symmetric  modes 
Fepeated  for  J=1,NSYM 


Pecord  type  11  is  omitted  if  NASYM=0 


11  1-NSBE  HZA  deflection  in  z-direction  at  midpoint  of 

slender  body  elements  due  to 
antisymmetric  modes 
Fepeated  for  J=1,NASYM 


12  1-NSBE  HYS  deflection  in  y-direction  at  midpoint  of 

slender  body  elements  due  to 
symmetric  modes 
Fepeated  for  J=1#NSYM 


Record  type  13  is  omitted  if  NASYM=0 


13  1-NSBE  HYA  deflection  in  y-direction  at  midpoint  of 

slender  body  elements  due  to 
antisymmetric  modes 
Repeated  for  J=1,NASYM 


TABLE  20  (CONT'D) 


AERODYNAMIC  (AERO)  FILE  DESCRIPTION 


Record 

type  Word  Item  Description 


Record  types  14  thru  31  are  repeated  for  K=1,NK 
Record  types  14  thru  31  are  all  complex  type  arrays 


14  1-NBOX  DPOS  forces  on  liftinq  surface  boxes  due  to 

symmetric  modes 
Repeated  for  J=1,NSYM 


15  1-NBOX  FGPS  forces  on  lifting  surface  boxes  due  to 

symmetric  gusts 
Repeated  for  J=1,NGUST 


Record  types  16  thru  19  are  omitted  if  NB=0 

16  1-NSBE  DZOS  forces  on  z-oriented  slender  body 

elements  due  to  symmetric  modes 
Repeated  for  J=1,NSYM 


17  1-NSBE  FGZS  forces  on  z-oriented  slender  body 

elements  due  to  symmetric  gusts 
Repeated  for  J=1,NGUST 


16  1-N3EE  DYOS  forces  on  y-oriented  slender  body 

elements  due  to  symmetric  modes 
Fepeated  for  J=1,NSYM 


19  1-NSES  FGYS  forces  on  y-oriented  slender  body 

elements  due  tc  symmetric  gusts 
Fepeated  for  J=1,NGUST 


Record  types  20  thru  25  are  ommited  if  NAS YM=  0 


20  1-NBOX  DPOA  forces  on  lifting  surface  boxes  due  to 

antisymmetric  modes 
Fepeated  for  J=1,NASYM 


21  1-NBOX  FGFA  forces  on  lifting  surface  boxes  due  to 

antisymmetric  gusts 
Repeated  for  J=1,NGUST 


Record  types  22  thru  25  are  omitted  if  NB=0 


22  1-NSBE  DZOA  forces  on  z-oriented  slender  body 

elements  due  to  antisymmetric  modes 
Repeated  for  J=1,NASYM 
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TABLE  20  (CONT'D) 


AERODYNAMIC  (AERO)  FILE  DESCRIPTION 


Record 

type 

word  Item 

Description 

23 

1-NSBE  FGZA 

forces  on  z-criented  slender  body 
elements  due  to  antisymmetric  qusts 

Repeated  for 

J=1,NGUST 

24 

1-NSBE  DYOA 

forces  on  y-criented  slender  body 
elements  due  to  antisymmetric  modes 

Repeated  for 

J=1,NASYM 

25 

1-NSBE  FGYA 

forces  on  y-oriented  slender  body 
elements  due  to  antisymmetric  gusts 

Repeated  for 

J=  1 ,NGUST 

26 

1-N  DPOS 

generalized  forces  on  lifting  surfaces 
due  to  symmetric  modes 
(  (DPOS (I ,  J)  ,I=1,NSYM)  ,  J=1 ,NSYM) 

FGPS 

generalized  forces  on  lifting  surfaces 
due  to  symmetric  gusts 
((FGPS (I, J) ,I=1,NSYM) ,J=1,NGUST) 

Record 

types  27  and 

28  are  omitted  if  NE=0 

27 

1-N  DZOS 

generalized  forces  on  z-oriented  slender 
bodies  due  to  symmetric  gusts 
(  (DZOS  (I,J)  ,1=1 rNS¥M)  ,  J=1 ,NSYM) 

FGZS 

generalized  forces  on  z-oriented  slender 
bodies  due  to  symmetric  gusts 
(  ( FGZS ( I, J)  , 1=  1  ,NSYM'  , J= 1 , NGUST) 

28 

1-N  DYOS 

generalized  forces  on  y-oriented  slender 
bodies  due  to  symmetric  modes 
((DYOS(I,J) ,I=1,NSYM  , J=1 , NSYK) 

FGYS 

generalized  forces  on  y-oriented  slender 
bodies  due  to  symmetric  gusts 
( (FGYS (I,J) , 1= 1,NSYM)  , J=1 ,NGUST) 

Record 

types  29  thru  31  are  omitted  if  NASYM=0 

29 

1-N  DPOA 

generalized  forces  on  lifting  surfaces 
due  to  antisymmetric  modes 
( (DPOA (I , J)  ,1=1, NASYM)  ,J=1,NASYM) 

FGPA 

generalized  forces  on  lifting  surfaces 

due  to  anti  symmetric  g usts 
(  (FGPA ( I, J) ,  T= 1 , NAS YM)  ,J=1,NGUST) 


TABLE  20  (CONT'D) 

AERODYNAMIC  (AERO)  FILE  DESCRIPTION 


Record 

type 

Word 

Item 

Description 

Record 

types 

30  and 

31  are  oiritted  if  NB=0 

30 

1-N 

DZOA 

FGZA 

generalized  forces  on  z-oriented  slender 
bodies  due  to  antisymmetric  modes 
( (DZOA (I ,  J)  *1  =  1, NAS YM)  ,J=1,NASYM) 
generalized  forces  on  z-oriented  slender 
bodies  due  to  antisymmetric  gusts 
(<F G2P (I,J) ,1=1, NASYM) ,J=1,NGUST) 

31 

1-N 

DYOA 

FGY? 

generalized  forces  on  y-oriented  slender 
bodies  due  to  antisym  etric  modes 
( (D YOA(I , J)  ,1=1, NAS  YM  ,J=1,NASYM) 
generalized  forces  on  y-oriented  slender 

bodies  due  to  antisymmetric  gusts 
(  (FGY£  (  I,  J)  ,  !•=  1,NASYM)  ,J=1,NGUST) 
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TABLE  21 


UNIT  LOAD  (UNIT)  FILE  DESCRIPTION 


Record 


type 

Word 

Iteir 

Description 

1 

1 

TYPE 

data  set  type  (=4HUNIT) 

2 

IDENT 

identification  number 

3 

NINTLD 

number  of  integrated  loads 

a 

NSTPSS 

number  of  stresses 

5 

NSNGS 

number  of  engines  for  thrust  calculation 

6 

NS  YM 

number  of  symmetric  modes 

7 

NASYM 

number  of  antisymmetric  modes 

8 

NK 

number  of  reduced  frequencies 

9 

NABGRF 

number  of  aero  box  groups 

10 

NSBGRF 

number  of  slender  body  groups 

11 

NBOXES 

maximum  no.  of  boxes  per  aero  box  group 

12 

NAERSE 

maximum  no.  of  slender  body  elements  per 
slender  body  group 

2 

1-N 

STALDS 

Integrated  load  definitions 
((STALDS (I, J) ,1=1, NINTLD) ,J=1,8) 

3 

1-N 

STRESS 

Stress  definition  matrix 
( (STFESS  (I,J) ,I=1,NSTRSS) ,J=1 ,  NINTLD) 

u 

1-N 

THRLOC 

Integrated  loads  due  to  thrust 
(  (THRLOD  (I,  J)  ,1=1  .NINTLD)  ,J  =  1,NENGS) 

5 

1-N 

THPGNF 

Generalized  loads  due  to  thrust 
( (THFGNF (I, J) ,1=1 ,NSYM) ,J=1,NENGS) 

6 

1-N 

PIO 

Integrated  inertial  loads  due  to  unit  modal 
amplitudes 

(  (PIQ  (I  ,  J)  ,1  =  1, NINTLD)  ,J=1,NM) 
where  NM=NS YM+NAS YM 

record 

types 

7  and  8  are  repeated  for  K=1,NK 

7 

1-N 

PPQS 

Integrated  aero  loads  due  to  unit  medal 
amplitudes  for  symmetric  modes 
( (PADS (I,J) , 1= 1, NINTLD)  , J  =  1,NSYM) 

8 

1-N 

PfiVA 

Integrated  aero  loads  due  to  unit  modal 
amplitudes  for  antisymmetric  modes 
((PAQA(I,J) ,1=1, NINTLD) ,J=1, NASYM) 

9 

1-N 

SYMCOD 

Symmetric  centerline  load  modifier 

(SYMCOD  (  I)  ,  1=1 .NINTLD) 
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TABLE  21  (CONT*  D) 

UNIT  LOAD  (UNIT)  FILE  DESCFIPTION 


Record 

type 

Word 

Iteir 

Description 

1C 

1-N 

ASNCOD 

Antisymmetric  centerline  load  modifier 
(ASMCOD(I)  ,  1=  1  ,NINTLD) 

11 

1-N 

NFNLAE 

Aero  box  group  definitions 
( (NFNLAb (I, J) ,1=1,3) ,J  =  1,NABGPP) 

12 

1-N 

PINTP 

Integrated  load  due  to  aero  boxes 
(( (PINTP (I,J,K) ,1=1 , NINTLD) ,J=1,NABGFP) , 

K=  1 ,  NEOXES) 

13 

1-N 

NFNLSE 

Slender  body  element  group  definitions 
( (NFNLSB (I , J) ,1=1 ,3) ,J=1 ,NSBGFP) 

1U 

1-N 

PINTZ 

Integrated  loads  due  to  slender  body 

Z-force 

(( (PINTZ(I,J,K) ,1=1, NINTLD) ,J=1,NSBGRP) , 

K=  1  ,  NAEFSE) 

15 

1-N 

PINTY 

Integrated  loads  due  to  slender  body 

Y- force 

( ( (PINTY (I, J, K) ,1=1, NINTLD) ,J=1,NSBGFP) , 

K=  1  ,NAEFSB) 
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TAELE  22 


FPEQUENCY  F ESPONSE  (FKSP)  FILE  DESCFIPTION 


Peccrd 

type 

Word 

Item 

Description 

1 

1 

TYPE 

data  set  type (=4HFFSP) 

2 

I  LENT 

identification  number 

3 

NGUST 

number  of  qust  orientations 

4 

NFFEQ 

number  of  frequencies 

5 

NS  YM 

number  of  symmetric  modes 

6 

NASYM 

number  of  antisymmetric  modes 

7 

NK 

number  of  reduced  frequencies  used  in  the 
interpolations 

8 

9-10 

NTOTAP 

number  of  total  aero  points 
not  used 

11 

VFL 

true  velocity  (ft/sec) 

12 

SIGMA 

density  ratio 

2 

1-N 

CF 

gust  orientation  direction  cosines 
{(CF(I,J)  ,1  =  1,3)  ,J=1, NGUST) 

3 

1-N 

OMEGA 

frequencies  (rad/sec) 

(OMEGA  (I)  ,1=1,  NFFEQ) 

Feccrd 

types 

4  thru  9 

are  repeated  for  K=1 , NFFEQ  and  N=1, NGUST 

4 

1 

NOB 

orientation  number 

2 

IFFEQ 

frequency  number 

3 

FF  EQ 

frequency  (HZ) 

4 

IZFOS 

orientation  for  zero  symmetric  qust  loads 

5 

IZFOA 

orientation  for  zero  antisymmetric  gust 
loads 

5 

1-N 

COEF 

interpolation  coefficients 
(COEF  (I )  ,1  =  1  ,NK) 

6 

1-N 

IVBHI 

reduced  frequency  numbers  of  aero  used  for 
interpolation  of  this  frequency 
(IVBWI(I) ,1=1 ,NK) 

7 

1-N 

Q 

complex  generalized  displacement 
(Q (I) ,1=1 ,NM) ,  Where  NM=NSYM+NASYM 

8 

1-N 

FGAFOS 

complex  symmetric  aerodynamic 
element  gust  forces 
(F  ( I)  ,  I  =NTOTAP) 

9 

1-N 

FGAFOA 

complex  antisymmetric  aerodynamic  element 

gust  forces 

(F ( I)  ,I=1,NTOTAP) 
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TABLE  23 

UNIT  GUST  LOADS  (LOAD)  FILE  DESCRIPTION 


Record 


type 

Word 

Iteir 

Description 

1 

1 

2 

3 

4 

5 

6 

7-10 

11 

12 

13-50 

TYPE 

I  DENT 

NGUST 

NFPEQ 

NINTLE 

NACC 

VEL 

SIGMA 

data  set  type  (=4HLCAD) 

identification  number 

number  of  gust  orientations 

number  of  frequencies 

number  of  integrated  loads 

number  of  acceleration  mass  points 

not  used 

true  velocity  (FPS) 
density  ratio 
not  used 

4. 

1-N 

STALDS 

Integrated  load  definitions 
((STALDS (I, J) , I=1,NINTLD) ,J=1, 8) 

3 

1-N 

CP 

gust  orientation  direction  cosines 
((CR(I,J)  ,1=1,  3)  ,J=1,NGUST) 

4 

1-N 

OK  EGA 

freguencies  (rad/sec) 

(OMEGA (I) , 1= 1r  NFREQ) 

Fecord 

type  5 

is  oiritted  if  NACC=0 

5 

1-N 

INDACC 

acceleration  mass  point  and  degree  of 

freedom  numbers 

((INDACC (I, J) ,1=2,) ,J=1,NACC) 

F ecord 
Fecord 

types  6 
types  6 

thru  9  are  repeated  for  K=1,NGUST 
and  7  are  omitted  if  NACC=0 

6 

1-N 

AS 

complex  symmetric  accelerations 
((AS(I,J) ,1=1, NACC) ,J=1,NFFEQ) 

7 

1-N 

AA 

complex  antisymmetric  accelerations 
( (AA  ( I,  J)  ,1=1,  NACC)  ,  J=  1 ,  NFF  EQ) 

8 

1-N 

PS 

complex  symmetric  integrated  loads 
((PS  (I,J)  ,I=1,NINTLD)  ,J=1,NFPEQ) 

9 

1-N 

PA 

complex  antisymmetric  integrated  loads 

(  (PA (I* J) , 1= 1,NINTLD) ,  J=1 ,  NFPEQ) 
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TABLE  24 


MASS 

STORAGE  FILE  DESCRIPTION 

Record 

type 

Word 

Item 

Description 

Record 

types 

1  thru  4 

are  repeated  for  K=1 ,NK 

1 

1-N 

DS 

complex  generalized  forces  doe  to  symmetric 
modes,  ( (DS  (I,J) ,1=1, NSYM) ,J  =  1 ,NSYM) 

2 

1-N 

DA 

complex  generalized  forces  due  to 
antisymmetric  modes 
( (DA(I,J) ,1=1, NASYM) ,J=1 , NASYM) 

Record 

types 

3  and  4  are  repeated  for  N=1,NG 

J 

1-N 

FS 

complex  gust  forces  due  to  symmetric  modes 
(FS (I) ,1= 1,NT0TAP) ,  where  NTOTAP=NBOX+2NSBE 

4 

1-N 

FA 

complex  gust  forces  due  to  antisymmetric 
modes,  (FA (I) , I=1,NTOTAP) 

5 

1-N 

SPLHS 

aero  force  integration  matrix  for  symmetric 
modes 

((SPLHS (I, J) ,1=1 ,NTOTAP) ,J=1,NSYM) 

fc 

1-N 

SPLHA 

aero  force  integration  matrix  for 
antisymmetric  modes 
((SPLHA (I, J) ,1=1, NTOTAP) ,J=1, NASYM) 

7 

1-N 

GEOMBX 

AAS  coordinates  of  inboard  and  outboard 
edge  of  1/4  chord  of  aero  boxes 
( (GEOMBX (I , J) ,I=1,NB0X) ,J=1,6) 

8 

1-N 

GEOMBC 

AAS  coordinates  of  leading  and  trailing 
edge  of  slender  body  elements 
( (GEOMBD (1 , J) ,1=1 ,NSBE) ,J=1,6) 

NSYM  *  number  of  symmetric  modes 
NASYM  =  number  of  antisymmetric  modes 
NEOX  =  number  of  lifting  surface  boxes 

NSBE  =  number  of  slender  body  elements 


TABLE  25 


SEGLOAD  CIFECTIVE  LISTING 


A 1  TPEE  GEOM 

A2  TFEE  SPLINE 

SPLINE  INCLUDE  SOFT 

A?  TFEE  DOTP 

DOTP  INCLUDE  OPGN 

AU  TFEE  GEND 

SEND  INCLUDE  DPPS, CPZ ¥, D YPZ, DZPY , SUBB , SUEP 
A 5  TPEE  SP 

SP  INCLUDE  DOMULT, DZ YMAT , FO WD YZ ,MUZ YC 

A 6  TPEE  BFSMAT 

BFSMAT  INCLUDE  FWMIJ,FZY2 

A 7  TPEE  WANCWT 

WANCVT  INCLUDE  GO  ST, FHSIDE, SOLVIT, MATMUL 
A 8  TPEE  GENF 

A 9  TFEE  AEFO 

A10  TPEE  NEWH 

All  TREE  PISTON 


Ml  TPEE  CSDLM-(A1,A2,A3,A4#A5, A6, A7,A8,A9,A10, A1 1) 

CSDLM  INCLUDE  SDLM 
CSDIK  GLOBAL  KDS ,  ELM 
M2  TFEE  INEFTN 

M3  TpEE  ACS 

MU  TREE  LOAD**  (BEAM, GEOMME # MSSPHI  , GEOMABrGEOMSB, AFOLOD) 

LOAD  INCLUDE  TFAEPC 

M5  TPEE  CFPEQF 

M6  TFEE  CFFLOL 

M7  TFEE  CGUST- (CTPIM,GUSTDR) 

GUSTEF  INCLUDE  GSTHST 
M8  TREE  CPIGIC 

M9  TPEE  CMEFGE 

MAIN  TPEE  EACGUST- (Ml ,M2,M3,M4,M5,M6,M7,M8,M9) 

GLOBAL  DISK2 ,ZZZ,XTF,AEFOMX,DCTELS 
END 


SECTION  VII 


EXAMPLE  PROBLEM 

The  sample  problem  chosen  for  analyses  demonstration  consists  of  a 
large  twin  engine  transport.  The  number  of  degrees  of  freedom  for  the  model 
are  sufficient  for  demonstration  of  the  large  analyses  capability  of  the 
program,  though  fewer  aero  strips,  mass  points,  integrated  loads,  and 
elastic  modes  were  used  than  necessary  for  a  complete  analyses  of  an  air¬ 
craft  of  this  size;  116  boxes,  19  slender  body  elements  and  28  masses. 

Figure  23  shows  the  aero  node  points  and  aero  strips.  Figure  24 
shows  the  mass  points  and  beam  network  used.  Figure  25  defines  the  inte¬ 
grated  loads.  The  modes  consisted  of  nine  rigid  body  modes  and  trim  modes 
plus  3  symmetric  and  4  antisymmetric  elastic  modes.  The  elastic  modes  are 
simple  symmetric  and  antisymmetric  surface  flapping  modes  of  the  wing  and 
horizontal  and  vertical  tail  surfaces  and  simple  in-plane  and  out-of-plane 
fuselage  hinge  modes. 

Details  of  the  model  data  may  be  found  in  the  sample  data  input 
listings  and  sample  run  data  which  are  given  in  Tables  26  and  27  through  30. 

Table  27  consists  of  the  printout  of  a  pass  through  the  aerodynamic 
module  for  a  reduced  frequency  of  zero.  Table  28  illustrates  the  printout 
of  the  unit  load  module.  Table  29  is  for  the  Frequency  Response  and  Unit 
Gust  Load  Modules.  Only  selected  orientation  output  has  been  included  for 
illustrative  purposes.  Table  30  details  the  output  for  a  single  iteration 
pass  through  the  trim  and  blast  portion  of  the  program  for  one  orientation. 
The  plotted  data  are  shown  for  the  first  iteration  only. 
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EXAMPLE  PROBLEM  INPUT  DATA  LISTING 
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TABLE  2a 


EXAMPLE  PROBLEM  OUTPUT  LISTING 
UNIT  LOAD  MODULE 
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TABLE  29 

EXAMPLE  PROBLEM  OUTPUT  LISTING 
FREQUENCY  RESPONSE  AND  UNIT  GUST  MODULES 
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